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BINDING 
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PRO) GEEDINGS 


Binding cases and cloth for Volume 60 
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It is suggested that in general Volume 60 
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ber) should be bound together to form one 
volume for 1948. This is particularly advis- 
able in view of the extreme dearth of the 
special green binding cloth used for the 
Proceedings. 
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At last! * Electronics 


—a really 


By F. G. Spreadbury, ~ 
A.M. Inst.B.E., 


guide to the Lecturer in Physics and Mathematics at the 
Working Men’s College, London. 


authoritative 


subject 


This is the most up-to-date, comprehensive and 
reliable guide yet published to electronic theory and 
Its applications. It is invaluable to all workers in 
this field, including radio and electrical engineers 
and technicians, and all students of modern physics, 


700 pages. Illustrated. 55/- net. 


Pitman 
Parker Street, Kingsway, London, W.C.2 
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GO-AX CABLES 


ATTEN 
db/100#. 


LOADING 
AW 


Bending Radhius 5° 
FPhotocell Cable, 


TRANSRADIO LTD. 138, CROMWELL ROAD, LONDON, S.W.7. 
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one step ahead! 


Our experience in the manufacture of toroidal- 


wound resistors goes back such a long way 
that we are confident of always maintaining 
the lead in design and efficiency. As the pioneers of this type of unit 
in this country, we can offer our customers a product which is the 
result of long research and practical experience, and we are also able to 
supply more alternative arrangements than any of our competitors. 


TOROIDAL WOUND POTENTIOMETERS 


BRITISH ELECTRIC RESISTANCE CO. LTD., Queensway, Ponders End, Middlesex 
BR4013—EHI Telephone: Howard 1492. Telegrams: Vitrohm, Enfield 


HIGH VAGUUM GAUGES NOW AVAILABLE 


“5G THE B.A.R. THERMOCOUPLE-IONISATION GAUGE 

CONTROL. The newly developed Thermocouple- 
lonisation Gauge Control, complete with Thermocouple and 
lonisation Gauges, covers the pressure range from 2x 10-7 mm. to 
Imm. Hg operation is dependable and simple. Ranges are 10—3— 
Imm. and 0-5, 0-1.0, 0-0.1, 0-0.01 microns. The unit is removable 
from its cabinet for incorporation in a central panel. 


THE B.A.R. ALPHATRON is an_ ionisation type Se 
vacuum gauge using the ionising power of alpha 
particles from a radium source to measure total pressure 
of any gas, vapour or mixed atmosphere from Ip to 10 mm. 
Hg with instantaneous linear response. The B.A.R. Alphatron 
is quickly available from batches now in production. 


DESIGNERS AND MANUFACTURERS OF HIGH VACUUM GAUGES VALVES SEALS” 
OIFFUSION PUMPS STILLS FURNACES COATING EQUIPMENT AND DEHYDRATION PLANT | 


BLOCK €E2 - HILLINGTON NORTH GLASGOW ~< SW Hy 
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etd. . 


DEVELOPMENTS 


iOkV DC 


from two rectifiers, type 
SOEI 45 anaes 2502 C.p:s: 
voltage-doubler circuit. 


bkV DC 


from three rectifiers, 
S6eimio5i,. ina pulse 
voltage tripler circuit. 


5kV DC 


from 350 volts A.C. using a 
Westeht EHT unit. 


type 


Delivery ex stock. 


METAL RECTIFIERS 


Write for literature to Dept. P.S. 5. 


Interested manufacturers may obtain 
small supplies as samples. 


* Each only 7” Xx 7}” 
+ Each only 7” Xx 23” 
Nestinghouse Brake & Signal Co., Ltd. 
2 York Way, King’s Cross, London,N.! 


Sind 


} 


For 


MEN ONLY 


interested in 


GLAND SEALS REFRIGATION 
CONTROL : THERMOSTATIC and 
PRESSURE OPERATED DEVICES, etc. 


Seamless one-piece metal bellows . 
produced by a process unique in this country 
. . no spinning, no annealing, no localised 
strain or thinning. Every bellows uniform in 
life and performance, absolutely reliable in 
operation, pretested during forming. 


Root diam. 3” to 3”. Outside diam. {” to 43” 


Send for List No. V800-1 


Drayton Hydroflex 


METAL BELLOWS 


DRAYTON REGULATOR & INSTRUMENT Co. Ltd. 
WEST DRAYTON West Drayton 2611 MIDDX. 
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but there's nothing 


more attractive than 


“TECONAL"” PERMANENT 


REGD. TRADE MARK 


MAGNETS wn: ov Mara 


(co) 0 MULLARD ELECTRONIC PRODUCTS LIMITED, MAGNET DIVISLON, 


CENTURY HOUSE, SHAFTESBURY AVENUE, LONDON, W.C.2. 
(iT234)) 


THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


Vii 


Quartz crystals enable astronomers to detect small 
changes in the earth’s rate of rotation immediately they 
occur. Numerous observations on the moon and stars 
were previously necessary to detect such changes, 


GREENWICH OBSERVATORY uses the quartz crystal clock to help 
maintain the nation’s standard of time. The earth itself, which rotates on 


ie its axis at a remarkably constarit rate, is our fundamental standard of 
BS pal time. The modern crystal oscillator, however, has a frequency stability 


better than the earth itself, and has enabled small irregular changes in the earth’s rate of 
motion to be detected. 


For generation of stable frequency there is nothing better than the quartz crystal 
_ oscillator. 


QUARTZ CRYSTAL UNITS 


Om CROs Maeve Un Nes inGs Ami linOm Nips ERORURIS re MalcaaNe hi 


ALFORD Fale Cc) t, Ro 1 C cA LE INSTRUMENT S Bay Phe 
=EL WORKS, SALFORD 3 Phone: 6688 (6 lines) Grams and Cables: Sparkless, Manchester 


tcoprietors: THE GENERAL ELECTRIC CO. LTD. of England 
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There is a G.E.C. electronic device to me 
mercial, industrial or scientific need. 


receiving types illustrated ab 
U52, KT66 ana PXz5. 


in very short supply, 
on request. 


et every com- 
The particular 
ove are Z14, X61M, KT6l, 
Valves are still 
dly be supplied 


Osram r5.C Osram 
PHOTO CELLS CATHODE Raya VALVES 
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VITREOSIL MERCURY 
VAPOUR PUMPS | PHOTOELECTRIC EQUIPMENT 


| THE DENSITOMETER 


THE EEL P.R.U. Densitometer, 

designed in conjunction with 
the Pneumokoniosis Research 
Unit of the Medical Research 
Council, is a simple, compact, 
direct reading photoelectric 
instrument, to measure the 
optical density of stains on 
filter paper. 


This instrument may be used 
for many other purposes where 
it is desired to obtain a meter 
reading indicating the com- 
parative densities of materials. 


PRICE 24 Guineas 


with supply of sample clips 
(excluding 6-volt accumulator). 


This new M.V. Fore Pump will operate from an i A.C. MODEL, incorporating 
ordinary water filter pump, and when used in con- i constant voltage transformer. 
junction with our Single-Stage or Two-Stage Pump, . 

pressures less than 0°00002mm Hg are attained. UG SOR 


Write for full particulars of this and other 


Write for descriptive leaffet EEL photoelectric equipment incorporating the 
famous EEL selenium cell. 
THE THERMAL SYNDICATE LTD. A product of 
Head Office: YVallsend, Northumberland EVANS ELECTROSELENIUM 
London Office: 12-14 Old Pye Street, S.W.| Harlow LTD. Essex 


FRELIIT 189 Pore 
3} 
Wl ditiies 


PIPE COUPLINGS 


—__> 


ELECTRICALLY HEATED PRESSURE HEADS 


=> 


FILM ASSESSORS AND SCANNING MICROSCOPES 


—_ 


ANN 


CONTINUOUS FILM RECORDING CAMERAS AND EQUIP- 
MENT FOR CATHODE RAY OSCELLOGRAPHY, ETC. 


HAL ATRIA 


We undertake the Design, Development 
and Manufacture of any type of Optical 


— Mechanical — Electrical Instrument. 
Including Cameras for special purposes. 


Avimo Limited, Taunton, England + Telephone Taunton 3634 


117909 
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NICKEL-CHROMIUM 80/20 
Used almost without exception for all high value wire-wound fixed 
and variable resistors. It possesses a very high resistance to corrosion, 
is non-magnetic and combines a high resistivity with a low tempera- 


ture coefficient. 


COPPER-NICKEL [6/44 
Well known as Constantan or Ferry, is characterised by a moderately 
high resistivity together witha very low temperature coefficient 
and is widely used in measuring instruments. 


MINALPHA 
A manganese-nickel-copper alloy superior to the older Manganin 
alloy in respect of both temperature coefficient and thermo-electric 
effect. Employed in standard resistances and for coils in measuring 
apparatus where resistance must remain constant despite fluctuating 


temperatures. 


MANCOLOY 10 
A copper alloy with a low resistance combined with a relatively low 
temperature coefficient. Its use is advantageous when a low but 
practically constant resistance is required in instruments and radio 


apparatus. 


Comprehensive information on these and other materials, 
together with tables of resistance per yard, and tolerances, 
are given in an illustrated booklet, ‘‘ Electrical Resistance 
Materials’, available as publication 1440. 


vonmane Johnson (Mp 
Matthey 


JOHNSON, MATTHEY & CO., LIMITED, HATTON GARDEN, LONDON, E.C. 


GD53 
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H iL PS ELECTRON MICROSCOPE 


This instrument —-entirely new in design — has a performance far in 
° advance of any Electron Microscope of the conventional type. Some of 
the important features include : 


ENTIRELY NEW FOCUSSING DEVICE - HIGH RESOLVING POWER - TWO EXTRA 


LENSES +» CONTINUOUSLY VARIABLE MAGNIFICATION PRACTICALLY 
UNINTERRUPTED OPERATION : LARGE SIZE IMAGE OF 20 CM. AT ALL 


MAGNIFICATIONS - ENTIRELY SHOCKPROOF AND RAYPROOF 


PHILIPS ELECTRICAL 


EU i a ie 


POLES ER LAL XA DIE P Ty, ;CEN TUR Y HOUSE, SHAFTESBURY AVENUE, LONDON; er 
XB3I5A)» 
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There’s. more in resistance work 
than meets the eye! : 


| 

Muy have remarked on the fine appearance 
of Muirhead Decade Resistance Boxes — 
but we are even more concerned with what 
we put inside them. 

Take, for instance, the A-25 series; all 
resistors are wound on cards by methods that 
ensure very low time-constants. The in- 
ternally mounted switches are protected from 
contamination and have low and constant 


contact resistance. 
Write for Publication C-102-A which gives full details and specification 


Sev. Se | 


Muirhead & Co. Limite, Elmers End, Becken tam, Kent. Tel. Beckenham 004] -2 
FOR OVER 60 YEARS DESIGNERS & MAKERS OF PRECISION INSTRUMEN' 


See i 


CRC: 
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COSSO 


The Model 1035 is a_ general 
purpose Oscillograph, consisting 
of a Double Beam Tube Unit, 
Time Base, Y Deflection Ampli- 
fiers and internal Power Supplies. 
The two traces are presented over 
the full area of a flat screen tube 
of 90 mm. internal diameter and 
Operating at 2 kv. Signals are 
normally fed via the Amplifiers, 
with provision for input voltage 
calibration. The Time Base is 
designed for repetitive, triggered, 
or single stroke operation, and 
time measurement is provided by 
a directly calibrated Shift Control. 


| MODEL 1049 
Industrial Oscillograph 
is designed specifically for indus- 
trial use where the main interest 
js in the observation and measure- 


ment of low frequency phenomena. 


Its presentation is generally similar 
to that of Model 1035 illustrated 


and a comprehensive specification For photographic recording, both 

| models have provision es bee 
i i attachment of a Camera, Model 
includes 4 ky tube Operation for 1428, which may be operated 
i H either manually or by motcr drive. 
transient recurding. A similar camera, Medel 427, 
is also available for use with 


Model 339 Oscillograph. 


Further details on application to :— 
A. C. COSSOR LED., INSTRUMENT DEPT., HIGHBURY, LONDON, N.S 


ci 3 


en ek, 
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F. A Rev Seas 


ELE@ ERO MEG ra 
FOR RADIOLOGICAL WORK 


A unique electronic instrument for 
research and routine testing in 
Hospital Radium and X-ray Therapy 
Departments. 


It has an input resistance of 101° ohms 
and an input capacity of less than 
1 ppl. Developed primarily for use 
in Radiological work, where small 
condensers of the Sievert type are 
used extensively for the measurement 
of gamma and X-ray intensities. 


Fully descriptive leaflet supplied on 
request. 


a BROOKLANDS WORKS, PRINCES ROAD, DARTFORD, KENT. 
oes Telephone: DARTFORD 2989 


BALDWIN INS 
VOLTAGE RANGES | ALDWIN INSTRUMENT COMPANY LTD 
1H11| 


| 
IM MICROID FLASK SHAKER 


A new variable speed “ wrist-action”’ shaker—silent, portable_and efficient 


This - handy — shaking 
machine will stand on 
the laboratory bench. 
Total load, four 500 ml. 
flasks each half full. 
Frequency of shaking 
continuously var.able 
up to 500 cycles/min. 
Power consumption 60 | 


watts. The Microid = : : L 
Flask Shaker is silent in : AES LE aie 


WZ, . . 
operation. It abolishes free on application 
the tedium of hand 
shaking. It can find a use in every chemical laboratory. 


Grirrin and JATLO oKLtd 


Established as Scientific Instrument Makers in 1826 


LONDON MANCHESTER GLASGOW E 
D 
Kemble St., W.C.2. 19 Cheetham Hill Rd., 4. 45 Renfrew St., C.2, 7 Tene Rae 


BIRMINGHAM : STANDLEY BELCHER & MASON LTD., Church Street, 3 


ooo 
i 
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‘Kodak’ 0.250 Metallographic 
Plate has a fine-grain high- 
contrast emulsion of high re- 
. solving power, especially colour- 
the plate for PHOTOMICROGRAPHY -ensiizeforphotomicrography 
and metallography. It has its 
highest sensitivity in the region 
of maximum visual sensitivity, 
where most microscopes yield 
sharpest focus, and compara- 
tively little speed is therefore 
lost when working with a green 
or yellow filter. 
0.250 is one of a long line of 
special-purpose materials pro- 
duced by Kodak for applied 


photography. For details write 
to: 


KODAK 


* DEPT. 540 - KODAK HOUSE - KINGSWAY - LONDON - W.C.2 


@ Cast commercially pure 
lead, electrolytically polished 
and etched. X25. 


The same, X1000. @ 
Photomicrographs made on 
‘Kodak’ 0.250 Plates. (Courtesy 
of Ford Motor Co. Ltd.). 


- KODAK LIMITED 


‘ WE , MEASURING INSTRUMENTS 


Decade Resistance Boxes 
201 Logarithmic Resistor 

202 Decade Voltage Divider 
210 Decade Condenser Unit 


VOUIMETER 230 Decade Inductometer 
Type 613A 300 Ratio Arm Unit 

erat ce 301 Limit Resistance Bridge 
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by baking for four hours. The resistance of the thermometer coil was designed 
to suit the Wheatstone bridge and galvanometer available and enabled the tem- 
perature to be measured to +0-02° below 20°K. The 308.w.c. Eureka leads 
passed through a heat trap connected to the helium expansion chamber to reduce 
heat conduction to the calorimeter. 

The second calorimeter was 3-0cm. in diameter and 3-8cm. long to accom- 
modate the larger specimens of mercury and sodium. The heater coil of 
40s.w.c. Eureka was of 180 Q resistance and the 47s.w.c. Eureka thermometer 
coil was of 485Q resistance. As measurements with the mercury were 
also required at temperatures well above 20°K., where the Eureka thermo- 
meter becomes insensitive, it was supplemented by a platinum resistance thermo- 
meter of 40 « bare wire wound on a thin paper base and attached to the calorimeter 
by bakelite varnish. This coil had a resistance of 322Q at 273°K., falling to 
4-5Q at 25°x. Arrangements were made to measure the lead resistance to the 
platinum coil. The temperature corresponding to the platinum resistance 
was determined from tables published by Clusius and Vaughen (1929). The 
wire used inthe present experiments followed closely the values for the sample 
Pt-23’-1915 quoted in those tables. The Eureka thermometers were calibrated 
in the hydrogen and helium regions by means of the vapour pressure/temperature 
data for these liquids, and in the intermediate region an interpolation was made 
on the basis of the information obtained in a previous experiment (Pickard and 
Simon 1948). 

The samples were cooled to about 10°. by condensing hydrogen into the 
calorimeter and then pumping it away. ‘They were further cooled by condensing 
helium into the calorimeter and then pumping it. When all this helium had been 
removed, a small quantity of helium gas was re-introduced to permit thermal 
exchange to take place between the calorimeter and the sample. When the 
temperature drift of the calorimeter had become steady the calorimetric 
measurements were begun. ‘The temperature was measured at 20-sec. intervals 
throughout this stage of the experiments, heat being supplied for periods of one 
or two minutes with intervals to allow the temperature drift to be measured. 
‘The experiments were carried out in such a way that the temperature difference 
between the calorimeter and its enclosure was always small and the thermal 
isolation of the calorimeter was good. It was therefore unnecessary to use any 
special method to correct for imperfect heat insulation. The rise in temperature 
due to the heating was determined by plotting a temperature/time curve before 
and after the heating and extrapolating to the middle of the heating period, 
thus correcting for the temperature drift. 

During the heating the temperature of the calorimeter rose above that of 
its contents, but thermal equilibrium was rapidly attained after switching off 
the heater current. It was found that the use of different heater currents intro- 
duced no systematic error into the results. The “overheating” was of the order 
of $° with the mercury sample; with the palladium it was less than this, and with 
the sodium it was so small as to be hardly noticeable. 

Correction was made for the heat capacity of the calorimeter using the values 
for the specific heat of copper and tin measured at Leiden (Kok and Keesom 
1936, Keesom and Ende 1932). The heat capacity of the resistance coils was 


negligible and that of the helium exchange gas could be neglected at all except 
the lowest temperatures. 
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. The metals will now be dealt with individually, any special details in technique 
being described with reference to the particular metal concerned. 


; §2. PALLADIUM 

Keesom and Clark (1935) have shown that the atomic heat of nickel at low 
temperatures contains a linear term 0:001744 7 which is many times bigger 
than the contribution from free electrons to be expected from the Sommerfeld 
theory. Mott (1936) has suggested that this term is due to the unfilled quantum 
states or positive holes in the d-shell of nickel, and that a similar term should 
occur for the other transition metals, whether ferromagnetic or paramagnetic. 

In order to test this theory and to make more quantitative data available on 
the subject, the specific heat of palladium was determined, palladium being chosen 
on account of its similarity to nickel and its high paramagnetism which, according 
to Mott, should imply a large electronic specific heat. Though the electronic 
component in these transition metals is larger than usual, it is nevertheless small 
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Figure 1. Atomic heat of palladium. Figure 2. Atomic heat of sodium. 
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in absolute magnitude. But since it varies in direct proportion to the temperature, 
whereas the specific heat due to the lattice vibrations varies near the absolute 
zero with the cube of the temperature, the two become comparable in magnitude 
at sufficiently low temperatures as pointed out by Simon (1928). Experiments 
were therefore carried out in the liquid hydrogen/helium region. 

The sample of palladium was kindly put at our disposal by Professor N. F. 
Mott and was a solid cylinder of the metal, 1-4cm. in diameter and 1-4 cm. long, 
weighing 25-21gm. The palladium sample was strongly heated im vacuo to 
determine if it contained any absorbed hydrogen. The amount obtained was 
scarcely detectable and certainly too small to affect the results at all. Its atomic 
heat was measured in the small calorimeter in the range 2—22°K, 

The values obtained for the atomic heat are plotted in figure 1. A mean 
curve was drawn and was analysed by trying a series of values of the Debye 


Lae 
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characteristic temperature (6p) to calculate the lattice atomic heat (Cr) and 
finding which gave the most nearly linear remainder for the electronic specific 
heat (Cy). The best fit was found with @)=275°, and table 1 gives values of 


Table 1. Atomic heat of palladium (in cal/deg. gm. atom) 


Temp. (° K.) Gare CiG—275=) Cy/T 
2) 0:0066 0:0002 0-:0032 
4 0-0132 0-0014 0:00295 
6 0-023 0:0046 0-:0031 
8 0-036 0-0114 0:0031 

10 0-054 0:0224 0:00315 
12 0:076 0-039 0-0031 
14 0-106 0-062 0-00315 
16 0-142 0:0915 0:00315 
18 0-186 0-130 0:0031 
20 0-240 0-178 0:0031 
22 0-306 0-237 0:00315 


the observed total atomic heat (C,,,) from the experimental curve, with values 
of Cy, for 9@=275° and of Cy/T. The values actually measured are of the atomic 
heat at constant pressure, but since the correction C, —C, is very small at these 
low temperatures, the observed values have been taken as Cy. 

The observed atomic heat can therefore be represented as the sum of a T? 
term (lattice atomic heat) and a T term (electronic specific heat) in the range 
2°-22° K. as 


C=0-0000224 T3 + 0-0031 T cal/deg. gm. atom. 


In figure 1 a curve of this form is plotted for comparison with the experimental 
points. ; 

The value of 6)=275° from the present results may be compared with the 
value of 6) =270° given by Meissner and Voigt (1930) from electrical measure- 
ments and with that of 260° from the Lindemann melting point formula. 

In the expression for the atomic heat at 2:5°K. the linear term is about 
twenty times as great as the 7? term while the two become of equal magnitude 
at about i2°K. ‘The value for the electronic specific heat is, as predicted by 
Mott, much higher than that of nickel (Cy =0-001744 T) while work by Kok 
and Keesom (1936) has shown that another transition metal, platinum, also 
exhibits a large electronic specific heat which, however, is not as great 
(Cy =0-001607 7) as that for palladium. 

A short note on these experiments has been published in Nature (Pickard 
1936). 

A short time ago Clusius and Schachinger (1947) measured the specific 
heat of palladium between 14°K. and room temperature and they found very 
good agreement with our experiments in the overlapping region. In addition 
these authors have now all the data available to analyse the whole specific heat 
curve (including measurements of other authors up to 1000°x.), and in combining 
their values with those of the paramagnetic susceptibility they are now in a 


position to estimate the electronic specific heat over the whole temperature 
range between absolute zero and 1000°. 
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§3. SODIUM 
The sodium used was from a sample of pure metal, hydrogen free, which 
had been prepared and used for experiments with neutrons in this laboratory. 
18-90 gm. of sodium in the form of a cylinder which just fitted within the large 
calorimeter was used and measurements were made in the range 2—25°K. : 
The experimental points are shown graphically in figure 2 while some values 
for the atomic heat are given in table 2. 


Table 2. Atomic heat of sodium 


Lemp: ©K.) Cy Op (°) Temp. (° kK.) Cy 8p (°) 

Dy, 0-0055 88 10 0-140 149 
3 0-0080 116 il? 0-250 147 
4 0:0135 130 14 0:385 148 
5 0:027 129 16 0-54 150 
6 0-059 120 18 0:71 152 
7 0-124 109 20 0:90 154 
8 0-093 ISI 25 1:45 156 
9 0-100 150 


The present values are slightly higher than those of Simon and Zeidler (1926), 
who measured the specific heats between 17 and 60°K. Their values, which 
are shown in figure 2 together with the new ones, correspond to a6 of 159°. It 
can be seen from table 2 that the characteristic temperature becomes progressively 
lower with falling temperature. 

Doubt had been expressed by Simon (1926) whether the value of 6 of about 
160° really represented the characteristic temperature due to lattice vibrations. 
He was led to this suggestion by the fact that the characteristic temperatures 
derived from Lindemann’s melting point formula and from electrical conduc- 
tivities were considerably higher. In particular, the latter values could be 
interpreted on the basis of a 6 of 208°. Simon had suggested that, as in the case 
of grey tin and other substances discussed in his paper, this specific heat curve 
should be interpreted by a Debye term (@=202°) plus a “Schottky” term 
(@=95°) due to an internal transition of as yet unknown nature. Later, when 
Blackman (1934, 1935) had shown that the Debye function was less well founded 
than had previously been assumed, the opinion had been expressed that there 
was no need to assume the existence of such internal transitions, and the incom- 
patibility of the specific heat values with the Debye curve was attributed to the 
defects of the latter theory. 

While it is certainly possible that some relatively small anomalies in specific 
heat curves are due to this, we maintain the opinion that in the majority of cases 
discussed in the paper mentioned above, we are faced with internal transitions. 
A systematic investigation of the whole question has now been started at the 
Clarendon Laboratory, including measurements of thermal expansion, which 
will provide an independent value of 6. In the case of sodium, recent very 
careful experiments by MacDonald and Mendelssohn (1948) had shown that 
the electrical conductivity down to 2°K. can be represented very accurately 
with a @ of 202°. The difference between this value and the one derived from the 
specific heats without the assumption of an anomaly is much too big to be 
explained by a deficiency of the Debye theory. 

The most striking feature in the specific heat curve is the anomaly at about 
7° (see figure 3). An anomaly of a similar type had been observed by Cristescu 
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and Simon (1934) in beryllium at 11°. In both cases the total energy due to 
this anomaly is not of the order RT, but in the case of sodium RT/70 and in the 
case of beryllium R7/700. These are of the same order as the anomalies in the 
specific heat observed in superconductors, and investigation of the electrical 
properties (see also Kurti and Simon 1935) may shed more light on the nature 
of this anomaly. I 
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Figure 3. Atomic heat of mercury. 
© AE] Pickard and Simon ; *+ Simon; X Pollitzer ; °< Onnes and Holst. 


§4. MERCURY 


The object of the experiments was to investigate the possibility that a poly- 
morphic transition, as predicted by Bridgman (1935), might occur in the region 
below liquid oxygen temperatures. It would manifest itself in the form of a 
discontinuity in the total heat/temperature curve. It was at the same time 
desired to extend the previous measurements (Pollitzer 1911, 1913, Simon 
1923 a, 1923 b) down to the helium region. 

Bridgman’s measurements were concerned with the polymorphic transitions. 
of elements at very high pressures. An extrapolation of the data obtained by 
him for the temperature variation of the transition pressure between the poly- 
morphic forms of mercury indicated that this transition might occur at normal 
pressure at a low temperature. From an inspection of Simon’s measurements. 


(1923 a, 1923b) he claimed that a distinct break in the tangent to the atomic: 


heat/temperature curve was evident between 53°-68°xK. and that this might 
be evidence of the occurrence of the transition. 


The rather complicated crystal structure of mercury makes it improbable. 


that its specific heat could be represented by a single characteristic temperature ;. 
Simon (1923b) had found good agreement with the experimental data, using a 
two-term expression consisting of a Debye-function and an Einstein-function, 


the latter having a very low characteristic temperature, and, therefore, being 
responsible for the hump in the solid hydrogen region. A change in the curvature: 
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was therefore to be expected. In view of Bridgman’s prediction of a transition, 
however, it seemed desirable to check the specific heat again in the region between 
hydrogen and air temperatures. 

The method adopted in the search for any thermal discontinuity was to heat 
the calorimeter slowly and continuously by passing a constant current through 
the heater coil, and to measure and record the temperature at the usual 20-sec. 
intervals. A graph of temperature against time was then plotted and examined 
for evidence of any sudden change which would have been indicated by a dis- 
continuity in the curve. 

The sample of very pure mercury (75-80 gm.) was kindly lent to us by Dr. H. 
London and was contained in a vessel of thin copper, varnished to prevent amal- 
gamation of the mercury with the copper. his vessel was placed in the large 
calorimeter. ‘The measurements of the atomic heat were extended over the 
range from 3:5° to 95°K. The slow continuous heatings extended over the whole 
region from 25°-80° k. 

The experimental points for the atomic heat are given in figure 2 and a mean 
curve 1s drawn through them. Some values from this mean curve are given 
in table 3. 

Table 3. Atomic heat of mercury (Cy) 


ANssenvon (CF iy) (Cx Temp. (° kK.) Or 
329 0-12 25 3:18 

+ 0-20; 30 3-68 

6 0-68 40 4-4] 

8 0-87; 50 4-95 

10 1-06 60 5°30 

2 ies 70 Soy 

15 1-72 80 5-67 

20 2:46 90 5:80 


The present values for the atomic heat between 10° and 14°xK. and above 
25°x. are slightly lower than the previous ones of Simon (1923 a, 1923b) which 
are shown in figure 3 together with the results of Pollitzer (1911, 1913) and two 
early results of Kamerlingh Onnes and Holst (1914). The present measurements 
confirm that the atomic heat is anomalously high at low temperatures. ‘The 
results employing continuous heating were carefully examined but no evidence 
of any thermal discontinuity was observed within the limits of the experimental 
accuracy, about 0-02 cal/gm. atom in the region investigated. It must be men- 
tioned that although the temperature was only recorded every 20sec. it was 
actually measured almost continuously, and so the possibility of even a very 
small discontinuity passing unnoticed is very remote. 

It should be added that mercury becomes superconducting at 4:2° and that 
there is a small discontinuity in the specific heat at this point (Daunt and 
Mendelssohn 1937). 

§5. THERMODYNAMIC FUNCTIONS 

Our measurements at the lowest temperatures together with those of other 
observers at higher temperatures now permit calculation of the thermodynamic 
data of the three metals. Tables 4 to 6 present in the usual manner the specific 
heat (C,) entropy (S), enthalpy (J) and thermodynamic potential at constant 
pressure (F,) at intervals of 10°. Weare very much obliged to Mr. G. K. White 


for the trouble he has taken in calculating these data. 
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Table 4. Palladium (1 gm. atom = 106-7 gm.) 


TGR) a CPO 6 
OO De43ono? 
70! Beast S87 
its) Sog8}  OAY4: 
1908 35<7220"355 
200 02/22 40.6% 
210 5:86 6:93 
22 eo mee 20) 
DSN SA aed) 
240 6:02 7:72 
250 @-067 7-97 
260- 6107 “S-Z5 
270 6714" “Ste 
280 6:18 8-66 
290 6:22 8:88 
300° 6:27" 9-09 

gm. atom = 23-0 gm.) 


T Cm. Ce 
160 6-00 
170 6-07 
180 6:13 
190 6-18 
200 6-24 
210 6:30 
220 6:36 
230 6-41 
240 6:47 
250 6:53 
260 6:59 
270 6:65 
280 6-71 
290 6-77 
300 6-83 


S—So 
8:40 
8-76 
Qoilil 
9-44 
9-76 

10-07 

10-31 

10-60, 

10-87 

11-14 

11-40 

11-65 

11-89 

12-13 

12236 


Mercury (1 gm. atom = 200-6 gm.) 


S—Sy I-I, —(p—F»,) 
0:0 0-0 0:0 
0-038 0-110 0:27 
0-121 1-40 1:02 
0-30 5°71 Son7 
0:57 15:3 7°54 
0-93 Sil od) 15-15 
1-34 53:6 26°5 
iors, 81-9 42-1 
2°22 115 62:1 
2:66 153 86-1 
3:09 194 115 

3305) 238 148 
3:92 285 185 
4-31 334 226 
4-68 385 7A 
5:05 437 320 

Table 5. Sodium (1 

S=So I—Ih —(Fp—Fp,) 
0-0 0-0 0-0 
0-067 0-452 0-118 
0-38 5°23 2°31 
0:96 19-9 8°81 
1°67 45-0 21:9 
2-44 79-4 44-6 
3:19 120-6 70:7 
3-89 166-6 105-6 
4°54 215 148 
5-15 267 196 
5°71 320 251 
6:23 375 310 
6°72 431 375 
7:17 = 488 444 
7:60 546 518 
8-01 605 596 

Table 6. 

S-Sp J—-lh =F.) 
0-0 0-0 0-0 
0-70 4-62 2°26 
1:85 22-2 14-7 
3-11 53-6 39°5 
4:27 94-3 76:5 
5053) 141 125 
6:25 193 182 
7:08 247 249 
7°83 303 324 
8:51 360 405 
9-12 419 494 
9-69 478 588 
10-21 538 687 

OEY 599 792 

11-45 660 901 

11:58 722 1014 

11:98 784 1139 

* Solid. 


PACK) Cp 
170 6°31 
180 6:36 
190 6-41 
200 6:47 
210 6:53 
220 6:62 
230 6:74 
234:4* 6-78 
234-4+ 7:00 
240 6°95 
250 6-90 
260 6°85 
270 6-76 
280 6:70 
290 6:65 
300 6°65 


+ Liquid. 


S— So 
12:36 
12-72 
13-06 
1339 
I33e7/4 
14-01 
14°31 
14-43 
16-84 
17-01 
17-28 
17-56 
17-81 
18-06 
18-29 
18-52 


f-I, 
491 
Sy115) 
601 
658 
715 
ges 
833 
892 
952 

1012 

1073 

1133 

1196 

1258 

1320 


=I, 
665 
i225 
786 
848 
910 
972 
1036 
1100 
1164 
1229 
11294 
1360 
1437 
1505 
1573 


847 

Dit 

974 
1039 
1104 
1169 
1236 
1265 
1820 
1854 
1924 
1991 
2058 
PVA) 
2192 
2258 


— (Fp —Fp,) 
372 
428 
486 
548 
613 
681 
752 
825 
901 
980 
1061 
1144 
1230 
1317 
1407 


es Caen ain 
679 
764 — 
854 
947 
1043 
1142 
1233 
1338 
1445 
1556 
1668 
1784 
1892 
2012 
2134 


1254 
S79 
1507 
1640 
1775 
1914 
2055 
2116 
2116 
Drea 
2418 
2570 
- 2746 
2939 
3107 
S297 
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The Theory of the Propagation of Second Sound in Helium II 
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Communicated by N. F. Mott ; MS. received 20 January 1948 


ABSTRACT. In this paper the conduction of heat in liquid helium II will be discussed, 
and the velocity of second sound, the attenuation, and reflection coefficients under various 
conditions will be obtained. The shape of the velocity—temperature curve will be discussed 
on the basis of Landau’s theory of helium II, and 'Tisza’s theory will be criticized. 


§1. INTRODUCTION 

HE experiments of Peshkoy (1946) in Russia and Lane and co-workers (1947) 

in America have demonstrated beyond doubt the existence of a new type of 

sound wave in liquid helium II, one in which temperature fluctuations are 

the important phenomena and not pressure fluctuations as in the case of ordinary 
sound. The essential assumptions of the theories developed by ‘Tisza (1938, 1940 
and 1947) and Landau (1941, 1944) which account for second sound are as follows : 
(i) Helium II may, at least from a phenomenological standpoint, be divided into 
two parts. One of these parts, the so-called superfluid, is probably to be explained 
as some kind of condensed phase, perhaps the condensed phase to be expected 
below a certain temperature in a substance obeying Bose-Einstein statistics 
(London 1938, 1939). (ii) The entropy of one part, the so-called normal fluid, is 
greater than that of the superfluid. Kapitza (1941) has shown, by experiments in 
which the two fluids are separated by the use of a very narrow gap, that the entropy 
of the superfluid is very small: it may be zero. (it) To these assumptions may 
be added a third, which is important in the theory of the flow of heliumII. ‘The 
normal fluid obeys the hydrodynamical equations appropriate to any other normal 
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fluid; thus its equations contain a coefficient of viscosity ; but the superfluid obeys 
hydrodynamical equations appropriate to a fluid in which the curl of the velocity is 
everywhere zero, so that the motion of the superfluid is always irrotational. In 
particular, the tangential velocity of the normal fluid must vanish at a surface, 
whilst the tangential velocity of the superfluid need not vanish there. 

In this paper we shall consider several problems in the theory of the propagation 
of heat and of second sound in helium II using the above assumptions ,, 
the problems are chosen both because of their importance in experimental investi- 
gations and also as illustrations of the use of the equations and boundary conditions 
which will be established in this paper. 


§2. THE SHAPE AND SPEED OF ORDINARY 
“CONDUCTED? RUS E'S 

Throughout this section, the word ‘“‘conducted”’ will mean a flow of heat not 
involving fluid movements. 

When Tisza (1938, 1940) predicted the existence of “temperature waves”’, 
he pointed out that these waves could easily be distinguished from ordinary 
“conduction waves”’ by the fact that the latter would be greatly attenuated, and 
would also have a velocity which would be proportional to the square root of the * 
frequency, unlike the temperature waves, which should have a velocity 
independent of the frequency. 

For suppose that 7 is the temperature excess over the initial uniform value, ¢ the 
time, K the ordinary thermal conductivity (due to heat flow without fluid move- 
ments), and C the specific heat at constant volume. ‘Then elementary arguments. 
give the equation of heat conduction as dT/dt =(K/pC)V?T, where p is the density. 
Suppose for simplicity that we are dealing with plane propagation along the 
zg axis. Then this equation becomes d7/dt=Dd?T/dz* where D=K/pC. If 
the source consists of a plane at z=0 undergoing sinusoidal temperature 
fluctuations of frequency w/2z, then dT/dt=iwT; d?T/dz*=iwT/D; and 
T cce'exp { ~2(1 +2)/(w/2D)}. 

Hence the distance in which the amplitude of temperature falls to 1/e of the 
temperature of the source is just one wavelength, and the apparent velocity is 
4/(2Dw). 

If on the other hand the ‘‘conduction wave”’ is emitted from the plane z =0 at 
the time t=0 as a sharp pulse, the appropriate solution of the equation is 


T ocexp ( —2?/4¢D)/+/t, 


where D=K/pC. If the conduction effect is the only one of importance, or if all 
the temperature waves of the Tisza-Landau type are reflected back without loss, 
then the proportionality constant may be found by integration of T (the excess 
temperature) over all space. The result must clearly be independent of time 


except during the initial sharp emission of the pulse. Thus if QO is the quantity 
of heat in the pulse, 


00) 


D ANC 10 0) 
Q=C] Tas=— { exp (—22/4tD) dz = ACy/(4nD). 


where A is a proportionality constant. Thus finally 


7: ee 


Sei ae 
Ca/(4ntDy ©? | 27/4tD). 
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The maximum of the pulse, for a given value of g, is given by (07 /dt),=0; 
1.€. tax = 8"/2D and Ty, =0-2420/Cz. Thus the apparent velocity of the pulse 
depends on the distance z; this is simply because the conduction effect is dispersive. 
The shape of this ‘‘ conducted pulse” may be deduced from the above equations. 

In problems of this type it is only necessary to obtain solutions of the physical 
problem involved for the simple case of pure sine waves; all other cases may be 
deduced mathematically from this solution by means of Fourier analysis. 


$3. DHE VELOCITY OF SECOND SOUND 


A satisfactory deduction of the velocity of second sound was first given by 
Landau (1941). Recently a simplified approximate treatment has been proposed 
by Gogate and Pathak (1947). They obtain the same formula for the velocity, 
but some of their arguments appear unsatisfactory to the present author, especially 
as in the deduction of one equation second-order quantities (involving the product 
of the two first-order quantities dT and v, in their notation) are used. An alter- 
native simplified theory is therefore given below, and a more detailed discussion of 
the problem in §9. We consider that the treatment given by Meyer and Band 
(1947) is illogical, as they commence by writing down equations appropriate to a 
single fluid, and finally divide it into two effective parts. 

We may form the following physical picture of the wave motion associated 
with second sound: due to the relative motion of normal fluid which carries 
entropy, and superfluid which does not, we obtain a state of the liquid such that 
alternate layers contain an excess or a deficit of entropy. It follows that these 
layers are at an effectively higher and lower temperature respectively. In the 
following simplified deduction of the velocity of such waves, we shall neglect the 
correction caused by the small pressure differences which will be set up in the 
liquid due to the fact that the layers at the higher temperatures will contract 
slightly, and the layers at the lower temperature will expand a little, helium II 
having a negative coefficient of thermal expansion. We shall also neglect the 
difference between total and partial differentiation with respect to time, since this 
difference is quadratic in the velocities—that is, we neglect the difference between 
differentiation with respect to time whilst remaining stationary in space, and 
differentiation whilst following the movements of the particles. 

Let T be the temperature excess over the mean undisturbed value, p the density, 
S the entropy per unit mass of the liquid considered as a whole, p, the density of 
normal fluid, p, the density of superfluid, ¢ the time, V,, the velocity of the normal 
fluid, V, the velocity of the superfluid, and ¢ the thermodynamic potential of the 
fluid per unit mass. 

Then p=p,+pz, and for the current 7 we have j=p,V,+p,V,. If we assume 
for the present that entropy is carried only by the normal fluid (we return to 
this question in §9), we have the equation of conservation of entropy, 


A(pS)/dt + SV.(pV,,) =0. 
Since C= T’0S/OT, this equation may be rewritten (provided the fluctuations are 
small and density changes negligible) : 
OViot+(T S/C)\V.V.=0) ee sae: (1) 


Here 7” is the actual temperature of the fluid. 
Confining attention to waves in which the fluctuations in pressure and density 


are negligible (these pressure fluctuations being associated with ordinary first 
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sound), we see from the conservation equation 0p/dt+V.j=0 that j=9, as any . 
other constant would imply a constant net flow of the fluid in one direction. 
The final equation required is provided by an argument similar to that given by 
Landau (1941), deducing the equation giving the relation between force and acceler- 
ation for the superfluid alone. We imagine two thermally isolated vessels, 
connected by a capillary so narrow that only superfluid may pass along it. Then 
the net work required to transfer unit mass of superfluid from one vessel to the other 
will equal the difference 0¢ in the thermodynamic potential (Gibbs function) per 
unit mass between the two vessels. This work will give the change in potential 
energy. Thus if we consider again the bulk liquid, the gradient of the thermo- 
dynamical potential with respect to distance (taken with a negative sign) will give 
the effective force acting on the superfluid alone. Therefore, for the super- 
fluid, the equation relating acceleration to force is given by OV,/dt+Vd=0. 
Since Vé= —SVT+VVp= —SVT for displacements where Vp is negligible, 
éV,/at —SVT=0. Equating 7 to zero, we have p,V,= —p,V,, and hence 


OV,j0t+ Sp,V1j/p,=—0- a ee (2 
Combining this with equation (1), we obtain a wave equation 
OL LS ape 
ee VT, “SS on See eee 3 
of? GC Pn ~ ( ) 


with which there is associated a velocity u,=(T’S*p,/Cp,)?. 

We shall now show that in a travelling plane wave of second sound the heat 
flow at any point is proportional to the temperature at that point, and not to the 
temperature gradient as in a conducted wave, and also that in a stationary plane 
wave of second sound the heat flow at any point and at any time is $7 out of phase 
both in time and in space with the temperature at that point. These results 
follow from equation (2). The rate of heat flow per unit area is T’pSV,, and thus 
the rate of heat flow is proportional to the velocity of the normal fluid. Ina travel- 
ling wave we may take T proportional to cos (wt —yz), where y is the wave number 
w/u,. Thus 0V,/dt « dT /dz csin(wt —yz), and hence V,« cos (wt —yz) a T. 
In a stationary wave we may assume J proportional to coswtcosyz with 
suitable origins of time and space. Vis then proportional to sinwt sin yz. Thus 
the velocity of the normal fluid, to which the rate of heat flow is proportional, is 
37 Out of phase with the temperature both in space and in time. 


§4. THE ENERGY FLOW IN SOUND 


Energy considerations are valuable in the deduction of reflection coefficients 
and allied problems, and it is therefore important to determine the rate of energy 
flow per unit area in sound waves. 

In a travelling wave of first sound, the kinetic energy per unit volume is given 
by 3p, and if the wave is sinusoidal, the average density of kinetic energy is given 
by tp(V°)?, where V° is the peak velocity. It is easily verified that the density of 
potential energy is the same as that of kinetic energy in a plane wave of first sound 
($10); it follows that the total energy density is }p(V°)?. Therefore the rate of 
energy flow of first sound per unit area is given by 3p(V°)?u,, where u, is the 
velocity of the waves. Since in first sound the maximum excess pressure p° is 


given in terms of the maximum velocity V° by the equation p = pu, V°, we may also 
write for the rate of energy flow 4(p°)?/pu,. 
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In a travelling wave of second sound, the kinetic energy per unit volume is 
2(PxVn+psV3). Since pyV,= —p,V, this may be rewritten as 4p,pV/2 p, or 
tPnp(Vn)?/ps, where Vz is the peak value of the velocity of the normal fluid. It 
will be proved later (§ 10) that the density of potential energy has the same value, 
and hence we may deduce that the total energy density is given by 4,p(V/°)?/p,. 
Thus the rate of energy flow of second sound per unit area is dpnp(Vp)*U9/ psy 
where up is the velocity of this sound. Since in second sound the maximum excess 
temperature is given by T°=p,u,V;/p.S, we may also write for the rate of energy 
flow ppsS°(T°)?/2pyuy = puyC(T°)?/2T’, where 7’ is the mean temperature of the 
liquid. 

§5. THE ATTENUATION OF SECOND SOUND DUE TO OTHER 
HEAT PROCESSES 

It is interesting to attempt to calculate the attenuation of second sound caused 
by non-reversible conduction of heat. For the sake of simplicity only processes in 
which heat flow is directly proportional to the temperature gradient (with a factor 
of proportionality K per unit cross-section per unit time) will be considered. 
Thus included in K will be, fer example, the ordinary heat conductivity without 
fluid movement. 

In the equation of entropy conservation we must now consider also irreversible 
processes. If we consider unit volume, the rate of loss of entropy is —0(pS)/dt. 
This must be equated to the sum of the reversible loss pSV.V, (assuming 
that only the normal phase carries entropy), and the irreversible loss 
(1/T’)(eO/ot) = —(K/T’)V?T, where K is the irreversible conductivity defined 
above. Thus d(pS)/dt+pSV.V, —(K/T’)V?T =0 and, therefore, 

Clipe NaS i mle ee 

oe yg argh dies) 
Eliminating V,, between this and equation (2), which is unchanged to this approxi- 
mation, we find for the (complex) velocity 2, 
T'S" (1 Kiwpn ) 
Con pl'S*p,) ° 
Now the amplitude of the wave is proportional to exp (tw2/u.), of which the real 
part decays as exp ( —Kw*z/2u3pC). Thus the distance in which the amplitude 
falls to 1/e of its initial value is given by 2u3pC/Kw?. Since u,=0 at the A-point, 
we see that the attenuation increases rapidly in the region of the A-point. It also 
increases near the absolute zero since here C becomes very small, whilst w. should 
remain finite according to theory (§11). ‘Thus the attenuation should exhibit a 
minimum between 0°x. and the A-point. 

For the ordinary non-reversible thermal conductivity K it is reasonable to 
insert the value found for heliumI, where it is not masked by second-sound 
convection effects. Keesom and Keesom (1936) give K~6 x 107° cal/deg.cm.sec. 
at 3:3°x. For the specific heat we insert the smoothed values given by Keesom 
(1942). Then very approximately the distance in which the amplitude falls to 
1/e of its initial value is 2 x 10°°/f? cm. at 1:6°K. and 7x LO /f? cm. AtizeLS  K., 
falling rapidly to zero at 2:19°x., the A-point. It therefore appears likely that the 
attenuation due to irreversible heat effects is only of importance at very high 
frequencies except close to the A-point and probably, also, at extremely low 
temperatures. In the above, f is the frequency in cycles/sec. 


TB 
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§6. THE REFLECTION OF SECOND SOUND AT WALLS HAVING 
FINITE HEAT CONDUCTIVIRY 

The solution of this problem of reflection is important because one experi- 
mental method for measuring the velocity of second sound depends on the setting 
up of standing waves in a tube. Suppose that there is a source of temperature 
waves at one end of the tube, and a solid block of material (e.g. a thermometer) of 
thermal conductivity K at the other. As the heat conductivity, or rather heat 
convection, of helium II is so much greater than that of any solid, we shall assume 
that the solid block is effectively infinite in depth. In analogy with transmission- 
line problems we consider the problem of reflection at a junction of two “trans- 
mission tubes’’, the liquid being one and the solid the other. Let the temperature 
T be the counterpart of voltage, and the heat flux / the counterpart of current. 
The ratio 7/J across a junction will be called the impedance Z across that junction. 
Then the problem may be divided into three parts: first, the calculation of the 
characteristic impedance Z, of the liquid ; second, the calculation of the impedance 
Z, of the solid at the interface between solid and liquid; and finally, the deduction 
of the reflection and transmission coefficients from these impedances. 

For simple harmonic waves of frequency w/27, equations (1) and (2) of §3 may 
be written 0I/dx= —iwCpT and @T/dx= —iwI/pCus, where I is the heat flux, 
pol Vs 

Now for a loss-free transmission line of inductance L* per unit length and 
capacity C* per unit length, we have as equations linking voltage V and current I 

dLlox=—1aCVy eee (4) 
and GV jez = 1a L*ie = eee (5) 
Comparing these equations with those for our case, we may put C*=Cp and 
L* =1/pCuz. (From §5 we may also derive the value G=w?Kj,,/u3 for the 
effective resistance per unit length across the tube.) We may now take over the 
well known theory of loss-free transmission lines. For instance, we may write the 
characteristic impedance of the liquid as Z)=4/(L*/C*)=1/pCu,. The 
temperature 7(x) at a point x and the heat flux J(x) are given by 
T y= Tee Pee ere 
Zilia) = Le ap eae. Se fee (7) 
where T;and 7, are constants giving the temperature maxima of the forward and 
reflected waves respectively, and 6 is the wave number w/uy. These formulae 
satisfy equations (4) and (5), 

The standing wave ratio (i.e. the ratio of maximum to minimum) will be given 
by (7; | = | T, | iq T;| — | 7. | ). Using (6) and (7) to determine this in terms of Li 
and Z, (the impedance at the end x=/), we have for the standing wave ratio 
| Z,/Z,|. By a double application of (6) and (7) we may also show that the ratio of 
temperature to heat flux at a distance / from a junction, across which there is an 
impedance Z,, is given by 

iz a 6l+72Z, sin Bl 
Fo a2, sin pie Zecos ait: 
Next we calculate Z,, the impedance of the solid at the interface. It follows 
from §2 that the temperature in the solid at a depth z from the plane of contact is 
given by 
Toy = Toe! exp {—ay(1 +0}, 
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where y* = wp. Co/2K, psc, and C, being respectively the density and the specific 
heat of the solid. The rate of flow of heat into the block per unit area is then given 
by so = — KOT s01/02),-9 =4(1 +1) Tye’, where a = Ky = (wpyoiCyK)!. Thus the 
impedance at the interface is given by Z,=1/a(1 +2) 
| Z/Zo| =spC x2. 

Finally, we must find the reflection and transmission coefficients in terms of Ze 
and Z;. ‘These may be obtained with the aid of the two boundary conditions that 
the temperature must be continuous across the interface, and that the heat flow 
must be a continuous function of distance. 

Taking the origin of coordinates at the interface, the first boundary condition 
gives T;+ T,=Ty., and the second J,+J,=J,.. Using the relations T,;=Zpl,, 
T,= —Zol,, and T= Zl, we obtain a reflection coefficient 


1,/1,=(2Z, —Z)|(Z,+Zp), 


and a transmission coefficient 


. The standing wave ratio is 


Tyoilli=2Zo(Zy+ Zo) = 2Z of Z 


‘The (complex) transmission coefficient is thus given by 2«(1+2)/pCus. 


o/b en Pe CrLON TOR SECOND ESOUND Aiwa A 
LIQUID-VAPOUR INTERFACE 


In an experimental arrangement originally suggested by Onsager, and used in 
the experiments of Lane and co-workers (1947), stationary plane waves of second 
sound are set up in a tube of liquid helium II and, owing to imperfect reflection at 
the liquid—vapour interface at one end, first sound is transmitted in the vapour. 
The process may be pictured very simply as follows: the stationary waves of second 
sound in the liquid cause a periodic fluctuation in temperature at the surface of the 
liquid, since, as there is very little heat flow at the interface, this is a temperature 
antinode. ‘There is, therefore, a periodic fluctuation in the amount of vapour 
emitted from the surface; this produces first soundinthe vapour. In this section 
we demonstrate how the reflection coefficient for sucha process may be calculated. 

Suppose that a travelling wave of second sound of unit amplitude is incident 
on the liquid—vapour interface, and that a wave of amplitude R is reflected, the 
remainder being transmitted into the vapour as first sound. ‘The boundary 
conditions which will be assumed are: (1) the vapour pressure at the interface is the 
equilibrium pressure at the temperature of the surface of the liquid; (ii) the 
difference in the energies transmitted by the incident and reflected waves in the 
liquid is equal to the energy transmitted by the sound wave in the vapour. 

In the liquid the temperature at any point is given by Tyg = Tye'"(e** + RE**). 
Here w/2z7 is the frequency of the wave, k the wave number w/w, and K is the 
reflection coefficient. 

At the interface we have 7,=7)(1+R). Now from the first boundary 
condition we may deduce that the excess pressure in the first sound wave at the 
interface is given by ~py= 7y(1+R) 0p/0T, where dp/0T is the slope of the equili- 
brium vapour pressure with respect to temperature; this is given by the Clausius- 
Clapey1on equation 0p/0T = L/T(Vyap — Vig): 

The excess pressure at any point in the vapour ts given by p = py exp (wt —1hx). 
Using the results of §4 we deduce that the rate of energy flow per unit area in the 
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vapour is 4p,2/p,u,, and the difference in rates between the energy flow in the in- 
cident and reflected waves in the liquid is }p,u,CT3(1 —R?)/T’. From equation (5) 


ap\?(1+R)* _ piel, _ pe 
Ghewen 
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§8. THE PROPAGATION OF SECOND SOUND IN THIN TUBES 


For convenience in experimental technique, the investigation of second sound 
is carried out in tubes of diameter about 1-5cm. It is therefore important to 
enquire whether the velocity is appreciably affected by the narrowness of the tube. 
The physical reason of the expected effect is that the normal fluid has a viscosity 
and, therefore, sticks to the walls of the tube, whilst the superfluid has no viscosity 
and, therefore, does not stick. 

Consider a plane wave propagated down a circular tube of radius 7 with viscosity 
of the normal part of the fluid wu. Then the force acting on the superfluid is that 
due to the gradient of the thermodynamic potential per unit mass, assuming the 
entropy to be zero, and the force acting on the normal fluid is that due to the 
sticking effect. ‘The mass-acceleration of the fluid as a whole will be equal to the 
sum of the sticking force and the gradient of pressure. 

By an argument similar to that given by Rayleigh (1926), we obtain the equation 


Op (,,2 fw), 20 [yw _vo 
sa(1+> /oe)+ oe [MS =v, 


This equation should hold provided the tube is sufficiently wide for the thickness 
adhering to the walls to be small compared with the diameter of the tube. Pro- 
vided the correction terms we are calculating are small, we may rewrite this 
equation as 07p/dt? =(1/f)V2p, where 


f=142/r)V/(u/2ep) + (2/ieor) /(uco]p). 


This equation, with 0?,S/dt=S*p.vy?T/p, ft, gives the velocity of the two 
sounds. By a slight extension of the work of Landau (1941), we obtain 


1./dp\: TS) > <EeS2hpop 
oi) Se) BS 
F\e@p/s Cen Conf Op 


According to this equation, the velocity of first sound is given by 


n= /G2)=2 a ae 
; fop]/~ dp TA} 2wp]’ 


and the velocity of second sound by u,=(T’S?p,/Cp,)?. Thus the correction to 
the velocity of second sound, on account of the narrowness of the tube, is much 
smaller than that to the velocity of first sound. f 


+ This equation (cf. (3) § 3) is true even in the case where pressure differences Aas 
are not negl 
(Landau 1941). gligible 
{ Note added in proof. A more detailed consideration now shows that there is a very small: 
correction proportional to the product of the coefficient of thermal expansion and the square root 
of the viscosity. 
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§9. THE USE OF THE LAGRANGIAN AND THE TOTAL ENERGY IN 
Lite DEDUGCTION-OF THE VELOCIIYeOLresSECOND SOUND IN 
A FLUID CONSISTING OF PARTS WITH TWO DIFFERENT 


ENTROPIES 

‘There is some question as to whether the entropy of the superfluid in helium II 
is really exactly zero; we shall therefore discuss and apply a method which is an 
extension of that of Tisza (1947). 

Let the entropy of the normal fluid alone be S,,, and the entropy of the super- 
fluid is S,. ‘Then the measured entropy S of the whole fluid is given by 
pS=pyS,+p,S;. Now we introduce new coordinates € and 7 such that 
pf =PrV_+psV, and 7 = (PnVnSn t+ psVgS5)/(PnSn +pgsSs)—§. Since is defined 
as the excess over &, the velocity of the fluid as a whole, of the velocities of the 
normal and the superfluid, weighted according to the amount of entropy which they 
carry, the equation of the conservation of entropy is expressible in these new 
coordinates in the simple form 0S/d¢+ SV.»=0. Also in these coordinates, the 
equation of the conservation of mass becomes (1/p)(dp/dt) + V.€=0. After some 
reduction, using the fact that p=p,+p,, we obtain for 7, S2p2n =p,p.S*(V, —V.), 
where S*=S,—S,. Thus the transfer of entropy depends on the relative 
velocities of the parts with different entropies, as is clear physically. From these 
equations we may obtain V, and V, in terms of € and 7 as follows: 


Vr=E+(mpS/prS*); Vs=& —(npS/psS*). 
Hence the kinetic energy per unit volume is given by 
EpVit dpsVe = 3pS* + 3(p2S?n?/pupyS*”). 


The potential energy is due to fluctuations both in volume and in entropy, 
Now if 5U represents a small change in internal energy, we may expand in terms of 
dV, the small change in volume, and 4S, the small change in entropy. Retaining 
terms up to the second order, we have 


2 2 27] 
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Thus the minimum work to cause the fluctuations is given by (Landau and 


Lifshitz 1938) 
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Here » is the coefficient of thermal expansion, which is small, and may be neglected ; 
the inclusion of this term gives the coupling between first and second sound. 


Thus the potential energy 1s, per unit mass, 
2 92 2 
a) (3 « BES. 
Op/s\ p KOT 
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It is convenient to rewrite the equations of continuity and of entropy conser- 
vation in terms of displacements rather than velocities. If we suppose that 
is the displacement corresponding to the velocity €, and 7 is the displacement 
corresponding to the velocity 7, we may write the equations of continuity and of 
entropy conservation as 

dp/p= —V.&y and 6S/S= —V.n,. 


The potential energy per unit volume becomes in these coordinates 
p (oP 2, pls® 2 
5 (22) (va) Ets VE (V.71) ) 
and the kinetic energy per unit volume becomes 


: ly SCY NO, 
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where dots denote differentiation with respect to time. Thus the total energy 
density is 


oy Pan aida \caee 0 Tas 
sob (se) at +3 (SE) (wees oo oom 


and the Lagrangian density is 


. pe [S\%. ap pl’ S? 
wits (S) at-a(Z) wea my. 


PnfPs 
The appropriate Lagrangian equations of motion are 
dt\0&,) 0&, o(V.é1) dt\0n;] Om o(V.n,) © 


Here dL/c(V.&,) anddL/0(V.7,) may be interpreted as pressures, and the gradients 
of these pressures give rise to forces. We obtain then the wave-equations 
£,=(€p|ap)sV%, and (p%/pup,)(S/S*)? iy = (7 S2/CWVyy. 

The velocity of the first sound associated with the first of these equations is, 
therefore, wu, =(0p/dp)*, and the velocity of second sound as given by the second is 
Uy =(TpapyS**/Cp*)?. 

We must again emphasize at this point that S* =S, —S, is the difference in 
entropies between the normal and the superfluid when account is taken of the fact 
that the fraction of the normal fluid is p,/p and the fraction of superfluid is p,/p. 
If we could isolate the superfluid, and it had the same density p as the bulk fluid, we 
could measure its entropy and call it, say, Sj. Similarly, if we had only normal 
fluid of density p, its entropy might be called S,. We should then have p,S, =pSq 
and pySa=pS». Then mee 

S*=S, —S; =p(Sn/Pn Sapa: 
In this notation 


ua =(T’/C){p(Si)?/Pn + Pal Ss)?/Ps —2SnS3}. 


This formula reduces to the well known one u5 = T’'S’p,/Cp,, when any entropy 
associated with the superfluid is ignored. ; 3 


S10. THE AVERAGE EQUALITY JOEVE ODEN AAI AND! Kol N Baie 
ENERGY IN A PLANE WAVE OF SECOND SOUND 
We may now prove the proposition of § 4 that the time averages of the potential 
and kinetic energies are equal. ‘This proposition may be stated in the alternative 
form of the vanishing of the time-average of the Lagrangian. 
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The Lagrangian density is, as shown in the last section, 


+9 2) Se Ves dase 
teeta V,e0E | + | pe ee Pg nye 
EE 2P . at £1) ] sls E Ae (=) Nee C (V.71) iL 


the terms in the first bracket being contributed by first sound and those in the 


second by second sound. If we suppose that €, and y, are displacements in plane 
sine waves, we shall have 


€,=€7sinw(t—x/u,) and 7, =7)sinw(t —x/Up), 


and thus = 4(@64)* and ne = noni)”. 
Also (8&,/8x)® a 5(w u 1/4)” = Hie 
and (An, /Ax)? = $(w7}/ty)? = 73/U2- 


Using the values deduced for u, and uy in the last section, we see at once that 
each of the two brackets forming the Lagrangian density vanishes when averaged 
over time. Since by Fourier analysis any plane wave may be split into a series 
of pure sine waves, this result holds for all plane waves. The result is thus a 
special case of the general proposition of the extremum of the time-integral of the 
Lagrangian function. 


Sn Ee SHAPE OR vob VELOCITY IE MPBRATURE CURVE FOR 
SECOND SOUND ON THE BASIS OF LANDAU’S THEORY 

Landau (1941) used the formula uw3=7'S%o,/Cp, for the velocity of second 
sound with calculated as opposed to experimental (or empirical) values of the 
thermodynamical quantities involved. ‘These calculated expressions contain at 
least two adjustable constants (arising in the roton contributions in Landau’s 
terminology: in the case of a sharp energy gap as considered by Landau in his 
1941 paper, just two constants, A and «, are required, A giving the magnitude of 
the energy gap and » the effective mass of the roton). ‘The statement made origin- . 
ally by Lifshitz (1944) that Landau’s expressions give for the velocity of second 
sound a monotonic decreasing function of temperature is only true for certain 
ranges of values of A and yw. It will now be shown that if the réle of phonons 1s 
iznored (using Landau’s terminology: phonons are Debye waves), the velocity at 
about 1-5°K. should rise with temperature, and this value should be almost 
independent of the particular value of Achosen. It is therefore possible to choose 
wand A in such a way that the roton contributions overweigh those of the phonon 
contributions and yet still give the correct magnitude for the velocity of the sound. 
In this way it is possible to obtain a dependence of velocity of second sound on 
temperature far from monotonic, having a new maximum between the A-point and 
the maximum at absolute zero. [tis impossible as yet to deduce A except crudely, 
as inertia (Andronikashvilli 1946) and specific heat experiments are not sufficiently 
accurate; as has been pointed out already, the velocity of second sound is largely 
independent of A in the range of temperatures so far employed for these experi- 
ments. 

To deduce this from Landau’s theory (1941), we start with the formula 
uz = TS2p,/Cp,.. Now at about 1:5°K. ps=p, the factor ps/p, merely giving the 
sharp decrease of velocity near the \-point. From Landau’s paper we may deduce 
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Thus, when only roton contributions are of importance, 
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Since KT <A at T~1-5°K., we see that under these conditions u, must rise with 
temperature. 

Landau’s results that the maximum value of the velocity is at absolute zero, and 
is equal to w,/4/3, will hold provided (i) any roton contributions have vanished due 
to the energy gap in the roton spectrum, and (ii) one part of the liquid has entropy 
zero and the other has an entropy due to Debye waves (phonons). 

The temperature, according to Landau’s theory, at which the rapid rise in the 
velocity of second sound takes place (at these temperatures close to the absolute 
zero) is governed chiefly by the value of A, which, as we have seen, is very difficult 
to derive experimentally from other sources. A value as low as about half a degree 
would be quite a possibility without contradicting present experimental evidence. 
Even this assumes the existence of a well-defined energy gap A. In this latter 
connection, a recent note by Landau (1947) is interesting. : 
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We only consider those aspects of Tisza’s paper (1947) which directly concern 
the problem of the propagation of second sound. 

Firstly, there is the problem as to whether it is possible for Debye waves to pass 
through the superfluid. Unlike Landau (1941), who includes phonons (i.e. 
Debye effects) as part of the normal fluid, 'Tisza considers that Debye waves are 
associated with the liquid as a whole, and that their presence will cause an attenu- 
ation of second sound. However, in his quantitative theory Tisza neglects their 
effect. If'Tisza’s viewpoint is correct, his theory must break down at temperatures 
below about 1°x., and his derivation of the zero value of the velocity of second 
sound at absolute zero (p. 852) and equations (11) and (12)—dp,/p,=dS),/S);3 
Pn/p = S,,/So—appear therefore to be invalid. Further, in the deduction of these 
equations Tisza considers the entropy as trapped in a cell possessing a membrane 
permeable only to superfluid. ‘The argument of p. 845 is only true provided no 
heat is supplied from outside, since any heat supplied would create more normal 
fluid and, therefore, more entropy in the cell. Under the adiabatic conditions the 
equations are thus correct, but equation (12) will merely tell us the temperature in 
the cell, and will not, as Tisza claims, apply to equilibrium values at any fixed 
temperature we like to name. Hence we do not agree that Tisza’s equation (13) 
{Pn/p =(T/Ty)", where Ty is the temperature of the A-point} follows from (8) 
{S,,/So=(T/To)"}, assuming that S,, and p, refer to equilibrium values at any given 
temperature. Also it is not clear physically why there should exist any thermo- 
dynamic relationship between the equilibrium values of S, and p, for any given 
temperature; one would expect the actual relationship between them to be specific 
to the detailed mechanism involved. 

The values which are required in the formula for the velocity of second sound 
are ordinary equilibrium values of S, C and p, for a given temperature 7’. These 
mean equilibrium values,+ about which the adiabatic fluctuations comprising 
second sound occur, are not directly connected with the values relating to Tisza’s 

} Actually not all the variables in this formula for the velocity of second sound are independent. 


Since C=T’(0S/0 T) and ps= p— pn, we may transform the formula to become u3=(1— p/pPy)/O(A/.S)/OT. 
Thus the mean equilibrium values required for any given temperature are those of py and o(1/.S)/aT. 
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“‘isopycnic”’ compression experiment. Thus it is incorrect to apply equation (12) 
(as Tisza does) to obtain the value of the velocity of second sound at 7'=0, and it is 
also irrelevant that Landau’s phonon equations violate this equation, as Tisza 
points out on p. 852 of his paper. It must be added that the relationships derived 
by Landau (1941) for his roton effects also conflict with Tisza’s equation (12), as 
may be seen from the values given in § 11 of this paper. (The entropy and normal 
density of an ideal Bose-Einstein gas do happen to obey Tisza’s relationship, but 
it is well known that the values calculated on this very crude model are quite 
inaccurate, the entropy being given as proportional to T°, as opposed to the experi- 
mental relation Sa T*°.) 

‘The present situation with respect to the temperature variation of the velocity 
of second sound might be summed up by saying that Landau and Tisza now agree 
on the formula v5 = T\S%p,/Cp,, although Tisza believes that the value of S inserted 
in this formula should refer not to the total entropy, but to this value less the Debye 
contribution. Tisza also considers that the formula will become inadequate 
below about 1°K., where these Debye effects become important. This author 
relies on empirical formulae for the entropy, with the constants determined from 
Kapitza’s experimental values (Kapitza 1941): he derives his values of p, from this 
entropy formula by the application of his relationship p,/p =.S,/S>, which, as we 
have shown above, has not the thermodynamic justification he has claimed. 
Landau goes further and attempts a theory involving two parameters, A and p, 
capable of giving the values of S, C and p, required, A and p being determined from 
Keesom’s values of the specific heat (Keesom 1942). 
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aroinatic systems, such as benzene or the polyenes, the method of molecular 

orbitals (m.o.) and the method of valence-bond structures (v.b.) agree 
remarkably closely. This agreement includes bond orders, bond lengths, 
free valence and resonance energy. However, during a study of the fulvene 
molecule C;H,.CH, we have found a large and serious discrepancy between 
the predictions of the two methods of approximation. ‘These discrepancies 
are large enough to be used to distinguish between the reliability of the methods 
by an appeal to experiment. Unfortunately the structural data are not yet 
available. But in view of the desirability of obtaining as many comparisons 
and tests as possible between the theories, it seems good to place on record the 
calculations that we have made, in the hope that they may stimulate experimental 
work with which to compare them. ‘The calculations themselves are of quite 
standard pattern and we shall not therefore describe them in any detail. : 

In the v.b. calculations, resonance is included between five canonical struc- 
tures. In the first place these are taken to be A —E. 


5d6 64 a6 


Figure 1. 


I is a well known fact that in dealing with most common conjugated and 


The normalized wave function for the ground state is 

Y= 0-6324 ob, —0-2109 (bg +4) +0-4217 (tp tug). «20... (1) 
‘The resonance energy, 0-5, has already been given by Wheland (1941). Bond 
orders may easily be inferred by the method of weighted contributions, as suggested 
by Pauling, Brockway and Beach (1935), or by the method of Penney (1937), 
using the Dirac coupling formula. The values are shown in figure 2, with 
the labels P.B.B. (i) and P.D. respectively. 

The m.o. calculation of bond orders are made by the method of Coulson 
(1939). ‘These values are shown in the figure with the label C. The resonance 
energy, 1-478, has already been given by Wheland (1941). 

Several comments suggest themselves at once from a study of figure 2. Firstly, 
the differences in bond order are much greater than is usually found, and even 
the sequence of orders is not preserved in any two of the diagrams. So far as 
the P.B.B. calculations are concerned this may be attributed to the impossibility 
of making any really satisfying choice of canonical structures. An objection 
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to the choice A-E is that none of them shows the bottom bond as double, so that 
in the P.B.B. method it cannot have any double bond character. The only 
way to avoid this is to make a different choice of fundamental structures. Any 
independent set of five will do, even if, in the language of Pauling, they are not 
canonical. This suggests that we introduce F and G. The appropriate modi- 
fication of the wave function (equation (1)) is easily made for any independent 
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Figure 5. Free valence. 
(a) Molecular-orbitals F.V.=4-68—/. 
(b) Valence-bond, relative to an ethylenic carbon, so that F.V.=4— 


set, using the ‘‘un-crossing rule”. In this manner we have used three reasonable 


sets, viz. (i) A-E, (ii) A, B, C, F and G, (iii) A, D, E,F, and G. Of these only 
(i) is a canonical set, but (iii) has the advantage of not using structures with more 
than one formal bond. ‘The results of our calculations in case (i) are already 
given in figure 2 : the results for cases (ii) and (111) are shown in figure 3. Even 
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the most cursory glance at these diagrams shows conclusively that there is not 
even the remotest agreement between them. It is quite evident that in molecules 
of this kind, in which ‘chemical intuition” does not suggest one single sensible 
choice of structures, the whole idea of weighted contributions is inadequate. 
Indeed, it suggests that even in other molecules such as anthracene or chrysene, 
although the choice of structures is reasonably straightforward, we must treat 
all calculation of bond orders by the P.B.B. method with some reserve. ‘The 
same is not true for bond orders calculated according to the P.D. method, for 
these are invariant for any choice of five independent structures, and we obtain 
the same values (figure 2) whether we use the set (i), (ii), (iii) or any other inde- 
pendent set of five. 

There still remains the large discrepancy between the P.D. valence bond 
and the m.o. bond orders. It is here that the distinction between the two methods 
is most clearly seen. For the m.o. method predicts quite large charge shifts. 
In its conventional form, it discusses the motion of the z electrons in a potential 
field supposed to be equivalent around each nucleus. But if we calculate the 
resulting charge distribution we find (figure 4) a large drift of electrons into the 
five-membered ring. For that reason the potential field around each nucleus 
cannot be equivalent, so that, as has already been pointed out in general terms 
(Coulson and Rushbrooke 1940), the m.o. method is not genuinely self-consistent. 
It is true that the shifts are probably rather less than figure 4 would indicate, 
but they should be great enough to provide a reasonably large dipole moment. 
On the other hand there is no dipole moment whatever in the v.b. calculations. 
‘Thus, if ionic terms are as important as the m.o. method suggests, then the v.b. 
method is inadequate in its simple form. But if they are relatively unimportant, 
then the v.b. method is the more reliable. An experimental determination 
of the dipole moment would settle this issue. 

A further discrepancy between the two methods is found in the calculated 
energies. For a large number of molecules (Wheland 1941), the ratio of the 
resonance integral (6) to the exchange integral («) required to give identical 
resonance energies has the value B/x=0-54. But in the case of fulvene it is 
0-34. A really satisfactory heat of combustion or hydrogenation could settle 
which method was more nearly correct. 

Our final comment concerns the “free valence” (F.V.) at any particular 
centre (Coulson 1948 a). According to the method of m.o., this is given by 


F.V.=4-68 —Y, 


where ./ is the total bond order (o +7) of all the bonds leaving the given atom. 
Values obtained in this way are shown in figure 5(a). If we may suppose that 
a high free valence is indicative of reactivity, then by comparison with benzene, 
for which the free valence is 0-35, it is seen that the non-ring carbon atom is 
extremely reactive. Some measure of polymerization, with or without dismuta- 
tion, would therefore be anticipated. This end carbon atom may be compared 
with the corresponding atom of benzyl, for which the free valence (Coulson 
1948 b) is only slighter greater, at 1-04 instead of 0-92. 

In cases where the choice of canonical structures in the resonance method 
is fairly evident, it is possible (Pullman 1947) to use the method of weighted 
contributions from each member of the set to estimate free valence, in a manner 
similar to that previously described for bond orders. But in the case of fulvene, 
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the reasons which we have given for disregarding the bond urders thus calculated, 
apply equally to the free valence, and we shall not therefore use this method. 
A similar objection does not obtain (Moffitt 1948) with the P.D. type of calculation, 
for which a description of the free valence may be given, parallel to that implied 
in equation (2) above for the molecular orbital theory. But since the greatest 
possible value of total bond order at a carbon atom (corresponding to 4-68 in 
m.o. theory) has not yet been settled for v.b. theory, our best way is to measure 
free valence relative to an ethylenic carbon. The resulting values are shown 
in figure 5(b). The high reactivity of the end carbon compared with the other 
carbons is now not at all noticeable and the difference between figures 5 (a) and 
5(b) is marked. A more adequate experimental study of the chemistry of 
fulvene should help to resolve the question which of 5 (a) and 5 (b) better represents 
the molecule. 
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ABSTRACT. The a-binding energy as a function of the number of neutrons and protons 
is shown to exhibit sharp fluctuations which occur at more or less regular intervals. ‘There 
is evidence that these ‘“‘ kinks’? are caused by discontinuous changes in the density and 


radius of the nuclei. 
It is shown that a mathematical connection exists between “‘ kinks’ in the a-binding 


energy curve and the irregularities of the f-stability line (Gamow valley), so that both must 
occur in the same region of the nuclei. This applies also to the irregularities in the mass- 


defect curve. 
A study is made of the energies of B-decay (AM) of isobaric nuclei in the region of the 


radioactive elements. ‘This shows that nuclei of odd mass-number A have a slightly 
higher binding energy for even proton number Z (a lower mass by 0:2 Mey.) than for 
odd Z. The effect, though well known for nuclei of even mass-numbers, has not before been 
observed for odd ones, where it is much smaller. 


§1. REGULARITIES IN THE a-STABILITY OF NUCLEI 
HE purpose of this section is to draw attention to certain regular variations 
in the «binding energy of nuclei which are superimposed on the general 
trend of ‘decreasing «-stability with increasing mass and charge. Being 
related to the gradient of the mass-defect curve, the «-stability is better suited 
than the mass defect itself for exhibiting irregularities which point to structural 


features in nuclei. 
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We shall employ the usual terminology, denoting the atomic weight by A, 
the atomic number by Z, the number of neutrons by N(=4 —Z), and the isotopic 
number by J (=N —Z=A —2Z). 


(i) Data for heavy nuclei 
One of the well known features of all the radioactive series is the extremely 
marked change in the «-energies between the isotopes of atomic number Z = 82, 


Figure 1. «-energies of large nuclei (Mev.). 
Heavy circle: Stable nuclei. Light circle: B-active nuclei. 
Light circle surrounding figures: «-active nuclei with or without B-decay. 


which are, within the limits of experimental accuracy, completely stable against 
a-decay, and those of atomic number Z=84, which have the most short-lived 
bodies and emit the most powerful «-particles. The average rate of change 
of the a-energy FE, in this region is (AE,/AZ)y— const, = 2:1 Mev. per proton. 
After element 84, the a-energy increases only slightly, [AE,/AZ]y, varying 
from 0.25 to 0 to 0.5 Mev. per proton. 

This striking discontinuity at the proton number Z=84 has given rise to. 
much valuable discussion (notably by Feather 1945), during which, however, 
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it has apparently been overlooked that a discontinuity of similar magnitude and 
importance exists also at the neutron number N=128. If we plot the energy 
of the emerging «-particles against the neutron number N and draw lines of 
constant Z, we get a sharp rise in EZ, between N= 126 and N=128 with a gradient 
of [AEF,/AN]z~ const, = 1:9 Mev. per neutron, followed after N= 128 by a constant 
decrease [AE,/AN],= —0-5 mMev., which lessens gradually after N=138 until 
it reaches.a value of —0-1 Mev. at N= 144. 

No systematic diagram of the «-energies is therefore quite complete which 
does not show this dual system of discontinuities at Z=84 and at N=128. It 
would appear that the best way of representing the facts is shown in figure 1. 
This represents a map of nuclei, using Z and N as abscissa and ordinate res- 
pectively, and the a-energy is indicated for every nuclear species. We can then 
draw contour-lines of constant «-energy, treating that quantity as if it were a 
continuous variable of Z and N. he diagram is ideally suited for the inter- 
polation and extrapolation of unknown «- and f-energies. 

To give an example, it is possible to predict with a fair degree of certainty 
that the only f-stable «-emitter of element 97 will have mass 247, an E,, of 6:0 
to 6-1 Mev. and a half-life of the order of one year. As the highest mass number 
so far reported in a trans-uranic element is 242, it is thus not likely that this new 
species will be discovered in the near future. 

Though in figure 1 the discontinuities at Z=84 and N=128 are very clearly 
noticeable, the general impression is that these discontinuities are not as important 
as they appear if Z or N alone is chosen as variable. It rather appears that 
the lines of constant #, form an irregular wave, with its axis somewhat steeper 
than the band of the known «-emitting nuclei which follow approximately the 
line of maximum f-stability (dotted line). The question immediately arises 
whether this E, ‘“‘wave’’ is confined to the radioactive elements or whether it 
can be traced also in the binding energies of the «-stable nuclei. 


(ii) Data for light nuclei 


For the lightest elements up to a mass-number A of about 40, the nuclear 
masses VV are well enough known to permit a fairly accurate calculation of the 
binding energy E,, for the last «-particle (which is the a-energy with opposite 
sign), according to 

EF, =$-3M+35M -7M. 


With reduced accuracy these calculations can be extended to A=60. ‘The 
mass data have been taken from an isotope chart compiled by Perlman, Seaborg 
and Segre in May 1947, and the «-energies are given in milli-mass-units 
(m.m.u.). 

Plotting again lines of constant £,, in a (Z—N)-diagram (figure 2), we find 
that, here too, the E,-lines perform irregular but more or less parallel wave motions, 
with the wave-axis approximately parallel to the Weizsacker line of stable 
elements. The bulges near A = 19, 39 and 59 represent minima of the a-binding 
energy differing by as much as 8m.m.u. from the normal levelomel.. The 
gradients [AE,/AZ]y and [AE,/AN], near the minima of EF, are actually of the 
same order as those observed with the radioactive elements near Z7=83, 1.e. 


about 2 Mev. per added nucleon. 
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Figure 2. «-binding energies of small nuclei (in m.m.u.). 


Double circle: Stable isotopes. 
Single circle: Isotopes unstable against B-, B+ or K-decay. 


(111) Discussion 
There is a considerable temptation to interpret these recurring fluctuations 
‘of the «-binding energy, which must be expected also near 4=90 and 140 (viz. 
the «-activity of 8Sm), in terms of a “shell”? model, though this need not 
necessarily be of the types considered by Elsasser (1933, 1934) and by Margenau 
(1934). With a shell structure it appears plausible that the addition of both 
neutrons and protons to a complete shell has at first an unstabilizing effect. 
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The «-stability would decrease until the newly added nucleons are numerous 
enough to interact with each other with a sufficiently large binding energy. 
From this point onwards the «-stability should increase rapidly to its normal 
value and remain there—with minor fluctuations——until the shell is filled, after 
which the process repeats itself. 

For small nuclei the electrostatic forces acting on a proton are small, and there 
is thus not much difference between the addition of a neutron or of a proton. 
The maxima and minima of E, should, therefore, be reached at a given number 
of nucleons independent of whether N or Z or J is kept constant. This is borne 
out well in figure 2. A similar picture can be applied to the elements of high 
atomic weight, but here the electrostatic effect of the high nuclear charge, resulting 
in a tilt of the curve, comes significantly into play; consequently minima of the 
binding energy (i.e. maxima of the «-energy) are obtained only for lines of 
constant Z or J, and no longer for constant N. 

While this ‘“‘shell’” model must be considered only as tentative, it should 
be pointed out that it does explain, at least qualitatively, the very puzzling 
feature that, at the proton numbers Z=8, 9, 17, 18, and 26, the addition of suc- 
cessive neutrons reduces the «-binding energy, and that, at the neutron numbers 
N=10, 11 and others, the addition of successive protons—in spite of their 
positive charge—increases the binding energy of the last «-particle. 


(iv) The case of argon 39 

The «-binding of ?*A, which in figure 2 is given as approximately 5 m.m.u., 
has been derived in the following way. It follows from the mass difference 
of 38S (34-9788) and 33Ca (43-9723) that (334, + 73E,) = 14-3 m.m.u. 

From its position in figure 2 one would expect 33E{, to be about 9m.m.u., 
which leaves for 33E;, only about 5:3m.m.u. But even if this value of EF) for 
839A is not accepted, it is apparent from the f-energy of {3S of only 0-11 m.m.u., 
and from the fact that 399A must be unstable, that }2h,<{8F,. This fact is of 
considerable interest, since all theories of the nuclear masses based on a nuclear 
radius r which grows smoothly with increasing mass (ro A?*), require that 
isobars have an «-binding energy which increases with the isotopic number J 
(Weizsacker (1935), Wigner (1937, 1939), or Feenberg (1947)). The case of 
389A obviously contradicts this assumption of a smoothly growing radius and 
requires that a sharp expansion of the nuclear radius takes place in the region 
of 2°A. 

To a minor extent this phenomenon occurs also at the other minimum, at 
A=17 and 19. Here the differences between the £;s of isobars for J= +1 
and J= —1 are only about half of the average AE,, between A=21 and 37, 
while according to all the theories they should be larger than the AE;s of subse- 
quent nuclei (AE; approx. oc A~*). uly 

In addition to this it is known from the application of the Gamow equation 
to the data of the decay constants and «-energies of the radioactive. elements, 
that a sharp expansion of the nuclear radius takes place in the region of the 
elements between Z=81 and 83 (Feather 1945), which is exactly where the 
biggest rise in the «-energies takes place. We therefore conclude that the 
minima of the «binding energy shown in figure 2 are also caused by a sudden 
expansion of the nucleus with increasing mass. While the nuclear radius. 
seems to grow, on the average, with A, the possibility must be considered that 
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it grows only little, if at all, in 
the interval between two [ 
“bulges” of the a-binding 
energy. ‘This means that the 
density of the nucleus must be 
expected to undergo regular 
‘fluctuations, minima of the 
density coinciding with 


minima of «-stability. 100 


T (sec) for even Z elements 
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There is some _ other 
evidence which supports the m1 
conception that the fluctua- va ea eee 
tions in the «-binding energy 00 # oe ges 
are not a purely energetic (? PA: 
feature but are due to changes Ue aS rine “i 


of the internal structure. ; Z E,* (in Aan 


If the Sargent curve of Figure 3. Deviations from Sargent curve near d=19 and 39. 


In Eg+ against Inv is plotted (Logarithm of half-life plotted against In Ext). 
for the ft-transitions from 
I= -—1 to J=+1 (figure 3), there is a marked deviation from the otherwise 


smooth line in the region of d=19. (Such a deviation can also be detected in 
the function Fr calculated by Konopinski (1943), which represents essentially 
the same feature.) It follows that the decay of 1%Ne is 2-2 times faster than 
one would expect from its 8+-decay energy. ‘There is also a small displacement 
in the same direction for 1’F, 44Sc and most likely also for 4Na. 

A similar phenomenon occurs for fS- transitions from J= +3 to J= +1. 
In the case of 19O the -decay is.about three times slower than one would expect 
from the extrapolation of the Sargent curve drawn through the points of ?3Ne 
and 27Mg. In either case the nuclei near the minimum of E,, behave as if the 
nucleus with the higher isotopic number had a smaller mass than it actually has. 
The effect reminds us of the similar deviations in the «-decay constants of element 
83 in the Geiger-Nuttall plot, which are caused by a comparatively large change 
of the nuclear radius and density. 


§2. THE INTERDEPENDENCE OF MASS DEFECT, c- AND 
B-STABILITY 


Feenberg (1947) has shown by semi-empirical considerations that a relationship 
exists between the mass-defect curve and the f-stability curve of nuclei, and that 
the irregularities of the one correspond to deviations of the other curve at the 
same values of the mass-number A. This 
follows directly—as was shown independently i 
by the author—from the fact that the sum 
of the a- and. f-energies in a cycle must be 4 
Zero. ; ‘ 

It may be seen from figure 4 that 


4[AE,/A A], const. —2[AE/AT| 4 const. — 0. Be Lae 
= pote (1) Figure 4. and f decay cycle. 
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This can be transformed into 


<a | --|33| /2| == — eA > AEs ; 
AA J}r,=const.=0 ATA AT lie lel ir (2) 


where E, is the binding energy of the last «-particle in the nucleus, and where 
[A7/AA]p,—0 is the gradient of the line of maximum f-stability, i.e. the bottom 


of the Gamow valley. [AE,;/A/] , =[A?M/(AZ)?], is known to be a comparatively 
smooth and very slowly changing function of A for nuclei of odd mass-number A, 
which is always positive. 

Normally E, increases with increasing isotopic number and /°= [Tes=0 


increases with increasing A. 

Equation (2) thus requires that in the region of the well known kinks of the 
B-stability line, where /° decreases with increasing A, the «-binding energy 
E,, must also decrease with increasing isotopic number J. Normally E/ does 
not change very much for neighbouring nuclei with the same isotopic number, 
but in order to fulfil this condition it is necessary that #, undergo a considerable 
change within a series of nuclei of the same value of J. 

The mass-defect D is related to the «-binding energy by the equation 


[AWD NA | conte Fi AD =i), BOS O0-C (3) 


where D, is the mass-defect of 3He. Equation (3) shows that the mass-defect 
curve is approximately linear only when £; is reasonably constant, and that a 
marked fluctuation of E, must cause a similar deviation of the mass-defect 
curve from near-linearity. 

It thus follows from equations (2) and (3) that in the regions where the 
B-stability line shows pronounced kinks, as, for example, near the elements 
igh, ag lc, 5gCe and 5,1, there must be equivalent kinks in the «-stability lines 
and in the mass-defect curves. Vice versa, in such regions where the «-stability 
is known to show marked fluctuations, e.g. near ,,F and ,,B1, similar deviations 
must also be expected for the f-stability line and, obviously, for the mass-defect 


‘CUTIVE. 


§3. THE ISOBARIC CROSS-SECTION THROUGH THE NUCLEAR 

MASS SUREACE IN THE REGION OF (GHE HEAVY ELEMENTS 

The preceding study of the «-energies of the radioactive nuclei of large mass 
has resulted in a diagram (figure 1) which renders possible the interpolation 
of unknown «a-energies with an accuracy of approximately +0-1Mev., and a 
similar, if not better, accuracy for the differences of #, for isobars (nuclei of 
constant mass). This makes it possible t~ calculate also the f-energies (Hz) of 
unknown nuclei from known f-energies by the relation 

gE p= 72k p+ GE Spee 

(E, stands here for the total energy change in a £-transition equal to AM). 

The accuracy of the f-energies so obtained depends on the number S' of 
such steps required to obtain the £, in question, but it should be better than 
10-1, Smev. If it were a purely statistical.error it would be 0-1 ./S Mev. 


This procedure makes it possible to calculate the energies of B-decay (AM) 
of a much larger number of isobars than is obtainable in any other region of the 
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system of nuclei, not excluding the region of the fission products, and to investigate 
how far our assumptions concerning nuclear masses are correct, as the f-energy 
is the gradient of the nuclear mass surface along a line of constant mass-number. 

Of particular interest are the isobars of odd mass number A for which the 
current theories predict a linear variation of AM with the isotopic number J. 

The calculations have been carried out for the largest and the smallest mass- 
numbers of the two types of odd A, the mass differences of which can be obtained 
safely, i.e. 233 and 221 for the series (4n+1) and-239 and 219 for the series 
(4n+3). These two extremes give a very similar picture and, whilst in the one 
case the small values of Z and in the other case the large values of Z show the 
greatest uncertainty, one can be sure that the extrapolation has not caused any 
significant disturbance. 

The results are shown in tables 1 (a) and (8). 

The B-energies used for the extrapolations are (in Mev.) 
for the (42 +1) series : 


eee ile 6, oe 0: 8, con oe Ale 05, ca ile 3, = 0: 7; 
for the (47 + 3) series : 
93 116, °53 0-67, 7S5-~0-2, 745 0-22, 743 1-4, 73] 1-47. 


Table 1(a).. Type 4n+1 
ee A=221 


Z If AM AM Number G i AM AM Number 
of steps of steps 
90 53 84 53 
pe 0 5) 3 down 
91 51 1-0 : 85 51 0-9 
roe 0 24 3 down 
92 49 2:4 86 49 {los 
re 2 up 0-6 1 down 
93 47 0:5 87 47 0-6 
= Pr 5 up ~0 tan 
PEE i 88 45 1-2 
= 2346) 6 up ; es) arp 
95 43 89 43 
Table 1(8). Type 4n+3 
pes A=219 
j 2a, Number ; Nia 
é 4 Me! of steps Z it AM A&M of ster 
92 55 82 55 
ee oan 329 5 down 
oe 0-5 $353 on 
oA v 3°5 5 down 
94 51 1:8 94 51 are 
ae ue 19 3 down 
95 49 0-9 85 49 oe 
a 2p tht 3 down 
96 47 ‘les 86 47 fee 
= Be5 7 up : 0-4 eis 
ak) xe Os S7 ia 0-9 
mas oat 0:5 3 up 


"2 aes 88 43 


eee t= 
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As A*M within the same nuclear type is approximately the same for given 


values of J, the average of the values of A?M are likely to give a more accurate 
picture; these are shown in table 2. 


Table 2. Average A2M 


ik 33 Syl 49 47 45 
A=4n--1 — 0:95 ies 0:55 NODS 
A=4n--3 0-45 ily 0:7 laws C $5 


Near A=229+10. 


It is apparent that A? is by no means constant, as required by the various 
theories (Weizsacker 1935, Wigner 1937, 1939, Feenberg 1947), but that it varies 
in such a way that nuclei with an even value of Z always have a larger value of 
A?M. ‘This means that the mass-defect values of the isobars of odd mass- 
number plotted against J do not lie on a smooth parabola but on two parabolas, 
for odd and even values of Z respectively, the latter values being on the average 
0-2Mev. smaller than the former. The effect is similar to, but not quite so 
marked as, that of the nuclei of even mass-numbers, where the greater stability 
(lower mass) of nuclei with even values of Z is well known. 

The trend in the variations of A?/V seems to be similar for the (47 +1) and the 
(4n +3) types of nuclei, and approximately the same values of A?M occur in the 
case of the (47+3) types at values of J, which are two units higher than those 
for the (47+1) series. It would appear that this effect is another evidence of the 
formation inside the nucleus of “‘complexes”’ consisting of two protons and two 
neutrons which, in the presence of excess neutron pairs, would result in a some- 
what greater stability for nuclei with all protons paired (even Z) than for nuclei 
with all neutrons paired (odd Z). 
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ABSTRACT. Astudy has been made of the results obtained in different parts of the world 
for the effects of ground temperature on cosmic-ray intensity. It is found that the great 
discrepancies between the values for the temperature coefficient as given by different 
observers, as well as the seasonal change of this coefficient, far from being a difficulty in 
Blackett’s interpretation may be taken as further evidence of the instability of the meson. 
The 12-monthly variation in cosmic-ray intensity as obtained by Forbush can be also 
accounted for as en effect of meson decay if one assumes that the bulk of the penetrating 
component originates at 7-5 cm. Hg pressure-level. This assumption has been put to the 
test and the results obtained appear to confirm that the layer of meson formation lies at, or 
little above, the 7-5 cm. pressure-level, at the height of the maximum of the Pfotzer curve for 
the intensity of the total radiation. 


Sue ENR OD Wie On 


N the last fifteen years, it has been reported by a number of investigators 

that the intensity of cosmic rays at constant atmospheric pressure shows 

from day to day and from month to month a variation inverse to that of the 
air temperature near the ground. It is well known that Blackett (1938) attempted 
to relate this temperature effect to the instability of the meson, and that in this 
way he explained the decrease of cosmic-ray intensity in warm weather as due 
to the greater height of the pressure level at which mesons originate. As Blackett 
pointed out, the coefficient of the temperature effect at the equator should be 
lower than at moderate latitudes since the incident primary rays are more 
energetic there. 

However, when the experimental results obtained in different parts of the 
world are closely examined, several difficulties in this interpretation appear 
to arise. We have first the fact that the values of the temperature coefficient 
for stations at moderate latitudes, as obtained by different observers, show 
discrepancies among themselves which are much greater than would have been 
expected from the accuracy of the measurements. These discrepancies are 
particularly noticeable if we compare the values obtained by following one method 
of computation with those obtained*by another. Some observers have correlated 
the monthly or seasonal means of cosmic-ray intensity with monthly or seasonal 
means of ground temperatures. Others have computed the coefficient for 
each month separately by using the daily means and then averaging the twelve 
coefficients. ‘The former method appears to lead to values which are about 
twice as great as those obtained by the latter. 

In the second place, when computed for each month separately, the temperature __ |} 
coefhcient shows a most unexpected seasonal variation. | 

Forbush (1938) has shown that the annual variation in cosmic-ray intensity 
consists of two components—one of a world-wide character related to geo- 
magnetic activity, and the other a 12-monthly wave inverse to that of the 
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temperature. ‘l'here might be added as a further difficulty in the interpretation 
by meson decay the fact that the comparison of this 12-monthly wave with the 
annual change in ground temperature leads to a temperature coefficient which 
varies from one station to another. 

In a former paper (Duperier 1941) the author showed, though in a manner 
which could not be considered as conclusive, that some at least of these difficulties 
would disappear if the mean temperature of the atmosphere, instead of the 
temperature near the ground, had been used for the correlation with cosmic-ray 
intensity. 

The more complete information concerning the upper air, which has been 
gathered in the last few years in England by means of four daily radio-sondes 
from the same station, has made possible a further study of the extent to which 
the temperature effect may be related to the instability of the meson. 


SOs Dobe TEMPERATURE ERPREBCT 

In a previous publication (Duperier 1945) it was reported that the intensity 
of cosmic rays, as measured by a counter arrangement registering triple coinci- 
dences, appears to be closely related to the height of the 7-5 cm. Hg pressure-level. 
It was found that the correlation of cosmic-ray intensity at constant atmospheric 
pressure with the height of lower pressure-levels is rather low, though gradually 
increasing with height. When the 7-5 cm. pressure-level is chosen—the highest 
for which sufficient meteorological data was then available—this correlation 
teaches the values —0-67. By admitting the principle of the instability of the 
meson and assuming that the bulk of the penetrating component originates at 
this pressure-Jevel, the value 18-6km. was obtained for L, the mean range of 
mesons. 

Now, if 64 represents the change in height corresponding to the variation 
dt of ground temperature, the coefficient « of the ground temperature effect 
will be given by the expression 

et LL MOmiOr)s Bs las © 0 bien ee (1) 

As.the heating of the atmosphere by conduction and convection, and so 
the expansion 6A corresponding to a certain variation dt of the ground temperature, 
is the greater the longer the period during which this variation occurs, it is to 
be expected that the value of the ground temperature coefficient obtained by 
using variations from month to month will be greater than that obtained by the 
use of variations from day to day. 

In the last few years it has been the practice in England to make radio-sonde 
observations every day at 00, 06, 12, 18 hours G.m.T. For the present study, 
the upper-air data obtained at Larkhill (about 100 km. to the south-west of London) 
during the period November 1943—October 1945 have been used, and the heights 
of the 7-5 cm. pressure-level determined. 

The mean heights at 00, 06, 12, 18 hours for each season have been plotted 
in figure 1, together with the mean temperature near the ground. ; 

It can be seen that, contrary to general belief, the mean height for each season 
changes quite appreciably with the time of day. As the number of days used 
for the computation of the means (more than 2000 in all) is sufficiently high 
practically to eliminate the effect of change in ground pressure, the variation 
in height of the 7-5 cm. pressure-level has to be ascribed to the thermal structure 
of the air below it. 


3-2 
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To see, however, how closely the diurnal variation in height is related to the | 
diurnal variation of ground temperature alone, departures from daily mean — 
height for each hour have been correlated with similar departures of ground 
temperature. By doing this, the seasonal effect is obviously eliminated, and | 
we find a value of 0:92 +0-02 for the correlation, and 18+1 metres per °c. for | 
b,, the regression coefficient. 

If we substitute 5, for 52/d¢ and 18,600 m. for L in (1), we have for the tem- 
perature coefficient «, = —0-10°% of mean intensity per °c. This value can be | 
compared with those given by Hess et al. (1940) and by Hogg (1947), since these _ 
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observers have followed a method of computation by means of which the seasonal |} 
effect is also eliminated. From series of observations of cosmic-ray ionization 
covering periods of several years, they obtained «, for each month separately 
by using daily means, and then averaging the twelve coefficients. As table 1 
shows, the three values are the same. 


‘Tablevt 
Locality a; (%) Observer 
The Hafelekar (Austria) —0:09 and —0-107 Hess et al. (1940) 
Canberra (Australia) —0-11 Hogg (1947) 
London —0-10 Duperier 
Mean —0-10 


Hess et al. and Hogg have also reported the existence of a marked seasonal 
change of the temperature effect. Hess has found that «, is nearly twice as great 
in winter as in summer. From Hogg’s results it can be seen that the ratio of 
the two coefficients is 1-5. 
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That the change in height of the 7-5cm. pressure-level is more pronounced 
in winter than in summer, compared with that in ground temperature, is already 
apparent in figure 1. From the regression coefficients for the winter and summer 
groups separately, the value 1-6 is obtained for the ratio %1w/%g, which is quite 
consistent with those of Hess and Hogg. 

If we correlate monthly 
averages of the means of the “ 
four heights and temperatures 16480 + 
for each day (plotted in figure 
2) instead of correlating daily 
departures, then we obtain for 
the correlation and regression 
coefficients the values 
9-98 +0-01 and 45+2 metres 
per °c. respectively. By 
applying the equation (1), 
we have now for the tem- 
perature coefticient by = Jan. Mar May July ~ Sept. Nov. 
—45/18,600 = —0-24% per °c. SIS 
which is about double the 
value of «, obtained before. The temperature effect obtained from monthly 
means should therefore be about twice as great as the value from daily means. 

Table 2 gives the coefficients which have been obtained in different parts 
of the world by using monthly or seasonal means of cosmic-ray ionization and 
ground temperature. 


16320 + 


16160 | 


Height of 7-5 cm Hg pressure -level 


16000 | 


y, aunzesadway 


43276 


dtablesZ 
Locality Ga Yo) Observer 
Voyage Vancouver (Canada) —0-18 Compton and Turner (1937) 
to Sydney (Australia) —0-19 Gill (1939) 
Cape Town —0-12 Schonland et al. (1937) 
Amsterdam —0-21 Clay and Bruins (1939) 
Mean —0-18 


It can be seen that all the values are greater than those given in table 1 and 
that the mean is about twice as great, in agreement with the results from the radio- 
sonde observations at Larkhill. The relatively small value for Cape ‘Town, at 
32:7° S. geomagnetic latitude, may be partly due to the primary radiation being 
more energetic there. 

On the hypothesis of the instability of the meson, the use of daily means 
results in a smaller value for « on account of the lag in the warming of the atmo- 
sphere relatively to the warming of the ground, as previously stated. As the 
height of any pressure-level actually depends on the spatial temperature defined 


by i t(z)dz/z of the air below it, it is clear that had the cosmic-ray intensity 


been compared with this temperature, the result would have been the same 
independent of the length of the period over which the data are averaged. 

In conclusion, the discrepancies between the temperature coefficients as 
obtained by different observers are due at least for the most part to the methods 
of computation having been different. Far from being a difficulty in the meson 
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decay interpretation, these discrepancies, together with the seasonal change 
of the ground temperature effect, may be taken as further evidence of the instability 
of the meson. 


§3. THE 12-MONTHLY VARIATION 


Figure 2 shows the annual change in height of the 7-5cm. pressure-level. 
The harmonic analysis leads to a 12-monthly variation with an amplitude of 
307 metres and the maximum in mid-summer. If again we assume that the 
bulk of mesons originate at this pressure-level we should have at the latitude 
of London a 12-monthly variation in cosmic-ray intensity with an amplitude 
of 307/18,600 = 1-65°% of mean intensity, and the maximum in mid-winter. 

Forbush has shown from the analysis of the observations made at five widely 
separated stations that the annual variation of cosmic-ray intensity consists of 
two components—one of a world-wide character corresponding with geomagnetic 
activity, and the other a 12-monthly wave inverse to that of the temperature 
near the ground. Table 3 gives the amplitudes and times of maxima of this 
12-monthly wave as obtained by Forbush for The Hafelekar, Cheltenham and 
Christchurch, the latitudes of which are comparable to that of London. : 


able 3 
12-monthly variation 
Station Latitude Geomag. lat. Ampl. (%) Maximum 
Cheltenham (U.S.A.) Steg INIs 50:1° N. 1-6 . January 1% 
Christchurch (N.Z.) 43-5° S. 48:0° S. 0:8 July 28 
Hafelekar (Austria) 47:°3°N. 48-4° N. 1:9 January 15 
-Mean 1-43 
London (predicted) LS? IN], 54-0° N. 1:65 January 18 


It can be seen that the predicted amplitude and phase for London is in fairly 
good agreement with the mean value of the amplitudes and phases at the other 
three stations.* 

As the upper air observations show that for stations at similar latitudes the 
range of the annual variation in atmospheric temperature, in contrast with the 
range of the variation of the ground temperature, is but little affected by the climate, 
it may be taken that the variation in height of the 7-5 cm. pressure-level in London 
is comparable with the mean of these variations at the other stations. If this 
is so, the agreement shown by the table may be taken as indicating strongly 
that the 12-monthly variation in cosmic-ray intensity is solely the effect of meson 
decay. 

This conclusion seems to rule out the possibility of explaining the observed 
12-monthly variation as due to the heliomagnetic field, as suggested by Vallarta 
and Godart (1939). If this field does not give rise to such a variation, the reason 
is perhaps that the numbers of positive and negative particles in the primary 
radiation are equal or, alternatively, it may be that the properties of the permanent 
magnetic field of the sun are such that this field does not change the intensity 
of the incoming radiation in the course of the year. 


* Yearly variation not determinable by counter-apparatus. 
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§4. THE HEIGHT OF MESON FORMATION 

We have seen in the preceding sections that the temperature effect as well 
as the 12-monthly variation of cosmic-ray intensity, as obtained in different 
parts of the world, are well explained by assuming that the bulk of mesons 
originate at the 7-5 cm. Hg pressure-level. 

In addition to this we have the fact that the value of 18-6km. for the mean 
range, obtained by the same hypothesis, leads to a rest life which agrees with 
the one generally accepted (Nereson and Rossi 1943). If we take for the mass 
180 m,, and for the relevant average momentum of the mesons along their path 
down to sea-level the value 2800 Mev/c. given by Rossi (1939) (derived from 
Blackett’s spectrum and assuming a momentum loss of 2 x 10%ev/c. per gm/cm?), 
then from L=pr /u we have +) =2-04 x 10-8 sec. 

We might therefore conclude from these results that the bulk of mesons is 
generated at the 7-5 cm. pressure-level. 

To see, however, what the results would be if a higher pressure-level had 
been assumed for the layer of meson formation, all the preceding computations 
have been repeated taking the 3 cm. pressure-level, the highest for which sufficient 
meteorological data are now available. ‘Table 4 compares the results so obtained 
with the previous results. 


Table 4 
GA) asco) 12-monthly var. ampl. (°%) 
3:0 cm. Hg pressure-level —0-15 —0-43 2:40 
7-5 cm. Hg pressure-level —0-10 —0-24 1:65 
Observed values —0-10 -—0-18 1-43 


We can see that the hypothesis of meson formation at the 3 cm. pressure-level 
leads to values which are roughly 50°, higher than those for 75cm. Hg. ‘To 
reconcile them with the experimental results, we should clearly have to take for 
the mean range of mesons a value of 28km., a 50% increase, and this value, 
together with the fact that in this case we should have to take for the average 
momentum of the mesons along their path 2100 Mev/c. (Janossy 1948), would 
lead to a rest life of 4:0 x 10-®sec., a 100% increase over the accepted value. 

The conclusion appears therefore to be that the bulk of mesons is generated 
at, or little above, the 7-5 cm. pressure-level, at the mean height of 16km. 

If the average layer of meson formation lies at that height, it follows that the 
intensity of the meson component should be decreasing at 3cm. Hg (22km.). 
On the other hand, the experimental result of Schein, Jesse and Wollan (1941) 
shows that the hard component increases gradually up to the height of 2cm. Hg 
reached in their experiment. The two results, however, are not necessarily 
mutually exclusive if, as is generally supposed, the mesons are formed by primary 
protons the penetrating power of which is not too different from that of the hard 
component. The intensity of the primary protons should increase gradually 
with height. 

In connection with the origin of the soft component, it may be significant 
that the layer of meson formation appears to be located at the height of the 
maximum of the Pfotzer curve for the total intensity (1936). 
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ABSTRACT. Radiation of frequency 200 Mc/s. was transmitted through the covering 
rocks of certain railway tunnels from a radar transmitter feeding a special aerial. The object 
was to measure the attenuation of rocks in situ to see if there might be practical applications 
to geophysical prospecting. The attenuation in dry sandstone was found to be 3 to 4 db. 
per foot, which was much greater than had been hoped. It thus appears that the pene- 
tration is insufficient to be of practical use. 


§1. INTRODUCTION 


HE use of electromagnetic radiation for geophysical prospecting, as opposed 

to the use of induction methods, does not appear to have been extensively 

discussed. A great deal of literature exists, however, on the effect of ibe 
surface conductivity of the ground on the propagation of surface waves. Recently 
in America attention has been focused on the possible use of surface wave measure- 
ments for detecting dykes and faults or other concealed geological discontinuities 
which have a considered effect on the signal strength measured near them. In 
Germany Grosskopf and Vogt (1940, 1942, 1943) have been principally responsible 
for the interpretation of effective surface conductivities in terms of a structure 
consisting of layers of different conductivity. In Russia, Mandelstam and 
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Papalexi (1938) have developed a technique which was used by Al’pert (1945) to 
investigate the refraction effects near a horizontal conductivity boundary. 

By contrast, little attention has been drawn to the possibility of penetrating the 
earth by suitably directed radiation. It has usually been regarded as impossible 
because of the conductivity of the rocks; nevertheless, consultation of available 
figures suggested that some rocks might be found with sufficiently low conductivity. 
Jakosky (1937), dealing quantitatively with the effect of moisture content on the 
conductivities of different rocks, found that the conductivity depended principally 
on the amount of moisture present, and very little on the material of the rock. 
A sandstone with 2-3°, water was found to have a conductivity of about 10-14 
E.M.U.* Figures for effective surface conductivity given by Grosskopf and Vogt 
ranged from 9.10" to 4.1048, and they suggested that the rocks below the upper 
layer might have smaller conductivity. Fritsch (1943) performed some experi- 
ments on radio communication in mines at frequencies between 4 and 30 Mc/s. and 
concluded that the ranges achieved accorded with a smaller attenuation than was 
commonly expected. 

The simplest way of using electromagnetic radiation for prospecting would be 
to measure the time required for a pulse to be transmitted to a buried object-and 
reflected. For such a “radar’’ technique the frequencies used must be high 
enough to enable sufficiently short pulses to be generated. Setting the pulse 
length arbitrarily at 3 usec. it is seen that frequencies much below 200 Mc/s. are 
useless. The U.H.F. conductivity loss is nearly independent of frequency, but the 
dielectric losses, notably of water, increase rapidly with increasing frequency: 
at 200 Mc/s. these losses are still negligible in comparison with the conductivity 
loss, and the losses of other constituents in the main body of rock—mica, 
quartz, calcite, felspar, etc.—are still small compared with the loss in the water 
present. 

A calculation given later shows that, for a conductivity of 10-14 and a frequency 
of 200 Mc/s. reflections might be obtained from within about 50m. Losses due 
to scattering, imperfect reflection, aerial mismatches, etc., would affect this range, 
but since it is determined mainly by an attenuation it would not be reduced 
appreciably provided these losses did not exceed a few db. On the other hand, 
the range would be almost exactly inversely proportional to the conductivity. 
Nevertheless, in view of the lack of data and the importance of the possibilities, an 
experiment at a selected dry site was considered practicable. 

To determine the attenuation without using reflections, a radar transmitter was 
mounted above a railway tunnel and provided with a special aerial for directing 
radiation into the ground after the removal of the surface soil. An exploring 
receiver was used to measure the field in the tunnel, the output being displayed ona 
cathode-ray tube with a high-speed trace. The tunnels investigated were: 
Colwall, Worcs., keuper marl; Alderton, Glos., great oolite; and Pontypridd, 
Glam., coal measure sandstone. 


ce 


§2, EXPERIMENTAL METHOD 


Essentially the method consisted in transmitting a known power downwards 
into the ground above the tunnel and measuring the field strength with a 
receiver in the tunnel after passage of the radiation through a given thickness of 


rock. 


* p.M.U. are used throughout the paper unless otrerwise stated. 
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The field strength E at the receiver is given approximately by 

E(v/cm.) =(7:0\/(PGfafr)/1) x e7 Kt, Bites ee (1) 
P=power of transmitter (watts); fp=fraction of power entering ground; 
fr=fraction of power emerging from rock to reach receiver ; G=aerial gain* ; 
t=thickness of rock (cm.); r=distance from transmitter to receiver (cm.); 
h is defined as the attenuation constant and 1/k is seen to be the distance in which 
the field strength is reduced to 1/e of the value it would have in the absence of the 
rock (k=0 in vacuo). 

It can be shown that 2k? =(472/A?)[4/(e2 + 40?/f?) —e], where A= wavelength ; 
«= dielectric constant ; o = conductivity in E.s.U.; f=frequency. 

At 200 Mc/s. we find, for reasonable values of « and a, that « >20/f and 

likmev/e2ae; 2.) eee (2) 
with «=9 and c=10-“E MW. 
1/k ~ 16 metres. 

The exponential factor in (1) is thus seen to be the governing term in deter- 
mining the-range of penetration to be expected and that, therefore, when solving 
for r, orders of magnitude only need be inserted for quantities in this expression 
other thank. jf, may be assumed 0-5 since we can control the method of direction 
into the ground, but fg may be only of order 0-1 since for practical reasons it is 
impossible to provide a special receiving aerial inside a railway tunnel. Putting 
E=10-§*, PG=10° and r=t, wethen find r~150 metres. This distance will, how- 
ever, be almost exactly inversely proportional toc. Inthe event of the attenuation 
being low enough for reflections to be obtained, the expected maximum depth of 
detection would be rather less than half this figure. It is not expected that e will 
show very wide variations and the expected range depends in any case only on its 
square root. : 

In order to determine k, the variation of E with ¢ should have been studied. 
This would have involved serious practical difficulties. Either the rock would 
have had to be removed from the site to reduce the thickness between transmitter 
and receiver or else new sites would have had to be occupied with a thinner cover 
to the tunnel. The former alternative would have been prohibitively laborious 
and the latter procedure would have been invalid if the rocks were inhomogeneous 
in a horizontal direction; the terrain might also have proved to be unsuitable. 
‘Therefore it was resolved to measure each of the quantities in (1) separately at any 
site which proved suitable. It must be borne in mind that the main aim of the 
experiments was to discover rapidly whether there was any chance of the proposed 
technique being successful: had it so proved, more accurate measurements would 
have been undertaken. 

Measurement of F, t andr was simple, and a theoretical figure was used for G; 
the measurement of P, fp and fp presented difficulties. It is clear that in en- 
deavouring to direct power from a transmitter in air into the ground there will be 
loss due to reflection at the interface. For the purpose of these experiments the 
magnitude of this loss was immaterial so long as it could be measured. Therefore 
a special aerial was constructed (see §3) from which power could not escape except 
either by entering the ground or returning down the lead to the transmitter. 
The power flowing both to and from the aerial was measured and the difference 


% This is the ratio of the power per unit solid angle sent in the direction of the receiver to that 
which would be sent if the power were uniformly radiated in all directions. 


Ultra-high-frequency attenuation by rocks 43 


therefore represented the amount entering the ground. In this way the values of 
P and fy were found. 

fx was simply the transmission coefficient from the rock into the air in the tunnel 
since no special aerial was used. It was assumed to be the same as that at the 
surface where the transmitter was situated, and was deduced from measurements 
made on the transmitting aerial. An overall test was also made in free space by 
transmitting across a stretch of water and comparing the observed value of E with 
that calculated from P, Gandr. This test gave satisfactory agreement. 


§3. DETAILS OF THE EXPERIMENTAL EQUIPMENT 
A block-schematic diagram of the equipment is given in figure 1. 


TRANSMITTER SITUATED ABOVE TUNNEL 


100 kc/s TRANSMITTER 
CRYSTAL 
OSCILLATOR 


100KW 219 Mc/s 


SIYLN3ZONOD — O1—NIML 


M 
AERIAL 25 YDS EASURING LINE 
| FLEXIBLE 20db 
CONCENTRIC OIODE RF 
CABLE 
| OETECTOR | 


PEAK PULSE 
DETECTOR 


RECEIVER SITUATED INSIDE TUNNEL 
DIPOLE AERIAL 


q < By 
F 100 KC/s 
“O. 


CAIN CONTROL DISPLAY FINE ADJUSTMENT * 2 c/s 


Figure 1. Block schematic diagram of apparatus. 


(1) Transmitter 


Type T3154B; frequency 219 Mc/s.; nominal peak power 100kw.; pulse 
length 2 ysec.; recurrence 3-400 c/s. 

It was necessary to measure the peak power being actually transmitted into 
the ground and to monitor this continuously in the field. Since no wattmeter 
operating at this frequency was available, a concentric-line directional coupling 
unit was developed. ‘This type of measuring line has only been described in 
Service reports ; asmall amount of the power passing along the line in one direction 
only can be coupled into a suitable detector. ‘The transmitter output was con- 
verted from balanced twin feeder to concentric line by a suitable transformer 
before connection to the measuring line: from the latter, 25 yds. of flexible cable 
(Uniradio 34) led to the aerial. ‘The measuring line was double-ended; it could 
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thus be reversed, so that the power flowing both to and from the aerial could be 
read. The peak power at the transmitter was usually about 75 kw. 


(1) Transmitter aerial 


This consisted of an open-mouthed waveguide of section 45 x 15in. about 3 ft. 
long and provided with a movable reflecting piston and adjustable probe, connected 
by concentric cable (Uniradio 
34) to the end of the measuring 
line (figure 2). The cable 
attenuation was 1-8db. The 
best dimensions of the probe 
were first found by pointing 
the mouth upwards into free 
space and altering the probe 
length and piston-probe spac- 
ing to achieve the best match. 
As the dimensions affected 
only the matching of the cable 
into the impedance of the 
guide, they were not sub- 
sequently changed. The 
mouth was then rested on ROCK SURFACE 
the ground and the positions Figure 2. Waveguide aerial. 
of both probe and reflector 
adjusted to give the best match. Measurements showed that usually 75% 
of the power available at the probe entered the ground, and this was considered 
adequate. 


THREADED SPROCKETS 
COUPLED BY CHAIN. 


SPRING STRIP, SLOTTED. 


FRAME SUPPORTING PISTON FACE. 


i ADJUSTABLE PROBE. 


THIRD DIMENSION 45° 


FLANCE 6 
oS 


(ii) Receiver 

Type R3170; frequency 219Mc/s.; bandwidth 1-4Mc/s.; noise level 
a LO Aw: 

It was necessary to provide a rapid means of monitoring the receiver noise- 
level in the field. For this purpose a noise-source was constructed using a CV 172 
diode and run from accumulators (B.T.H. report). Measurements showed that 
its output near 219 Mc/s. was 3 db. less than the theoretical value. 


The receiver gain was calibrated by well known methods (Cooper and Freeman 
1943, Varley, et al: 1943). 


(iv) Receiver aerial and display 


The receiver aerial was a simple dipole mounted on a bakelite rod and fed from 
a few yards of flexible cable so that it could be used to explore the field. Its 
matching into the receiver was determined with the aid of the measuring line. It 
was found that at most 10% of the power was reflected. 

In order to display the received pulses without introducing a connection 
between transmitter and receiver down which radiation might leak, and at the same 
time to have a time-base of sufficiently high speed to enable small time delays to be 
detected, it was necessary to control the recurrence of both transmitter and 
receiver from crystals. The transmitter recurrence was derived from a 100-ke/s. 
crystal by repeated division with multivibrators, and the receiver time-base - 
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consisted of a triangular waveform with a forward strioke of 7 psec. directly 
generated from another 100-kc/s. crystal. A small trimmer enabled the display 
to be synchronized with the transmitter to within 75 ¢/s. In this way the time- 
distribution of events could be studied sufficiently exactly to distinguish the 


direct signal from that leaking round the tunnel mouth. 


§4. RESULTS 


The selection of the tunnels was a matter of some difficulty. The main 
consideration was that they should be dry, but they also had to have convenient 
access by road, both for the transmitter in its vehicle above the tunnel and for the 
receiver and its generator, which had to be run in ona trolley from the nearest 
station. ‘Traffic density was also a factor in the selection. 

The dryness of the tunnels could only be judged by inspection of the walls and 
the ground above, but this did. not necessarily mean that water was not present in 
the intervening rock. The experiments were conducted in August 1947, which 
was the driest period experienced for several decades. However, neither in 
Colwall nor Alderton tunnels were any indications obtained of signals arriving 
other than by leakage round the mouth: it can only be supposed that the inter- 
vening rocks were wetter than expected. The cover in both cases was about 40 ft. 
thick. 

In Colwall tunnel sites were tried at both western and eastern ends; the western 
end was rather wet (a spring was later discovered near the tunnel mouth), whilst at 
the drier eastern end the transmitter could be mounted only on the road immedi- 
ately above the line of the tunnel and about 35 yds. back from the mouth, where 
considerable radiation leakage was experienced. This made the field in the tunnel 
very complicated. ‘The traffic-density in this single-track tunnel made thorough 
exploration of the field difficult, but no obvious indications of direct transmission 
were noticed. As the keuper marl at this end is at the foot of the steep Malvern 
Hills it would not be surprising if it were damper than superficially appeared. 

In Alderton tunnel the cover is constant and about 40 ft. thick, and consists of 
very inhomogeneous intercalations of limestone andclay. Although this appeared 
very dry, it is quite possible that much scattering and multiple reflection could 
occur. ‘Two sites were tried immediately above the line of the tunnel, and about 
100 yds. from the eastern end, to see if accidental imperfections were affecting the 
penetration, but no success was obtained at either. 

It is clear that the efficiency of transfer of the radiation into the air in the tunnel 
from the rocks is low. Not only is the transmission coefficient small, but the 
curved tunnel roof acts as a lens tending to spread the beam and reduce the field 
strength. From measurements on the transmitting aerial a transmission coefh- 
cient of 0:2 was obtained, and from this the refractive index of the ground could be 
estimated: it was thus possible to calculate the “lens effect” from the dimensions 
ofthetunnel. Inthe case of Alderton and Pontypridd tunnels the total calculated 
loss was about 25 db. 

At Pontypridd the site was considerably better: the cover increased gradually 
so that at the site 100 yds. inside the northern mouth it was only 33 ft. thick; it 
consisted of massive homogeneous dry sandstone and the tunnel was almost 
disused, so that it was possible to run the receiver up and down the line on a trolley 
and observe and plot the signal strength without waiting for the passage of trains. 
and clearance of smoke. 
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In figure 3 is shown the variation with distance into the tunnel of the signal 
strength, assessed as the height of the largest pulse visible. It is believed that the 
full curve shows the decay of the signals leaking into the mouth, and the dotted 
curve that of signals received directly through the roof. It will be seen that the 
observations show a greater scatter in the region below the transmitter and 
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Figure 3. Variation of signal strength in Pontypridd tunnel. 
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near the mouth: this may be because here there are strong interference fields, 
and small movements of the aerial, such as might easily occur in between runs, can 
cause large changes in signal strength. Nevertheless it is considered that the 
small maximum below the transmitter is significant (the distance of the transmitter 
from the mouth was not known until after the measurements had been completed). 
The aerial was mounted about 3 ft. above the offside rail: in the non-interfering 
zone the strength did not depend in a regular manner on position or orientation 
provided the aerial was not too close to floor or walls. Runs were made with the 
aerial polarized horizontally both perpendicular and parallel to the rails, but no 
outstanding difference was observable. 

The trace was also observed whilst the trolley was run back and forth under the 
transmitter site: a particular echo on the trace was observed to grow and decay 


during this test, which was very convincing. THIS ECHO REACHED A MAXIMUM 
The appearance of the trace is sketched in BENEATH TRANSMITTER 

figure 4. A small pulse is seen before the 

direct signal, which should in this case be the aS 

earliest arrival, but since the time-base runs < 
continuously this may in fact be an echo \ 


arriving on the next trace, 10 psec. later. This 

would be possible since many large reflecting wires 

objects were situated at the required range from ao So eee 
the tunnel mouth. Ignoring this pulse, there- TIME _ INCREASING 

fore, the next pulse seen in figure + is delayed Mecoresd., Appearance cairn 
by approximately 1 psec. on the direct pulse, with signals. 
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or about 350 yds.—a distance of the right order for a path round the tunnel 
mouth. 
The signals were strongly modulated by the approach of trains. 


§5. CALCULATION OF ATTENUATION 

The transmitter power is 7-5.104w. and aerial gain ~5; the power flux at 
lyd. distance is therefore 7-5 .104 (5/47) (91-44)?=3-56w/cm®. This is equi- 
valent to a field strength of 69-4 v/cm. 

The observed signal in the tunnel at 16 yds. range was 10 db. above noise. The 
noise-level of the receiver at the time of the experiment was 4.10-¥ w. 

A field of strength fv/cm. produces power A?f?/3207? w. in a receiver matched 
to a simple dipole. Hence in this case, since A\=135cm., we find f=0-838.10-® 
v/cm. ‘The field expected in free space at l6yds. range is 69-4/16 =4-34 v/cm. 
Hence the attenuation to be accounted for is 4:34/0:838.10-® or 134-3 db. 
Allowing 3 db. loss at the transmitter and 25 db. for the effects where the wave 
enters the tunnel, we find the attenuation of Pennant sandstone to be 3-2 db/ft. 
This is about 20 times greater than that expected assuming o= 10~!£.M.U., so that 
the conductivity must have been about 2 x 10-!E.M.Uv. (resistivity 5000 Q/cm.). 
According to Jakosky this could have been accounted for by the presence of 15° of 
water, assuming Pennant sandstone to be similar-to his reservoir sandstone. It 
seems therefore that even after prolonged surface dryness the rocks near the 


surface may be quite damp. 
§6. CONCLUSION 


The penetration of ultra-high-frequency electromagnetic radiation into rocks is 
insufficient to afford a practical method of geophysical prospecting. This is 
probably because of the great power of rocks to occlude and retain moisture even 
after prolonged drought. In cases where there was appreciable scattering due to 
inhomogeneous ground there would be even less chance of success. 

Increase of power and gain would not increase the range greatly since the atten- 
uation is exponential. Sufhcient penetration could not be obtained by reduction of 
fréquency since pulse widths would then be too great for the observation of 


reflections. 
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The Application of a Variational Method to the Calculation 
of Radio Wave Propagation Curves for an Arbitrary 
Refractive Index Profile in the Atmosphere 


By G. G. MACFARLANE 
Telecommunications Research Establishment, Ministry of Supply, Gt. Malvern 
MS. received 23 October 1947 


ABSTRACT. The problem of determining the field strength produced by a microwave 
transmitter beyond the horizon has been intensively studied in recent years. Charts and 
tables, which are based on approximate solutions of the wave equation obtained by the 
phase-integral method, have been produced (Booker and Walkinshaw 1946), from which 
the field strength can be calculated when the modified refractive index of the air (the 
M-curve) varies with height according to a power law. In many cases, however, the 
M-curve cannot be adequately represented by a power law, and serious analytical diffi- 
culties are encountered if the problem is tackled by the phase-integral method. 

The basic problem is to find the eigenvalues D,, and eigenfunctions Uj, of the wave 
equation d?U,,/ds? + {s-+f (s) +D»} Um=0, in which f(s), the M-anomaly, is a given function, 
which tends to zero as s tends to infinity. In this paper a method, previously described 
by the author (Macfarlane 1947), is applied to determine the first epee D, for the 


class of MM/-anomalies that can be represented by a series of the form, ZA, exp (—a»s). 
=1 


Two numerical examples are given to illustrate the method. ‘They were obtained by 
fitting one exponential and three exponentials in turn to an /W/-anomaly calculated from 
meteorological observations made at Kaikoura, New Zealand. In the first case the 
M-anomaly is represented by a surface duct and in the second case by an elevated duct. 
The exponential attenuation rate, which is determined by the imaginary part of D,, is 
calculated in both cases as a function of the wavelength. 

Finally, the theoretical height—gain function for a wavelength of 3 metres is calculated 
from the /M-curve observed 47 miles off shore from Ashburton, South Island, New Zealand, 
in the late afternoon of 4th November 1946. It is found that there is close agreement 
between theoretical and observed height—gain curves. 


SS UNCER OD WiC a rON 

HE theory of propagation in current use for describing the field beyond 

the horizon of a radio transmitter expresses the field as a sum of normal 

modes of the wave equation. ‘The independent variables of the wave 
equation, which are height and range, are separable, since the atmosphere is 
assumed to be horizontally stratified. ‘This simplification reduces the problem 
of solving the wave equation with two independent variables to that of ae: 
the eigenvalues and normal modes of the one-dimensional wave equation * 


a OF ade oF k®\ y(h) e AG Os = 0. eee eee (1) 


Each mode must satisfy the boundary conditions, that U. U,,(0)=0, and that as 
hoo U,, represents an upgoing wave. — In this equation 


yh) = N*(h)— T= 2 0 Mi) ee (2) 


where N(h) is the modified refractive index of the air and A,,, is the mth eigenvalue. 


* Since the subject matter of this paper is closely related to that of Pekeris (1946) we shall 
adopt a similar notation. 
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A number of different methods have been used to evaluate U,, and A,,. 
Booker and Walkinshaw (1946) applied the phase-integral method, originally 
discussed in this connection by Eckersley and Millington (1938), to find (Uo 
and A,, for the M-curve, 

M hy = Ko (h di) Aste ee 2s ee (3) 
where m=0-5. Here d is the duct width and K,, is the gradient of modified 
refractive index at infinite height. 

The phase-integral method leads to solutions which are only approximately 
correct but which, nevertheless, may be sufficiently accurate for most practical 
purposes; the approximation that the atmosphere is a horizontally stratified 
medium is indeed less satisfactory than the analytical approximation made in 
the phase-integral theory. However, the phase-integral method leads to trouble- 
some mathematical problems when the M-curve has a form more complicated 
than that of equation (3). Even when m=0-2 it was too complicated to evaluate 
analytically and it was found necessary to compute solutions by means of a 
differential analyser. ‘Techniques for evaluating the required solutions of the 
wave equation on the differential analyser were devised by Hartree, Michel 
and Nicolson (1946). 

Another method of evaluating U,, and A,, has been described in a recent 
paper by Pekeris (1946). ‘The M-curve is expressed in the form 

M(h)=K,h+ge~% 

and A,, are found by a perturbation method. The procedure is to 
express U,,, as a linear combination of the eigenfunctions U,,° of all the modes 
in the standard case, in which M(h)=K,h. The coefficients A,,,,, of U,,° and 
the eigenvalues A,, are found as solutions of an infinite system of equations. 
They are evaluated by an iterative procedure, which has been found to be rapidly 
convergent. Pekeris has also shown that a more general class of M-curves can 
be handled by a slight extension of the perturbation method. It includes 


M-curves of the forms 
MHA re EXD (a= 0) 2 6 ae (4.1) 
and VMN Keeh Pel ean” CX). (iC ee 6 | Byockh g (4.2) 


7 


and U 


m 


Since an arbitrary /-curve can be represented by either of those series to any 
desired degree of accuracy, it might be said that the general problem has been 
solved. On the other hand the convergence of the procedure is least rapid when 
the mode is slightly leaky and this is a case of great interest. 

The procedure we propose is based on a variational method in which the 
eigenvalue is found from the solution of a transcendental equation. It is not 
easy to estimate the accuracy of the eigenvalues obtained, but the eigenvalues 
found by it for the power-law M-curve of equation (3) agreed within 0-4% with 
those obtained by Hartree, Michel and Nicolson on the differential analyser 
(Macfarlane 1947). : 

We have developed the method for an M-curve of the form 


M(h) =K,h +80 a > Sn exp ( ao) Rose. (4.3) 


so that it can be applied to any M-curve. If it is desired, an M-curve expressed 
in the form of equation (4.2) can also be handled with the information given in 
this paper. 
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Following Pekeris (1946) we introduce natural units of height 
yrstoledl Haha 
qg=dN?/dh=2K,, x 107° 2s 1 0g cai 

Dy =Anlhigy 
so that equation (1) assumes the form 
@U,,|/ds?+{s+fls)+Dy}Um=9. va eee (6) 


In terms of s the M-curve, equation (4), gives 


f 
J 


fis)=A,+ = A, exp(—235); 2 ee (7) 
n=1 
where A =e (7.1) 
pa LOS ce cag, 
ne en (7.2) 
== 26°25.64) ts Can) re 
and w=2X10-(Rig) See (7.3) 


=3-696-%, Ain cm. 
In the work to be reported here we have confined our attention to the determination 
of the first eigenvalue D, and its eigenfunction U,(s). This has been done 
because, in many cases of importance, the contributions of higher modes to the 
field strength beyond the horizon is negligible. In terms of natural units of 
horizontal distance d, 

Nad. Lhe (8) 
and the field strength at range x and height s, due to a transmitter of power gain 
G and power P at height s, is 


By 
os ane x? E exp(tD,,%)Un(S1)Uimn(S2) see (9.1) 
1207PG 
=! = Pb expliDyx)U(S UG) eves (9.2) 


at sufficiently great range beyond the horizon. In these formulae the eigen- 
function U,, is normalized so that 


[ov * 2(s)de=1, ie, <n (10) 


From equations (7.3) and (8) the numerical distance is expressible in terms 
of the wavelength parameter w as 


«= 0°1423-wtd ed in milesn | eee (11) 
‘The exponential attenuation rate of the first mode is, therefore, 
8-686 xd-1Im (D,) = 1-236 wiIm(D,) db. per mile. ...... (12) 


§2. THE VARIATIONAL METHOD 


The variational equations of the wave equation, equation (6), in which the 
function f(s) is given by ewes (7), are ene 1947) 


Di$)= | (Fy) dy— | yé%dy - SA, _ exp(—a,y)b’dy ......(13.1) 
and ye i MONS oct a te (13.2) 


* We take the “‘ effective ” radius of the earth to be 8 x 10® metres. 
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It should be noticed that the constant term Ay is absorbed into D,. D,+A, 
is then the minimum stationary value of the ratio D/H for functions ¢ that 
satisfy the boundary conditions (i) 4(0)=0; (ii) 4(y)+0 as y>oe, where 
—/3 <arge <0, and that have no zero values on the contour from y=0 to y= ove. 
As explained by Eckersley and Millington (1938) a suitable choice for ¢ is the 
Airy function GAi(cy =j,)¢ || 5) ae (14) 
where —j,= —2-3381 is the first zero. 

Parameter c is determined so that the ratio D/H has a minimum stationary 
value and «=1/c. 

Substitution of function ¢, given by equation (14), into the variational 
equations (13) gives 


iC) peas ake , , 
Figs ~DLP(@~ 2) BAe vee ‘at 
where 
Pe : {Ai'(y J.) }Pdy=0-38321; Q=| {Ai(y—j)dy=0-49170 ....., (16) 
4 0 


F(z) = ; exp ( —2y){Ai(y —j,) dy. 


The stationary values of the ratio D/H occur at the saddle points of the function 
D/H in the complex c-plane; that is, where 


To (Dva 4 ie tee eee, 
5.(3) = O | 2P (c+ =) ae eer a (24/0) =—=\()) 


or Deere C048 A ae (ately © sucteses (17) 


n=1 

The root of equation.(17) lying in the fourth quadrant which gives the smallest 
value for the modulus of D/H, is the required value of c. On substitution of 
this value of ¢ into equation (15) the approximate value of D,+ Ag is obtained. 
It should be noted that the constant term A, affects only the phase of the first 
mode, since, by equation (12), the attenuation rate is dependent solely on the 
imaginary part of Dj. 

The application of the variational method to the calculation of D, for an 
arbitrary -curve involves therefore four steps: 

(i) representation of the given M-curve by an expression of the form 


Mh)=Kh+got+ i e,exp(—eylt) = snes (18) 
n=1 


with r, the number of exponential terms, as small as possible ; 
(ii) calculation of the Laplace transform of the square of the Airy function 
and of its derivative for complex values of the argument , 


P(2)= | exp(—2y Ally A)Pdy seen (19) 


and a) — -| yexp(—ay){Aily —p)Pay; eee (20) 
0 


(iii) solution of the minimal variation equation (17) to find c; 
(iv) calculation of A,+D,=D/H from equation (15) using the values of ¢ 
found from (ii). 
When D, has been found the corresponding eigenfunction U,(s) is obtained 
by numerical integration or by use of a differential analyser. 
4-2, 
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For (i) one method of fitting curves of the form of equation (18) to given 
M-curves depends on putting c,—mc and rewriting equation (18) in the form 


M(h) —K..h —gy= = £n exp ( —nch) = X g,2", where z= exp(—ch). 
n=l n=1 


The range of the variable z extends from 0 to 1. Since K,, and gy can be 
obtained from the asymptotic behaviour of M(h) for large values of h, the problem 
becomes one of fitting a polynomial to a given set of points, in the range 0<2<l. 
Gram polynomials, which have been extensively tabulated by Davis (1935), 
are admirably suited to this purpose. - 

For (ii) the functions F(z) and F’(z) have been calculated to five figures at 
intervals of 0-2 in the real and imaginary parts of z in the range x =0-0(0-2)4-0, 
vy =0-0(0-2)3-2, where z=x+zy. Both real and imaginary parts are given in 
a paper by Hay (1947) together with coupled differences and asymptotic formulae 
for extending the range of the variables. 


§3. NUMERICAL EXAMPLES 


In order to illustrate the technique we shall consider the case of an M-curve 
calculated from meteorological observations made at 1941 hours on 10 October 


1500 


O CALCULATED FROM OBSERVATION OF 
TEMPERATURE AND HUMIDITY 
—90035h 
—-=— o(h) =-0331h +3004 37-6 @ 
—-0035Wn 
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Figure 1, MM/-curve found at 1941 hrs. on 10 October 1944 
at Kaikoura, New Zealand. 


1944 at Kaikoura, New Zealand. This profile of modified refractive index can. 
be represented by the following cubic function of exponentials * 
M(h) =0-0381h + 300 + 37-4¢~00085% 4 47.Qe-0'0070h _5 4.4 g~0'0105h 


* The author is indebted to Mr. J. W. Head for this formula. 
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The original M-curve and the eco, —— ee 
approximation represented by | 
equation (19) are shown in | 
figure 1. In the same figure | 
the approximation using only | 
the first exponential term in 
equation (21) is also shown. | 


In figure 2 the M-excess, ta 
denoted by y and equal to & | 
M(h) —0-0381h —300, is plot- = | | 
ted together with the approxi- “7 | : 


mations represented by one 
and three exponential terms. / 
This shows more clearly than oe ae SA GLE DEENA 
figure 1 the degree of fit | Eee Soe ee 
achieved. An even closer fit or | 


using five exponential terms Og 3 ‘5 5 
is also shown in figure 2, but ji = M— -033ih—300 
we shall not make use of it. Figure 2. M-excess for M-curve of figure 1. 


Figure 1 shows how one 
exponential term can represent a surface duct and how the addition of other 
two exponential terms can give rise to an elevated duct. As it is interesting 
to compare the eigenvalues for these two cases we shall illustrate our technique 
by working out both. 

The first step is to find the constants A, and «, of the series for f(s). ‘These 
are given in terms of the wavelength parameter w by equations (7.1) and (7.2). 
Thus 

Ay= 300w 
Ay= “371-40 a, =0-0918/a Aig ae (3 435 \ (22) 
A,= 42:0w a, =0°1836/w A 5 tee PAD 
A,= —51-1w a, =0-2754/w A,%,= —14-073. | 
We now consider separately the cases of the surface duct and the elevated duct. 


(i) The surface duct 
For the surface duct 
M(h) =0-0381h +300 +.37-4e7 700% (23) 


and A, and «, have the values given in equation (22). 
Equation (17), which determines ¢, is 
. 2+1+4479F'(0-0918/ew)=0. — .seee. (24) 


To solve equation (24) put 


= O-09NSC = ee | eee (25) 
and rewrite it as (0-0918/w)®= —29{1+4-479F(z)}. nae (26) 
Since w is entirely real the problem is to find values of z which make the imaginary 
part of the right-hand side of equation (26) zero. Moreover, since arg = —arge 


and since only values of ¢ are allowed for which 0 <arg c<7/3 we need consider 


only values of z for which —7/3 <arg z<0. It is best to consider trapped modes, 
for which w, c and z are all real, separately from leaky modes, for which c and # 


are complex. 
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Trapped modes.—If we plot w as a function of z, a real variable, we find that 
w decreases from infinity at :=0 to a minimum value of 0-242 at g=0°518; 
then it increases to infinity again where 1+4-479F'(z)=0. For w>0-242 there 
are therefore two real values of z that satisfy equation (24), but the smaller one 
leads to a smaller value of D/H, equation (15), and is the correct one to take. 
As w decreases from oo to 0-242, c decreases from 1-345 to 0-732. This is 
indicated in figure 3, which shows the locus of ¢ in the complex plare. For 
w <(0-242 c and w are complex and the modes are leaky. 

Leaky modes. ‘The cal- 
culation of w and c is more 
complicated than for trapped 
modes. Figure 3 shows how 
c varies from 0-732 at w=0-242 
to exp (27/3) at w=0) “Ehe 
complex values of c are found 
as follows: We know that 
at w=0-242, zs =0-518 and for 
w<0:242, z is complex and 
—7/3<arg z<0. Therefore, 
calculate the realand imaginary 


parts of —23{1+4-479F'(z)} ce: 00 

at the four points (0-4, —0-2) VALUES OF W i 
(0:-5,—0-2) (0-6,—0-2) (0-7, LIMITING POINT 
—0-2). These points are Figure 3. Locus of parameter c in complex plane. 


equally spaced, they lie along 

a line parallel to the real axis, and they enclose the value 0-518 between their real 
parts. Inverse interpolation, as described by Davis (1935), is used to find the real 
part of z, x) say, at which Im[—2s*{1+4-479F’(z)}] is zero. The values of 
Re[ —2°{1 +4-479F'(z)}] are then interpolated to find the value at the pcint 
(x9, — 0-2), which, from equation (26), is equal to (0:0918/w)3. Thus w is evaluated 
and c is calculated from equation (25), viz. c=0-0918/wz. Table 1 gives details. 
of the calculation. 


Table 1. Evaluation of w and ¢ for =x, —0-27 


re F(z) =(13-4-479F (2)) —23(1+4-479F'(2z))y 
Oo y Re Im Re Im Re Im 
0-4 —0:2 —0-34328 —-0-12298 0:5823 0:5508  0:0578 —0-0424 
0:5 =) —0-30043 —0-10163 0:3456 0:-4552 0-0871 —0-0195 
0-6 =(-2 —0-25654 —0-08437 0:1490 0:3779 0-1001 0-0234- 
0-7 —(0-2 021998 —0-07036 Ve Masulsil 0:0863 0:0859: 
0:5512 —0-2 0-0956 0-0000. 


Therefore w=0-2008, c=0-7329 + 0-2659:. 

This procedure is repeated for z =x —0-47, x —0-67 and so on until the entirely 
locus of ¢ in the complex plane has been found. D, is evaluated for each trio 
of values w, c and z by substitution into equation (15), which becomes 


D, + 300w = D/H =0-7793(c? —2/c) -76-06F(z). ....., (27) 
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Next the exponential attenuation rate 1-236w!Im(D,) is calculated for each 
value of w. Finally the wavelength corresponding to each value of w is read 
off figure 4, which is a plot of the function 3-696A-?. 
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Figure 4. Parameter w as function of wavelength. 


Table 2 gives the final results. 


Table 2. The eigenvalues and the exponential rates as functions of the wave- 
length for the dotted M-curve of figure 1 


2 w é 300w+D,=D/H y A (cm.) 

Re Im Re Im Re Im 
0-0 0-0 ore) 1-335 0-0 —oo 0-0 0-0 0-0 
0-2 0-0 0:4117 1-115 0-0 —11-87 0-0 0-0 26:5 
0-3 0-0 0:3057 1-001 0-0 — 8:13 0-0 0-0 42 
0-4 0-0 0:2597 0-884 0-0 — 6:59 0-0 0-0 54 
0-518 0-0 0:2422 0-732 0-0 — 6:04 0-0 0-0 59-6 
0°5512 —0:2 0-2008 0:7329 0:2659 — 4-795 0-100 0-055 79 
0-6280 —0:4 0:1412 0:7365 0:4691 — 3:094 0-456 0-212 134 
0-7180 —0:6 0-1034 0:7284 0-6087 — 2-111 0-835 0-332 213 
0-8067 —0:8 0:0803 0:7144 0-7085 — 1:594 1-148 0-402 311 
0-9680 = {| 2 0:0561 0:6664 0:8262 — 1:198 1-556 0-456 . 533 
1-1139 —1-°6 0:0438 0-6177 0-8758 = SAS iley/7/5 0-459 775 
1-262 ==))-()) 0:0364 0:5694 0-9026 — 1:230 1-896 0-447 1022 

fore) = 65 0-0 0-5000 0-8660 — 1:169 2-025 0:0 ore) 


y =exponential attenuation rate 1:236z Im(D,) 


We shall now consider the better approximation to the experimental values 
of M that is given by three exponential terms. 


(ii) The elevated duct 
For the elevated duct 
HIG = 0-03818 +300 4-37-4e7 © 0? 4.42-e 00 — 1 Te eeee ee. (28) 
In place of equation (26) we now get | 


(0-0918/e0)? = —23{1 +4-479F’(2) -+ 10-061 F’(2z) —18-362F"(3z)}. ws. (29) 
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As before, we have to find values of z that make the imaginary part of the 
right-hand side of equation (29) zero. The procedure is the same as in (i) above 
and the final result is given in table 3. 


Table 3. The eigenvalues and the exponential attenuation rates as functions 
of the wavelength for the continuous M-curve of figure 1 


z w c 300w+ D,=D/H y A (cm.) 
Re Im Re Im Re Im 
0:177 0:0 ora) 0-0 0-0 = 0-0 0-0 0-0 
0:2 0-0 0-8085 0:5077 0-0 —26:28 0-0 0:0 9°8 
0-3 0:0 ~ 0:3663 0:8354 0-0 —11-20 0-0 0-0 32-0 
0-4 0:0 0:2669 0-8600 0-0 — 7:692° 0:0 0:0 51:5 
0:5 0-0 0:2254 00-8146 0-0 — 6-490 0:0 0-0 66-4 
0:6 0-0 0:2089 0:7324 0-0 — 5-672) 0:0 0-0 74-4 
0-625 0:0 0:2084 0-705 0-0 — 5-648 0:0 0-0 74-7 
0:6586 =0-2 0:1748 0:7300 0:2217 — 4-558 0-075 0-039 97:2 
0-723 —0:-4 0:1275 0-762 0-422 — 3:036 0-369 0-163 156 
0:7956 —0-6 0:0945 0-779 0:587 — 2:040 0-751 0-285 244 
0-867 —0:8 0:0749 0-764 0-705 — 1:479 1-081 0-366 347 
0-930 —1:0 0-0626 0°732 0-786 — 1:280 1-332 0-412 440 
1-107 —1:6 0:0430 0-624 0-902 — 1:178 1-742 0-447 790 
oo = ©) 0:0 0-500 0-866 — 1:169 2-025 0-0 oo 


| 


é 


FOR STANDARD ATMOSPHERE 
(M=0-O381K + 300) 


FOR M-CURVE REPRESENTED BY ONE 
4 EXPONENTIAL TERM + CASE (i) 


FOR M-CURVE REPRESENTED BY 
THREE EXPONENTIAL TERMS 
CASE (i 


EXPONENTIAL ATTENUATION RATE IN dbs PER MILE 
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Figure 5. Dependence of exponential attenuation 
rate on wavelength. 


Finally, in figure 5 we have shown the dependence of exponential attenuation 
rate on wavelength for the approximations to the experimental M-curve of 
figure 1 given by the surface and the elevated ducts. 

The curve representing a standard atmosphere is also shown. It is seen 
that the attenuation rate is greater in case (i) than in case (ii). For A10 metres 


Radio wave propagation curves for arbitrary refractive index profile 57 


there is very little difference between all three curves but for \<6 metres the 
attenuation rate for the standard atmosphere is much higher than for the atmo- 
sphere represented by the curves of figure 1. In considering the curves of figure 5 
it should not be forgotten that the total attenuation includes that due to the 
factor « * of equation (9.2) as well as that due to the exponential attenuation 
rate, which is shown. 

It must be emphasized that cases (i) and (ii) have been considered in detail 
primarily to illustrate the power and simplicity of this method for determining 
eigenvalues. In some cases as many as seven exponential terms might be 
required to represent adequately the experimental M-curve, but even so the extra 
labour involved is small. 


§4. COMPARISON OF THEORETICAL AND EXPERIMENTAL 
HEIGHT-GAIN CURVES 


The observed M-curve, shown by dots in figure 6, and the height—gain curve, 
shown in figure 7, were observed 47 miles off shore from Ashburton, South 
Tsland, New Zealand, 
in the late afternoon of 
4 November 1946. 
The radio field was 
produced by a set 
operating on a wave- 
length of 3 metres at 
a height of 35 feet at 
Ashburton. Meteoro- 
logical ascents made 
at different points in 
the path of transmis- 
sion showed that the 


2000 


atmosphere was very ) 

: 329 340 360 380 400 -20 10) 20 40 
nearly horizontally M (h) DECIBELS 
stratified, ee that this Figure 6. M]-curve found 47 Figure 7. Height gain on 
assumption in the miles offshore from Ash- A=3 m. The dots 

burton, South Island, New represent measured 


theory was justified. 
The M-curve re- 
presents an interesting 


Zealand, in late afternoon 
4 November 1946. 


values; the full line is 
obtained from single- 
mode wave theory. 


situation in which a 
duct just fails to occur; M is constant up to about 400 feet and then increases 


linearly. For this M-curve only the first mode should be important at a range 
of 47 miles below the horizon. The situation was therefore simple enough for 
single-mode wave theory to apply. The method of analysis follows the lines 
described above : the M-curve is represented by a linear term modified by two 
exponential terms: M(h)=310+0-0365h +46-4e~ 9008" 227 Sen ee This 
function is plotted in figure 6 together with the experimental points, and it will 
be seen that it represents a satisfactory approximation. ‘The eigenvalue D, is 
then calculated for a range of wavelengths (table 4). age . 

Lastly, the height-gain function is calculated by numerical integration of 
the wave equation using the eigenvalue appropriate to A=3 metres (table 5). 
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Table 4. The eigenvalues and the exponential attenuation rates as functions 
of the wavelength for the /-curve of figure 6 


g w c 310w+D,=D/H y A (cm.) 
Re Im Re Im Re Im 

0-4047 —0:2 0:5796 0:4134 0:2040 —10:933 0-152 0-130 16:1 
0:5685 —0:4 0:2488 0:5710 0:-4018 — 4°505 0-489 0-302 57-3: 
0-7106 —0-6 0:1551 0:6396 0-5401 — 2-666 0-817 0:398 116:3, 
0-8319 —0°8 0:1130 0:6674 0:6418 — 1-893 1-087 0-451 187-1 
0:9360 —1:0 . 0:0895 0:6731 0:7191 — 1°507 1-309 0:484 265-4 
1:0268 —{-2 0:0747 0:6654 0:7776 — 1:301 1-481 0-500 348°1 
1:1077 —1:4 0:0645 0:6507 0-8226 — 1:190 1-615 0:507 433-8. 


1-402 1-387 0-492 300 


| 


y =exponential attenuation rate in db/mile. 


Table 5. ‘Theoretical height-gain function 
h (feet) 159 318 445, 036 764 95455 1080559127254 00 a0) 
20 log |U —6:07 —O-17 2:72 5°87. 7-57 9177s 11-09 12-92 14-05 ieee 


Figure 7 shows a plot of the theoretical and experimental height-gain functions. 
It will be seen that the agreement is good. 


It is also interest- ‘ 
the exponential at- 


tenuation rate varies 

with wavelength. The 
results already given 
in table 4 are plotted 


1-0 


0-9 


in figure 8 with those 
for a standard atmos- 
phere and for the 350- 
foot surface duct of 
§3(i). It will beseen 
that. the attenuation 
rate is small, but not 
zero, even at very short 
wavelengths, which is 
accounted for by the 19 100 1000) 
fact that the M-curve ACM 

fails to show any duct. Kj : 

: igure 8. Dependence of exponential attenuation rate on 


This type of M-curve, wavelength for the M-curves of figures 1 and 2. 
with a vertical tangent 


at ground level and no duct, represents a limiting case which divides M-curves 
into two categories: those with ducts and those without ducts. 


FOR DUCT OF 350 FEET 


EXPONENTIAL ATTENUATION RATE IN DECIBELS PER MILE 
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Radio Shadow Effects Produced in the Atmosphere by 
Inversions 
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ABSTRACT. Reports from radar operators in Australia in 1940 revealed marked increase 
in range of surface objects on certain occasions apparently associated with distinctive 
meteorological conditions. It appeared likely that this phenomenon should be accompanied 
by the existence of ‘“‘ shadow ”’ zones in the atmosphere, in which aircraft would not be 
detected. A theoretical investigation by ray-tracing methods shows that such shadow zones 
should exist above certain inversions. Controlled experiments, in which measurements of 
all the relevant physical conditions were made, show that these shadow zones do exist. The 
position and extent of the zones agrees well with the theory. 


§1. INTRODUCTION 


HORTLY after the installation of coastal radar equipment on the east coast 
of New South Wales, in May 1940, operators began to report the 
occasional detection of surface objects at abnormally great ranges. ‘The 

Commander Fixed Defences, Major-General J. S. Whitelaw, pointed out to 
Dr. D. F. Martyn that this effect coincided usually with the existence of an 
extensive anticyclone system over the eastern part of the Australian continent. 
Theoretical and experimental investigation was initiated, at first by the Radio- 
physics Laboratory of the Council for Scientific and Industrial Research, and later 
by the Australian Army Operational Research Group (O.R.G.). 

It was early realized that the phenomenon was mainly due to the sharp 
temperature and humidity gradients associated with certain meteorological 
conditions, notably temperature inversions (Englund, Crawford and Mumford 
1938), and attention was concentrated on this aspect Information soon reached 
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Australia that a similar investigation in England, initiated principally by Booker 
and Walkinshaw (1946), was directed mainly towards the determination of the 
increased fields close to the surface of the earth, in the ‘radio duct”? formed 
within and below such gradients by the “trapping” of radiation. Now it seemed 
probable that under abnormal conditions of high surface fields there should 
occur a weakening of the field above the abnormal gradients. If sufficiently 
severe this might explain the occasional failure of aircraft warning radar sets 
to detect aircraft flying at certain heights when within the normal range of the 
set (Martyn 1943). Accordingly attention was devoted in Australia to this 
latter aspect, which was of considerable operational importance. 

In the early stages of the investigation an attempt was made to derive 
information from the log books of radar operators and aircraft pilots. However, 
this method proved unreliable and a series of carefully controlled flights off 
the Sydney coast was made during the summer of 1942-3. Steep climbs were 
made by an aircraft at a range of about 24 miles from a 200 Mc/s. transmitter, 
250 feet above mean sea level. Field strength and temperature were measured 
in the plane. Humidity was not measured, but use was made of hygrometer 
readings taken by the R.A.A.F. during their regular meteorological flights at 
Richmond, about 25 miles inland. At the same time propagation through a 
horizontally stratified medium was investigated theoretically by a geometrical 
optics treatment, this being considered best for the region near the level of 
inversion (Price 1943). ‘The experimental results were published in a service 
report by D. R. Brown (1943). Examination of the results showed that on days 
of temperature inversion the three lowest lobes of radiation were depressed 
from their normal position, while higher lobes were not measurably affected. 
Analysis by ray methods showed that the extent of the lowering was of the right 
order (Price 1943). The results indicated that further experimental work was 
desirablé, especially along the lines of combined air sounding and field strength 
measurements. Work was commenced on the development of accurate instru- 
ments for this purpose, and taken over by the Radiophysics Laboratory after 
the Operational Research Group was disbanded in 1944. 

Upon the completion of sounding equipment designed for use in aircraft, 
a series of flights was carried out off the N.S.W. coast near Sydney in January 
1946 by the Radiophysics Laboratory. The results proved of interest in connec- 
tion with (a) the development of inclined stratified air masses produced by 
off-shore streaming of hot dry air, modified to some extent by a sea-breeze type 
of circulation, and (6) the sharp decrease of field strength with height experienced 
on passing through a layer with negative refractive index gradient. 

This paper discusses the ray theory of propagation of ultra-high-frequency 
radio waves through such layers. ‘The general theory for inclined layers developed 
in §2 is applied in $§3 and 4 to an analysis of the field above the region of surface 
trapping caused by a negative refractive index gradient. Rays have been traced 
for one of the January flights and the results are compared in §5 with radar 
observations of echo strength from the aircraft. 


§2. RAY THEORY FOR A STRATIFIED ATMOSPHERE 
2.1. Scope of ray theory 


An essential feature of ray theory is that the radiation field can be described 
in terms of the propagation of energy along rays whose direction is given by 
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the laws of refraction, a given parcel of electromagnetic flux travelling within 
a tubular space bounded by the rays. Surfaces of equal phase, orthogonal to 
the rays, constitute the wave front. This description will be valid, provided 
that the laws of the electromagnetic field (i.e. Maxwell’s equations) are satisfied. 
It is known that there are limitations to the validity of ray theory, the main ones 
being :— 

(i) |grady| <1/A, where A is the wavelength. This means that the radius 

of curvature of the ray must be sufficiently greater than J. 
(it) The rays must not come to a focus. 


(iii) The region considered is not near the boundary rays produced by 
obstacles, nor within confined spaces whose dimensions are comparable 
with the wavelength. 

To some extent ray theory can be extended by the use of Huyghens’ Principle 
and the idea of interference. When this method becomes too complex other 
solutions of the wave equation are sometimes available which then enable the 
field to be described in simpler terms, such as, for example, by means of transverse 
electric (or magnetic) modes. Which is used depends on the nature of the 
problem and the completeness required in the solution. Treatment by ray 
methods, keeping their limitations in mind, sometimes provides a useful if 
approximate solution when a more complete quantitative answer is arithmetically, 
if not algebraically, impracticable. ‘This appears to be the case when dealing 
with the propagation of U.u.F. radio waves over the earth, just above a refractive 
index inversion, in the region where the lobe pattern begins to be sensibly 
affected by super-refraction. 


2.2. Equation of the ray 


Figure 1 shows a ray TR from a transmitter T, at height hp above the surface 
of the earth, which is assumed spherical with radius a, and centre.O. The 
position of any point R is given by polar coordinates (a+h, #) or by 
height # and reduced range d a 
measured along the surface of | 
the earth, | lhe contouts ef 
equal refractive index, ~=con- 
stant, may be curved with respect =o 
to the earth and inclined at small 
angle to the horizontal. For 
simplicity, cylindrical symmetry 
- about the vertical through the 
transmitter is assumed, so that 
the path of a ray lies in a vertical 
plane. In general, this would 
not be strictly accurate and 
lateral deviations (which could 
be important in bearing measure- 
ments) are possible, but since all 


s NY 
Vy XS HORIZONTAL, 


» 
“4 = CONSTANT 
RAY 


At ae 8 
deviations are small, it is con- 
sidered that vertical and horiz- re) 
ontal refraction can be treated Fiptire iete corde ete terete hae tee? ford Fay 


independently. passing through inclined layers. 
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We define the following symbols (figure 1):— 
«=angle made by a ray to the horizontal at the transmitter =8 +0 =y + pb. 


B=angle made by a ray to the horizontal. 

y =angle made by a ray to the stratification (angle of elevation). 
is = angle made by the stratification to the horizontal at the transmitter. 

§ =angle made by the horizontal to the horizontal at the transmitter. 

y= radius of curvature of the stratification, w = constant (figure 2). 

q=radius of curvature of the normal to the stratification in the plane of the ray. 
p = radius of curvature of the ray (figure 2). 

s = distance measured along the ray. 

Clockwise rotations are positive, so that downward curvatures are positive, 
in accordance with the usual convention for the earth’s curvature. 

Consider a ray passing from a point P, 
where the refractive index is p, to R, where 
itis w+dyu, while the complement cf the 
angle of incidence changes from y to y+dy 
(figure 2). The deviation of the ray is 
ba = dy + ow. 

As the ray passes from P to R the angle 
of incidence changes continuously, part of 
the change being due to the rotation of the 
contours. ‘To the first order the effect of 
this part on the deviation of the ray can be 
separated from the effect of the alteration 
in y due to refraction, ~vhich is taken care 
of in the formula by the changing value of 
y4, as in the case of parallel layers. Let the 
effective average change in the elevation Figure 2. Deviation of an element of 
angle be —8ys/2 so that the average value a ray produced by a thin stratum. 
of this angle is (y —Sy/2). 

Then by Snell’s law 

Su _ cos(y —54/2) —cos (y —8yh/2 + da) 
be cos (y —8y/2 + 8a) 
— sin (y —8y, 2 + 8x/2)da \ 
a cosy | 


\ZON 
g_HOR 


_sin fy +8y/2 ~@B/2 —Bp/2)}8a 
cosy 
If the interval 5s be chosen small enough for each contour curvature to be 


regarded as constant and for the change in curvature between the contours to 
be small, (Sy —8d:h) will be small and may be neglected. 

Referring again to figure 2, 3s is the element of the ray between the stratum 
boundaries PQ, SR, which have radii of curvature 7,, r, and centres L,, L, 
respectively. 5,7, 5, are the respective lengths of the normals Q’R, PS’, The 
dotted arc UR is drawn concentric with the first surface PQ. 

dy, the deviation of the ray relative to the strata can be considered to consist 
‘of two parts: dyy, the rotation of the ray relative to the first surface, in going 
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from P to R; dyy, equal to the angular displacement at R of the first surface 
relative to the second. ys may be similarly subdivided. Then 

oy =dyy + dy, dh = Sy + dep, Oey = ~-Soif, du=dy +dp=dyy +dyyp. ...... (2) 


5% is due to the displacement of the centre of curvature. 
Equation (1) may be written 


ba. 


We PS aoc © AR ae 
Oe oe (y+ dyy/2) at {sin (y + dy/2) —sin (y + 8,y/2)} 


7 


| ~& 
nS 
a | 8 


Cos y(du/ + doy . 5x/2) 
}yn 


1.e. 


: _ led 
Sia (y+ 8yy/2) = —1/p. ee cer, (3) 
1 


since cos (y+ dy/4 + 8yy/4) ~ cosy, 
1/p is here the average curvature of the ray PR. 
Equation (3) may be rearranged in the following form, suitable for com- 
putation: 
bu 1 Py Oo 
eee pon ¢ — COS (Soi /2)/sin (y + sm) ae) 5yn 
From the geometry of PQ’R, to the first order of small quantities, 
cosy/r,= —sin (y+ 5 ;y/2)ob/O42, na eee (5) 
remembering that 57 is positive when directed upwards and 7, is positive when 
the surface is concave downward. From equations (4) and (5), noting that y, 
can be put equal to y, we get the following relations: 


Ope 1 1 3 Sy 
Ene 7, NO) 07 eR UTZ) —- = Sy) ee 
i Le (; — COS 73 (doi/2)/sin(yy + eo) .: 4 suai") }4n 
a ae eel Se ro (6) 
=cosy,(1/r,—1/cosy;p)=cosyy/hy, ssa (7) 


which is the average curvature of stratification at P relative to a ray inclined 
to it at y. 1/k,* can be called the average “effective curvature”’ (Schelleng, 
Burrows and Ferrell 1933) of the stratification at P. 

When y <1 and p—1 we get 


byn 1 dyn i = 2 
50 bu (a) te = Fri = DUP), 


(y,+8sy/2) is the average slope of the ray relative to the first surface and is 
approximately equal to —8,n/8,d.  (825/2)(y1 + 91y/2) can always be made small 
compared with unity by decreasing 5p sufficiently. 

Since it can be shown that 


Ram =o /oy eee (9) 
from (5), (6) and (7), it follows that the range measured along the stratification 
at P will be given by 

Od =TOyb= —kyoyys are . (10) 


* Sometimes the symbol & is used for the ratio (“ effective radius ’’/actual earth’s radius), instead 
of for effective radius, as here. 
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The optical distance is 
(14 + 8/2) cos (7, + BSyy/2)81d — (+ 8/2) sin (yy + 8yy/2Z)0yN ee (11) 

and so can be separated into two components, one depending upon the distance 
along the stratum and the other along the normal. —9,i/3,n, which is a measure 
of the curvature 1/g of RQ, the line orthogonal to the strata, and the curvatures 
1/r,, 1/p, must be finite or zero for all values of y as du, da etc. approach zero, 
unless there is a discontinuity in the slope of the contours. Thus, in the limit, 
equations (6) and (7) reduce to 


COSY = cosy (HE. + 2) =siny 2 =siny (2 -*) veeee(12) TT 

when du>0. Also 
lR=U/r+dplndn)\=(/7=VWer), eee (13) 
dd=rdibt=—Rdy, ~~ 9) eee (14) 


In the limit the effective curvature of the stratification, 1/k, becomes equal 
to the curvature of this surface, relative to the tangent ray at P, whose curvature is | 


Llpr= —(1/q)(duidn). 
Numerical computation of the ray path. When y is small equations (8) and (10), 
giving finite differences, can be used for tracing the path of the ray. Ifa numerical 
procedure is adopted the path of the ray is determined in steps between con- 
secutive surfaces j.=constant, spaced so that the curvatures of each surface and 
of the ray are approximately constant, and the changes in their slopes (d%, dx) 
are small. 6,7, 1/71, 1/re, y, and {8% —(8d/r, —dd/r,)}, quantities which occur 
at the beginning of a step, can be read at each stage from a large scale plot of the 
contours, similar to figure 8(a). The last quantity is equal to the initial difference 
in % between the two surfaces. An approximate value of 5d can be obtained by 
inspection of the previous path of the ray. 
Values of 5,7, 5,% and (y+6,y/2) are then calculated from the relation 


(y +8,y/2) = —Syn/Sd + {8% —(8d/r, —8d/r,) + 4(dd/r. —Sd/r,)} = —8,n/Sd. 
eee (15) 
Using these values, (8) is solved for d,y by successive approximation. With 
the new value of (y+y,y/2), (15) is solved for dd leading to new values of 8,5 
and 6,. ‘The process is repeated in (8) and (15) until consistent values of 5,y 
and 6,d are obtained. Equation (10) should then check automatically. 
It is convenient to use the modified refractive index described in § 2.3.2. 


2.3. Transformation of coordinates 


The form of expressions (12), (13) and (14) shows that the relative curvature 
of the ray, at any point, viz. —(cosy)/R, is derivable from a “‘ modified” refractive 
index N whose logarithmic gradient dN/Ndn is equal to 1/k, the stratum at the 
point being regarded as flat. We may also regard the process as one of sub- 
tracting a curvature (cosy)/r from any ray making an angle y to the stratification. 
Inclination of the contours due to change in the position of the centre of curvature 
is not thereby eliminated. 


2.3.1. Concentric stratification. When the strata are concentric (13) becomes 
1/k=dyu|udr +1/r=tanydy/dr, 
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since now dy=d,y; dyy=0; din=dr. Integrating, we get the well known 
equation of the ray 


prcosy=constant. = 9 eee. CL7) 
2.3.2. Horizontal stratification. Further, if the stratification is horizontal 
1/k=du/(udh)+1/(a+h)=dN|(Ndh)=tanydy/dh ...... (18) 
=d(y*/2)/dh, <1, 
and dd = adys= —akdy|(a+h), 
where d is here the range, reduced to the earth’s surface. 
Integrating (18) with N=,.* when h=0 we get 
N-1=p(ath)/a-1=(u—-1)+h/a+(w—-1)h/a_ ...... (19) 
and Weosy=constant, «= cee (20) 


The following approximation to the “modified” index N is generally used 
(Freehafer 1943): 


Ni (ted ea we eae Sees (21) 
Also it is customary to define a quantity M, of convenient size, by 
10°M=(u-1)+h/a. ae (22) 


2.3.3. Plot in rectangular coordinates (Booker and Walkinshaw 1946). 
It follows from the last section that, in a horizontally stratified medium situated 
over a flat earth and having a refractive index N given by (18) or (19), a ray will 
have the same relation between y, / and reduced range d, as in a horizontal curved 
medium, refractive index w, over a curved earth of radius a. Thus a plot in 
rectangular coordinates of reduced range against height will give the same ray 
paths as would be obtained by straightening the strata and making their length 
equal to their projection to ground level. ‘These will not be quite the same 
curves as obtained using N over a plane earth, unless the ranges in the latter 
case are also reduced by the factor a/(a+h). 

It is convenient to eliminate the curvature of the earth and use a modified 
index of refraction even in the case of non-horizontal stratification. When 
h/a is now added to pu to obtain N, the lines of constant N will not remain parallel 
to those of constant p, but become inclined at an angle <, given by 

(—dp/dyn)sine=(1/a)sin(/—O+6¢), 4.2... (23) 
this being the condition for zero gradient of N in the direction inclined at € to 
the lines » = constant (figure 3). 

It will be shown that the path Ess 
of the ray will be approximately +| 
the same when N is taken as the 
refractive index and the earth’s 
curvature disregarded as when the 
actual refractive index and curva- 
ture are used, provided the angles 
of inclination are small. 

In figure 3 the angles and radit 
are defined as before, dashes dis- 
tinguishing those measured to 
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Figure 3. Transformation from contours of 
»=constant to those of N=constant. 


* The unnormalized value of N is here defined. 
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N=constant. Cosines of small angles can be put equal to 1. We have from 
equation (1) 


du du. cosy’ ae, 
—— =f —— |) (ES ei omer 
dn an? din (v ) cos y iE dn 
since angles RSP, RS’P each equals 90°. 
R d s (ob —8 cos (ys —8+ cosy’ eae 
Therefore & cose + Sy ~) oly eke cose-—— tany ie 
d\n a a cos y d\n 
. J ; 
But 1/ry= —tan po =¢curvature OPS 5 
1 


on putting dx=d(s+e+y’) and adding the last two equations we get 
(5 cost) a _ cos(# —8+¢) dy 
a 


Up 


1 
—— cose + — — —4__— ] tany’ —,;| 
d\n a a ) r dn’ 


or dN|djn'+1/r=tanydyjdi a) eee (24) 


since ‘1/r’ =1/ry —cos(%* —0+.«)/a=curvature of the N stratum relative to the 
earth and ((du/d,n) cose +cos(—O+e)/a) is the resultant gradient of N. 
Equation (24) gives the same form for the variation of y with height as when N 
is the refractive index and the earth is flat. 


$3, CHARACTERISTIC. OF IN AND VE CURVE SFE OR 
HORIZONTAL STRATIFICATIONS 


An advantage of the use of a modified index lies in showing quickly how the 
inclination of a ray to the horizontal changes with height. When the gradient 
of N(or M)is positive, y continually decreases with height and the ray is convex 
upwards relative to the earth; conversely for a negative gradient. The usual 
condition of the atmosphere is for grad N to have a positive value lying between 
30 x 10-°/ft. and 45 x 10~9/ft. and increasing with height. The standard value 
of 1/k generally used at low levels is 36x 10-°/ft. Temperature inversions 
and, more especially, abnormal hydrolapses producing negative gradients of 
do occur quite frequently, but they can only persist over limited intervals of height. 
The resulting structure of the atmosphere shows a series of maxima and minima 
of N with height, the greatest and most important N deficits occurring below the 
first one or two thousand feet. Higher level N inversions which occur can 
reflect radiation to an extent which is sometimes significant, the amount depending 
on the shape of the (N, /) curve at the inversion, and on the propagation frequency. 
If only the lower inversions are considered, it is found that the (N, h) curves fall 
into several definite classes as 
follows :— 


(a) Figure 4. The gradient 
of N is always pcsitive. -All rays 
are concave up, relative to the 
earth. Sub-standard propagation 
occurs when the value of the 
gradient of N is greater than M 
normal (36 x Cait.) Figure 4. Paths of a ray when gradient M is positive. 
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(b) Figure 5. The initial gradient is negative and a minimum occurs at B. 
Rays are concave down or up, according as they are below or above B. 
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Figure 5. Shadow regions caused by a surface duct. 


From a transmitter below B (figure 5 (a)) there will be a pair of rays leaving 
at angles +y,, which just reach the level of the minimum of N, that at +y, 
doing so after reflection from the earth. ‘The flux emitted between these rays 
is trapped between the earth and the horizontal through B, this region being called 
a surface duct (Booker and Walkinshaw 1946). 

The upper critical ray, initial angle —y,, reaches the level of B at a finite 
distance, whose angular range 0, is given approximately by 


0 
6, =| —kdy|(a+h)= —ky,/a; (R constant, h<a), ...... (25) 
= 


Yo= +V/2(Nr —N5)}. 0 eee Serta (26) 
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Then it returns to earth in a half-hop, is reflected, and continues between earth 
and level B in full hops around the earth. The lower critical ray, +e, reaches. 
the earth at an angular range 8, given by 


Za ie _kdy|(a+hVS kh —y)las (e constet <a) ee (27) 
Ye 
yas t+r*/2(Na—Ng) we wees (28) 
The angular range of the critical ray hop is 
61, = 2(0, + 02) ox —2ky,/a; (Rconstant,h<a). —.... (29) 


Each critical ray may be considered to divide at P, R, one branch continuing 
within the duct as described and the other proceeding into the region of positive k, 
above B. 

The angles of elevation at P, R are zero. Between the two upper branches 
is a laminar region of ‘‘semi-shadow”’ into which the direct rays do not penetrate, 
but which is traversed by reflected rays leaving the transmitter at angles greater 
than y,. The tilt of the lamina increases with height, but for a given height 
of transmitter the angular width of the shadow remains constant at 265. 

It is apparent that this width will be zero when the transmitter is at the surface 
but that the amount of trapped radiation is then a maximum, since y, is a maximum. 
As the transmitter is raised to the top of the duct (B), the trapped radiation 
decreases to zero, the width of the shadow increases to 0, but the angular distance 
to the centre, measured at the level B, remains constant at 04/2 (figure 5 (4)). 

The lower boundary of the semishadow PQR lies below the duct, the area 
PQR increasing with the height of the transmitter till the top of the duct is 
reached. 

If the transmitter is raised above the duct, trapping does not occur but the 
shadow zone is stillformed. ‘The lamina angular thickness now remains constant 
at 64, but its centre moves away from the transmitter with increase of altitude 
(figure 5 (c)). 

Beyond the lamina and above the duct is a full shadow zone not reached by 
direct or reflected radiation. 

The above description presents the ideal distribution given by ray theory. 
As stated in §2.1 the treatment can be made more quantitative to limited ranges. 
determined by the practical consideration of the maximum number of intersecting 
rays which can be handled by interference methods. At long ranges, within 
the duct, use has been made of the theory of guided-waves, the minimum range 
to which the field may be calculated being determined in this case by the maximum 
number of modes which can conveniently be dealt with. Considerable work has. 
been done in this way on surface trapping. Intermediate ranges are not amenable 
to either method and little work by way of calculation has been done for the 
semishadow region. 

It is clear that some radiation will be diffracted into the ‘‘semishadow”’ and 
“shadow” across the curved surfaces PU and RS, and that some will penetrate 
the boundary PQR, in accordance with wave theory, but the diffracted field 
will generally be weak compared with the radiation which would reach the space 
in the absence of trapping. Neglecting then the diffracted field in order to obtain 
an upper limit to the decrease in intensity, it is seen that the normal interference 
field at the maxima of the lobes will be decreased about 6db. in the semishadow 
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region, due to the absence of the direct ray. (This does not take into account 
increased divergence of the reflected ray due to the negative gradient.) Such 
a decrease would be important when searching for aircraft at long and medium 
ranges, since a radar set is then working near the limit of its ordinary detection 
range. 

Ray treatment predicts a second major decrease in field strength on passing 
into the full shadow across the false horizon RS—theoretically to zero. The 
problem of calculating the residual field beyond RS is similar to that for the case 
of diffraction around the earth’s surface near the horizon when superrefraction 
is absent, but is complicated by leakage across the diffracting boundary RV. 
A full wave theory is required. 

The (M, h) curves upon which the discussion has been based are drawn 
with discontinuities of gradient at the maxima and minima, and with constant 
gradient between. Actual soundings show some rounding of the turning points 
and curvature of the joining sections, but as will be seen from figure 8 (d), this 
is often small. ‘The effect of rounding of the minima is to make grad N-0, 
and doo. ‘The critical rays tend asymptotically to approach the horizontal 
through B, but the ray for which the value of y at level T is slightly greater than y, 
cuts the top of the duct at a finite range. ‘he region bounded by the two rays 
at +y will only contain a small amount of direct radiation above the duct, 
proportional to y —y,, and may be regarded as the semishadow zone (cf. figure 8 (c)). 
One result is that the shadow boundaries will be diffuse instead of sharp, a con- 
dition which is aided by diffraction. 
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Figure 6. Surface duct. Elevated layer. ‘Transmitter at max. M, 


(c) Figure 6. The initial gradient is positive and is followed by the duct- 
forming negative gradient with a minimum WN less than the surface value, The 
rays below D will be concave up, and maximum trapping occurs with the trans- 
mitter at D. Conditions otherwise vary much as in (b). 


(d) Figure 7. As in (c), a positive gradient precedes a negative gradient, 
but the minimum N now is greater than the surface value. This is an elevated 
duct. 

Here the lower critical ray descends only to the level A, where N has the same 
value as at B, so that the duct lies between the horizontals A and B. Interesting 
features are (i) the caustic passing through R, extending down to the level of the 
transmitter and up into the shadow; and (ii) the extension of semishadow 
beyond R, the full shadow boundary now being the continuation of the ray 
tangent to the earth beyond V. The rays whose intersections form the caustic 
lie between this tangent ray and the ray which leaves the transmitter horizontally. 
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When the transmitter is below the duct all radiation passes through it, so that 
there is no trapping or shadow. Ass the transmitter moves up through the duct, 
zero trapping occurs at A as well as at B, maximum trapping at D. The sequence 
of variation of shadow width and centre are similar to (b) and (c). 

h 


2(0,+92)2 
ao 


(6) ‘Transmitter above duct. 
Figure 7. Shadow caused by an elevated duct. 


§4. INCLINED LAYERS 
When the stratification is inclined, but concentric, the condition for possible 


trapping is that (N+ far'ir’) must have a negative gradient. We now have 


ve =V/{ —25(N + [dr'/r’)}, 
the increment 6 being taken from transmitter level to the top of the duct. 
If the position of the centre of curvature is changing, the negative gradient 
of (N + {dyn r’) must produce a change in elevation angle 5,y’ greater than the 


rotation d,y’ due to the shift in curvature, in order that trapping may occur. 
In this case there is the possibility that rays trapped at one section may later 
penetrate the trapping layer owing to a change in its tilt. Again, rays initially 
above the minimum of (N + | dn’,/r’) may be trapped on account of a favourable 


change in the tilt. 


§5. EXPERIMENTAL FLIGHTS IN THE SUMMER OF 1946 
5.1. Locality and equipment 


In January 1946, a series of meteorological soundings were made in com- 
bination with radar measurements on a 200-Mc/s. C.O.L. Mark V air warning 
set at the radar station at Collaroy, 400 ft. above sea level and 15 miles north of 
the Sydney Heads. The aircraft flew about 100 miles out to sea at right angles 
to the coast and back, in a sequence of ascents and descents, to heights of 2000 to 
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3000 feet. On three of these flights observations were made on the radar set 
of the time, and of the range, bearing and signal-to-noise ratio of the echo from 
the plane. Observations of range were made every 10 seconds and of bearing 
every 2-3 minutes. The power output was monitored but an absolute cali- 
bration was not attempted. 

In the aircraft a simultaneous record was made of time, number of the reading, 
air speed by air speed indicator, height, and hence pressure from an altimeter, 
and wet and dry wick temperatures, using thermocouples. The indicating 
meters were mounted on one panel and photographed at controllable intervals 
(3 to 6 seconds) by a camera with motor-driven shutter. The thermocouple 
and camera were designed and built by M. Iliffe and will be described elsewhere 
(Iliffe 1948). The hygrometer was mounted forward of the cabin and above 
the fuselage clear of slipstream disturbance and engine contamination. ‘Times 
were synchronized with the radar station. Air speed and rates of climb were 
maintained as constant as possible by the pilot. 

The flights were made on days forecast as favourable to strong M-inversions. 
At this time of the year the system of eastward moving anticyclones, which 
form the central motif of the Australian weather pattern, has a general track 
which crosses the eastern coast about latitude 40° S. and reaches the latitude 
of Sydney in the central Tasman Sea. The circulation associated with the rear 
portion of the anticyclone after it crosses the coast, creates strong off-shore 
streaming of very hot dry air. Modification by the cooler sea from below tends 
to produce the familiar surface advection duct with strong and deep temperature 
and humidity inversions. A further diurnal modification is produced on the 
circulation by the hot land—cool sea influence on the density and pressure gradients. 
This results in a large-scale sea breeze effect which produces a change in the shape 
of the temperature and humidity contours. The reversal of wind direction 
near the sea allows the surface layer of modified air to become cooler and moister, 
and more thorough mixing is achieved to greater height and closer inshore. 
The transition layer between modified and unmodified air becomes thinner 
and higher and an elevated duct is formed. The contours of constant temperature 
and humidity slope upwards out to sea, with a gradient which increases slightly 
from the land, producing an upward curvature. A further effect of the reversal 
of low level circulation is to slow down the eastward movement of the anticyclone, 
causing the inversion conditions to persist for some days. 

The detailed history of the development of these “coastal fronts”’ (Martyn 
and Squires 1944) is of considerable interest, and this aspect is to be investigated 
further by the Radiophysics Laboratory. In this paper, however, the interest 
is in the illustration afforded of the decrease in field strength across the transition 
layer between the lower and upper air masses. ‘The elevated radio ducts so 
formed have M-inversion layers of exceptional magnitude and_ steepness, 
changes in M of 30 units sometimes occurring in 150 feet. 


5.2. Experimental results of the soundings 


5.2.1. Experimental results for two flights made on the 10th and 11th January, 
are presented in figures 8 (a), 8 (b), and figures 9 (a), 9 (b) respectively and the 
ray analysis for the first flight in figures 8 (c), 8(d). Figures 8 (a) and 9 (a) give 
the contours of constant M for the outward run, and figures 8 (b) and 9 (d) give 
(M, h) curves at each 10 miles of this path. Thus (a) and (4) present the vertical 
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Figure 8 (a). Contours of //=constant for aircraft sounding, outward flight 10 Jan. 1946, from 
Collaroy, Sidney, N.S.W. 
Lines of constant 7 ———— 360. Lines of maximum M— - — - Lines of minimum M ———-—-—- 


The values of S/N are written alongside the track of the aircraft and also indicated by the 
distance apart of the lines drawn on each side of it. Where the aircraft is not seen the track is 
shown dotted thus 7... ..... 


T=position of the radar set at Collaroy. 
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Figure 8 (4). M,h curves at a number of ranges for the stratification of figure 8 (a). 
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Figure 8 (c). Calculated ray paths from Collaroy radar for the stratification of figure 8 (a). 
Direct ray leaving transmitter at given angle of elevation (x 107 radians) —2 ; reflected ray +2. 


Lines of maximum M—-—- Lines of minimum M—----—--— S/N shown as in figure 8"(a). 
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Figure 8 (d). Values of S/N calculated from divergence of rays in figure 8 (c). 

Figures in circles are calculated values of average S/N at range of height lines beside which they 
are written ; average calculated for the height between the short horizontal lines. These values 
are represented by the horizontal lines terminated by arrow heads to the same scale as the measured 
S/N ; height lines are dotted when total S/N <1. Figures at top of height lines are twice the 
calculated free space S/N for the given range. 
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DISTANCE FROM COLLAROY RADAR ON BEARING 120°T IN MILES 


Figure 9 (a). Contours of M=constant for aircraft sounding, outward flight, 11 Jan. 1946, from 
Collaroy, Sydney, N.S.W. For details see figure 8 (a). 
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Figure 9 (b). MM, h curves ata number of ranges from the stratification of figure 9 (a). 
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section through the atmosphere, at right angles to the coast, which coincides. 
with the axial plane of the transmitter. 

The track of the aircraft is shown in figures 8 (a) and 9 (a). _ The variation 
of the signal-to-noise ratio (S/N) is indicated graphically by drawing two lines, 
one on each side of the track, at a distance apart proportional to S/N. Values 
of S/N are written alongside. Where no echo was detected the track is shown 
dotted. For comparison, the positions of the first and second theoretical lines 
of maximum field for k=4a/3 are also drawn. 

The ray diagram 8 (c) was computed from figures 8 (a) and 8 (6) and from 
this were calculated the theoretical values of S/N shown in figure 8 (d). 


5.2.2. Accuracy of the data. The values of (4—1) were calculated from 

Stickland’s (1942) formula :— 

p—-l=(QPH lel +3'S x10 1 7)\10- <0) eee (30) 
the constants being rounded off to these values. P is the total pressure of the 
air, e the partial pressure of the water vapour (pressures in millibars) and 7” kK. 
the absolute temperature of the air. 

Temperatures were read to +0-1° cc. and corrected for plane speed. Apart 
from the question of the accuracy of the formula constants, it is considered that 
the overall accuracy in Mis + 1 M-unit, so far as individual readings are concerned. 
Accuracy of the contours is affected by the time factor involved in the duration 
of a flight and by the accuracy of the determination of position. The outward 
journey took about 50 minutes. A comparison made with the return flight 
(not given here) showed that the change in structure during the time of flights _ 
was small enough to justify using the contours for field strength calculations, 
so far as the time factor is concerned. Drawing of the contour lines by joining 
points, in some cases 20 miles apart, could not allow for small scale undulations. | 
and undoubtedly the curves could not be as smooth as shown. However, the fact 
that three consecutive points generally lie near a straight line would indicate 
that large scale undulations did not occur. Curves (6) were plotted from the | 
contours of (a), and so again are somewhat smoothed. Bearing on this point, | 
which involves the scale of the temporal and spatial homogeneity of the atmos-- 
phere is the fact that the original (MW, h) curves taken along the slant track of the 
plane were, in fact, comparatively smooth. At the lower levels readings were 
usually made at height intervals of 15-20 feet and departure fromi the degree 
of smoothing represented by (b) seldom exceeded the accuracy of the M values, 
+1M-units. A point of interest in the curves (4) is the sharpness of the turning | 
points where the gradients of 7 change sign. 


5.3. Analyses of the soundings 


5.3.1. Comparison of S/N with contours of M constant (figures 8 (a), 9 (a)). 
Examination of the figures reveals the following broad features also possessed 
by the other flights of the series :— 

(a) The field strength, as indicated by the signal-to-noise ratio has a value 
which remains high from a point about the middle of the first inversion down 
to near sea-level, where the aircraft becomes lost to view. A strong maximum | 
occurs near the first maximum of M (usually a little below it) so that the position. 


of maximum field lies along a line which follows the inversion upwards as the 
latter rises out to sea. 


Radio shadow effects produced in the atmosphere by inversions — 75 


(5) The field falls off rapidly shortly after the inversion is entered from below, 
Sometimes to a value below receiver threshold, through a height interval which 
Increases with range; above this S/N increases with height. This is shown by 
the values of S/N obtained during the second ascent and the first and second 
descents of the plane. The later climbs were not high enough to show the final 
increase of field. The upper boundary of the region of diminished field is thus 
an upward sloping line whose location is shown in §5.3.2 to be near the upper 
critical rays. ‘The lower boundary is within the M inversion. 


(c) The first lobe has disappeared in figure 9 (a). In figure 8 (a) there 
appears what could be an interference maximum representing the first lobe 
displaced upwards about 150 feet from its position with k=4a/3, and another 
at 57 miles 2650 feet, this time displaced upwards 600 feet. A neighbouring 
maximum at 55 miles would appear to be due to passage across the shadow 
-boundary at this point. 


(dq) Above the upper boundary the field tends to resume its normal lobe 
structure as indicated by maxima at 30 miles 2050 feet in figure 8 (a) and 29 miles 
1650 feet in figure 9 (a), representing the second lobe. In figure 8 (a) the position 
of the second lobe, like that of the first lobe, is high, while in figure 9 (a) the 
second lobe appears in its normal position. 


5.3.2. Comparison of S/N with ray tracing. Rays were traced for the 
first sounding, covering the range of angles of elevation measured at the trans- 
mitter from 10x10-% radians to —10x10-* radians. These are plotted in 
figure 8 (c). ‘The direct rays (negative «) are shown as full lines and the reflected 
rays (positive «) as dotted lines. The rays will be referred to by their initial 
angles of elevation in units of 10-° radians. The lines of maximum and minimum 
M and the variation of .S’/N along the aircraft track are repeated from figure 8 (a). 

It will be noticed that the elevated duct of figures 8 (a), 8 (d), is double, to. 
a range of about 75 miles. The first M inversion, which falls toa deep minimum, 
is the dominant factor affecting the ray paths. This inversion is followed by 
a smaller duct whose influence is mainly noticed on the reflected rays. ‘There 
is a further complication in the M-structure between 40 and 90 miles commencing 
at 2000 feet where another M-inversion occurs. On account of its height this 
has only a minor effect. 

Both the direct and the reflected rays are strongly trapped by the first inversion. 
Splitting of the direct radiation occurs at about 11 miles 700 feet, the —6:5 ray 
penetrating all the inversions while the —6-4 ray is trapped by the first inversion 
(see figure 8 (d)). 

The critical rays are not so clearly defined for the reflected radiation. The 
first inversion separates the + 6-5 ray from the +7 ray at about 40 miles 850 feet, 
and the critical ray would be about +6-8. However, rays up to +8 become 
trapped by the second inversion, trapping being aided by the rise in the position 
of the inversion. This results in the full shadow being divided into two regions 
of diminished field, separated by a caustic. ‘The upper region extends to about 
the +9 ray, but overlaps the semishadow, the two top edges crossing at 36 miles. 

An estimate was made of the limits of the field strength values above the duct 
in terms of S/N, from a consideration of the divergence of the rays. It is felt 
that the experimental figures obtained from the radar measurements did not 
warrant closer calculations by taking into account interference effects or by 
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using the extension of ray theory mentioned in $2.1, but at the same time the | 
results correlate well enough to suggest that more detailed field strength measure- | 
ments would justify fuller analyses. The expected S/N ratio for the direct | 
and the reflected rays was estimated from the approximate relations: | 


_ So%p 5 es Sorter 31 
SS Ts Bes 
Smax=Sp+ SR} Sn =| Spo ais A DOrC (32) 


where S=average S/N ratio at range d, between two rays initially inclined 
to each other at «, 
h = height interval between the two rays at range d and initial inclination «, 
S) = value of S in free space at unit range. 


The subscripts D and R stand for the direct and the reflected rays respectively. | 
In equation (31) the earth’s divergence factor is neglected and it is assumed | 
that its reflection coefficient is unity. The value of Sy was obtained by assuming _ 
that Sp, Sp are in phase at 30 miles 2000 feet, where the measured signal-to-noise | 
ratio passes through a maximum, with Sp + Sz=10. 

The results are shown in figure 8 (d). In this figure vertical lines are drawn 
at a number of ranges and on these lines are marked their points of intersection _ 
with the rays. The average values of Sp, Sz, calculated for the height intervals | 
between the intersection points, are written in circles each in its proper interval 
alongside the height line, and actual values of Sp and Sp are represented by | 
horizontal lines terminated by arrowheads placed against their written values 
-and drawn to the same scale as the representation of the measured values of S/N. 
Where Sp+Sp<1 the height line is shown dotted. In this way the two 
‘quantities, calculated maximum S/N(=Sp+ Sp) and measured S/N may be 
visually compared. ‘The number written at the top of each height line is twice 
the free space value for S/N at the particular range. The critical rays mentioned | 
in this section are shown by shading. The field at the caustic formed within | 
the shadow region cannot be determined by this method of calculation. 

The lower boundary of the shadow, which is also the upper boundary of the | 
radio duct, has been placed some 50-100 feet above the envelope of the rays | 
trapped by the first inversion. 

Comparison of the measured and calculated fields shows that 


(a) The position of the calculated shadow boundaries agrees quite well with 
the radar results ; 


(6) The measured values of S/N above the duct are, in general, in good 


accord with the requirement that they lie between Sp+Sp, an 
| Sp — Srl. 


Discrepancies which occur at 38 miles and to a lesser extent at 48 miles 
within the semishadow are larger than the experimental limits and deserve 
examination. ‘The recorded values of the minimum of S/N occurring at each of 
these places is 1, whereas the calculated values of (Sp — Sx), although passing 
through a minimum, only fall to 4 in the first case and 2 in the second. The 
field is almost entirely due to the reflected ray. A possible explanation of this 
discrepancy is that it occurs at the commencement of the semishadow where 
it is narrow, and that edge effects which tend to make the boundaries diffuse 
may extend to a big proportion of the total thickness. Thus leakage from the 
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direct rays across both boundaries and partial reflection of the reflected rays 
at the strong M inversion would both tend to decrease the value of S/N below 
the lower calculated limit | Sp —S,|. 

The methods used cannot of themselves give a quantitative estimate of the 
complex field distribution within the duct. However, the experimental fact 
that the signal is a maximum near the maximum of M is in accord with the 
bunching of rays and the occurrence of caustics at several places near these 
maxima. A full delineation of the caustics and their cusps would require the 
drawing of considerably more rays; the intersection of rays belonging to the same 
pencil may be observed at the following places, points in figure 8 (c): direct 
rays 20 to 35 miles, about 400 feet; 60 to 85 miles, 800-1100 feet; reflected 
ray, 55 miles (intersection of ray +6 and +6:5). 

Although a quantitative estimate of the field distribution within the duct 
is not feasible by these methods a calculation can be made of the average S/N 
across the duct at a given range on the assumption that the energy remains. 
trapped without appreciable leakage, but is attenuated through horizontal 
spreading. It is of interest to note that an average value calculated in this way 
does lie within the limits observed for S/N at the same range. For instance 
assuming the angle of the trapped radiation to be 13 x 10-* radians ( —6-5 to + 6-5); 
and the vertical spread at 110 miles to be 1300 feet (say from 300 feet to 1600 feet) 
the theoretical value of S is 2-5; the measured values vary from 4 to 1:5. 


§6. CONCLUSIONS 

The results of the last section appear to confirm, for frequencies around. 
200 Mc/s., earlier theories of the ‘existence of shadow zones during certain 
predictable weather conditions. ‘The results also indicate the utility and 
practicability of ray tracing methods for obtaining the broad outlines of the 
field distribution. 

It is suggested that those using high-frequency radio systems for the detection 
or guidance of aircraft in localities subject to strong inversions should keep in 
mind the possibility, at times, of failure in the system from this cause. 
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Equivalent Path and Absorption for Oblique Incidence on a 
Curved Ionospheric Region 
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ABSTRACT. ‘The problem of oblique incidence on a curved ionosphere with a Chapman 
distribution of ionization density is discussed. It is shown that absorption, equivalent 
path, and the projection of the path on the earth’s surface may all be expressed by equivalence 


theorems of the plane earth type with suitably modified parameters and functions. 
Numerical values of the functions involved are given. 


Say INT ROD UiC MiOmN 

N a recent paper (Jaeger 1947, referred to below as V) the equivalent path and 
absorption * for a wireless wave incident vertically on a Chapman region were 
calculated. This work is extended here to cover oblique incidence on a 
‘curved Chapman region in which the ion-density is a function of height only, 
and the effects of the earth’s magnetic field are neglected. Other studies of the 
oblique incidence problem have either assumed a parabolic distribution of ion-. 
-density (Appleton and Beynon 1940, 1947; for other references on the subject see 

the bibliography in these papers) or have assumed special properties of the region 
for example, that it be thin. 
The present work consists essentially of an exact solution of the problem of a 

ray incident obliquely on a curved Chapman region. Since, even for vertical 
incidence, a final stage of numerical integration is necessary, the accurate integrals 
for oblique incidence in any particular case are not much more difficult to evaluate 
than any approximations to them. An approximation is-introduced below in 
order to simplify the presentation of the work by reducing the number of para- 
meters in terms of which the solution is expressed. At all stages the validity of 
* Absorption in a Chapman region was first studied by Appleton (1937) who discussed the case 


-of penetration of a region in which p is effectively unity, and obtained the form i 
g yy ula V (26) which I 
attributed to Best and Ratcliffe (1938)—not (1928) as in V. He also considered the a of v §8 


in which v? may be compatable with f?, and ” 1 i 
5 gave the result V (45) with a graph of 
involved but not the explicit expression V (48) for it. ee 
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the approximation has been checked by evaluating the accurate integrals and the 
discrepancy found to be negligible. 

It is found that critical frequency, absorption, equivalent path, and distance 
traversed along the earth’s surface can all be expressed by simple generalizations 
of the equivalence theorems deduced by Martyn (1935) for a plane ionosphere. 


§2,. FUNDAMENTAL FORMULAE 


We have to consider the path of a ray in an 
ionized region whose ion-density is a function 
only of h, the height of the point considered above 
the earth’s surface. Let 7 be the angle of incidence 
of the ray at height f#, and « that at zero height, 
that is, at the earth’s surface (figure 1). Then 
if R is the radius of the earth, and p the refractive 
index at height h, the fundamental equation for 
refraction in a curved stratified region is 


u(R +h) sint=constant=Rsin«. 


All the integrals which will be required below 


involve the quantity Ficute 
ds dh sect _ dh e 
po | VP LRRFAP sin® of? 


where ds is an element of the path of the ray. We proceed to study the quantity 
(2) in detail; when this has been done the evaluation of the integrals is an easy 
matter. ; 

The refractive index is given by the relation 


AN a2 ran Os See we, (Pe cbeeyas (3) 


where f is the frequency of the wave, « and m are the charge and mass of the ions, 
and N is the ion-density. 
In a Chapman region N is given by 


N=W,expa(l=—s—secye) 0 ee (4) 


where x is the zenith angle of the sun (assumed constant over the region of the 
ionosphere concerned), and 


ot ,) ie ee eens (5) 


is the height above the datum level 4, measured in terms of the scale height H 
of the region (cf. V §2, or Chapman 1931). The quantities No, H, and hy may be 
obtained from analysis of vertical incidence data by a method outlined in V. No 
is the maximum possible value of N, and occurs at height /, at the equator at noon 
at the Equinox. The critical frequency at vertical incidence, f,, was found in 
V (9) to be 
ii = (6) 
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Using (3), (4), and (5) in (2), and making the substitution 


y=(scex)i cae ee (7) 
so that 
g=In(dsecy)—2Iny, = aaeeee (8) 
(2) becomes 
ds_ _2Hdy{, _, ; R? sin? « i 
pee by exp (—9")— eae dn Seo ac 
ee (9) 
where, as in V (14), we write 
b=(foF7\/ Ze) ee (10) 
To reduce (9) to a more manageable form, we introduce new parameters f, J, d, 
and a, defined by 
B=T-+A,/R+ (AiR) laGseey), a eee (11) 
sin =(sin)/B, eer (12) 
d=bsec* d= (f2)/) 2 2e) sec ee (13) 
a=(4H)/BR)tan2 1,” ee ee ee (14) 


In terms of these, (9) becomes 


ds 2H dy sec I 
pb es ae See ee i ee we Sc (15) 
where 
A=1—dyexp (—y") + tan? {1 —[1 ~@A/ BR) lng eels eae (16) 

The expression (15) is still exact; to simplify it further, the orders of mag- 
nitude of the parameters involved must be considered. Taking h, to be of the 
order of 120 and 300 km. for E and F regions, respectively, and H to be of the order 
of 15 and 60 km. in these regions, we find that as y runs from 0° to 80°, 8 runs from 
1-017 to 1-021 for E region, and from 1-040 to 1-057 for F region. If we allow 
values of « between 0” and 90°, the maximum values of J given by (12) will be 
79-5° and 74° for E and F regions respectively; also J increases steadily to these 
values as « increases from 0° to 90°. 

If (2H/BR)Iny<1, we may use the binomial theorem in (16) and obtain 

A=1 —dyexp(—y*)—alny{1+(3H/BR)Iny+....3, 7... .- (17) 
where a has been defined in (14). 

Now in all cases (H/6R) <0-01, so we may neglect all terms of the series in (17) 
after the first if 0-03a (In y)? is negligible to our order of approximation. This is 
the case in all the calculations made here, since in the integrals which have to be 
evaluated there are additional factors (such as yexp(—y?)) which are small for 


the large and small values of y for which Iny is large and the approximation 
above might be inadequate. 


Thus for purposes of calculation we may replace (17) by 
A=1—dyexp(—y?)—alny. 


This has the advantage of expressing the whole of the phenomena in terms of the 
parameters d and a. a is zero for normal incidence or for a plane ionosphere: 
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it increases from 0 to about 0-27 for E region, and from 0 to about 0-44 for F region, 
as % increases from 0° to 90°. We shall study the values 0-2 and 0-4 of a which may 
be regarded as fairly extreme values suited to E and F regions. In §3 the para- 
meter d will be related to the critical frequency f, for oblique incidence at «, and 


the results can thus be expressed in terms of the two parameters f/f, and a. 


S32) DEES CRITICAL EREOUBNG Y 


For a ray which is transmitted through the ionosphere, ju cosz in (2), and 
thus A in (15), never falls to zero. The critical frequency f, for a ray leaving the 
earth’s surface at « is the greatest value of the frequency f for which A vanishes. 
Thus for a given value of «, and therefore of a, the critical frequency is to be 
determined from the least value d, of d for which A vanishes. Using the form (18) 
of A, the condition for this is that the equation in y 


ia liuVvevarexp(—=V)eQe sn a t (19) 


have a double root. ‘This requires that in addition to (19) its derivative shall 
vanish, that is, 


Gi Von toy jexp =v") =O) ar ae (20) 

From (19) and (20) we find 
re CHEESE nae: (21) 
d,—|yvexp(—y7)((2y2=1) Inge) ke (22) 


These are a pair of parametric equations for a and d, in terms of the parameter y. 

If a=0, d,=+/(2e), and y=1/,/2, the values for vertical incidence. For 
oblique incidence only a small range of values of y greater than 1/1/2 is needed to 
cover the required range 0 <a <0-5 of a, and thus the replacing of (16) by (18) is 
justified. For the relevant range of 8 and a the deviation of results calculated 
from (21) and (22) from those from (16) is very small. For example in the rather 
extreme case a = 0-4, the value of d, calculated from (21) and (22) is 2-578, while the 
accurate formula derived from (16) with a=0-4, and (4H/6R)=0-03628, gives 
2-5/0. 

Values of [4/(2e)/d,]', calculated from (21) and (22), are shown in figure 2. 
The critical frequency f, now follows from d, by using (13); it is 


fu=helw/ (2e)/4a]* sec I. (23) a i | 7 Tal J 
ee : 
mE 0198 | 
For a plane earth this becomes <= 
f, =f. sec x, the plane earth equivalence So ahas 
theorem; this is seen to generalize \ 
immediately to the curved earth on . 
replacing « by J and multiplying by ost AS a wD ar 
the factor s/[1/(2e)/d,] of figure 2 d 
which is nearly unity. Figure 2. 
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§4. ABSORPTION IN A DOUBLE PASSAGE OF THE LAYER AT 
FREQUENCIES GREATER THAN THE CRITICAL FREQUENCY fu 


The total absorption is 
2 feds. ibe pe aoe (24) 


taken through the layer, where, as in V, 


De 1—p? 2 Nig CR : ; a 
ne ina: DrcpHifene 34(1—z-—secye*), ...... (25) 
V=Voel™, = eee (26) 


and c is the velocity of light. It is assumed that (27f)?> v?. 
Using (25) and (26) in (24), with the value (15) and (18) of ds/u, we find 


2/xdi=Hye-\cosxcosl FG @), awe (27) 
where 
Ice) 2 perry a 
F(f If; @ =4d{ __.2 ye : 
ails O84 J, dy exp (—9")—a lng) ) 


Alternatively, writing 
D(falfs =a *V/(2e)F(falfs 2), weve (29) 
(27) may be expressed in the form 
2{xds=Hye cos xsecT (f/f O(f,if; @) 
= Hye~(d,//(2e)) COs X cos I (f,/f)? TUG a). 
For vertical incidence a=0, and F(7,/f; 0) and ®(f,/f; 0) are the functions 
F(f,/f) and ®(f,/f) defined in (23) and (25) of V and tabulated there. 


For a plane ionosphere a=0, [=«, and (27) becomes 


2fuds=Hyye cosjecos x K(f, 170) es) ma ee (31) 


10 a 
This is just cose times the | [| 


absorption for vertical incidence 


— 
at the frequency fcosa as re- : 4 | } 
aL 


quired by Martyn’s equivalence 
theorem (1935). The modifica- 
tion of this result for a curved ¢ ° PERE aE) I 
ionosphere consists of replacing « < i : | 
by J, and the function F(f,/f; 0) > 
by F(fif; 4). 

The function F(f,/f; @) is 
plotted in figure 3 for the values 


F(f&/F; a) 
aie 
L 
N 


0, 0-2, and 0-4 of a, correspond- 7 

ing to vertical incidence (or plane se* 

ionosphere) and rather large ie — Be 
values for E and F regions ® Ne ie, 0-6 0°8 1-0 


respectively, F 
Pp . ; Figure 3. The absorption function for penetration. 
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§5. ABSORPTION FOR A RAY REFLECTED FROM 
THE LOWER SIDE OF THE LAYER 


In this case we require 2 \icds over a range from ground level to the highest 


point of the path. This point is determined by the value z, of z corresponding to 
the greater root y, of 


f—dyexp( —y?)—aluy=0, > eee (32) 


The values of y, for the values 0, 0-2, and 0-4 of a are shown in the dotted curves 


of figure 6. Using (5) and (8), the height h, of the highest point of the path above 
the earth’s surface is found to be 


h,=ho+HIn($secy)—2HIny,  ...... (33) 


Proceeding as in § 4, we find for the absorption 


2{x ds— ii. cos x cos/ F, (f/f, ; a), Ae (34) 

‘where 
F 5 a)=4d{ > ep yy _ Ar 
iF lfe ) y, V —dy exp CPL yates ee (35) 


In figure 4 the function F,(f/f,; a) s 
is plotted for the values 0, 0:2, and 0-4 of 
a. Small values of f are not included, 
since for these the simple approximation 6 
(18) begins to break down. ‘There is 
no difficulty in extending the curves to 
small values of f, but there is not much 
point in doing so because we have also 
assumed hat (27f )? > v?. 

It should be remarked that the fact 2 
that the curves of figure 4 cross and are 
rather more widely separated than those 
of the other figures is caused by the 


Fy(fifa; 4) 


multiplying factor d in the definition * ne fl fe ce se 
ay) of Fi(fH3 a). Curves of [/(2e)/ Figure 4. The absorption function for 
d|F,(f/f,; a), analogous to ® of (29), do reflection. 


not cross and lie close together. 


§6. THE EQUIVALENT PATH 
This is given by 
ds 
fe 


taken over the whole trajectory of the ray. In V the procedure of evaluating the 


excess of equivalent path over distance traversed was adopted, and in this way a 
convergent integral was obtained which was easy to evaluate. A similar simplifi- 
cation can be made for oblique incidence by calculating the excess of the equivalent 
path over the distance measured to the same height along the tangent to the path at 
ground level. This distance can be calculated by simple trigonometry. d 

In figure 1 let OB be the path and OA the tangent to it at ground level. Letz 
be the angle this tangent makes with the radius vector at height /, so that 


(R+A)sind’=Rsine nae (37 
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The excess of equivalent path over distance measured along OA up to any 


height / is now ; : | 
| a dh sec’. ee (38) 4 
0 (pt ke 
Using (15) and (18), and (37) with the approximation leading to (18), this. 
becomes 


p rf : | \2 
sige | (ian se (I—aln yy" | 


For a ray transmitted through the region (which by (1) and (37) emerges with 
i=i') the excess of equivalent path over distance measured along OA for a double 


traverse of the region is 
Hsecl?G Jf: -2) >) ee nee (40) 


where 


25, e 1 1 E 
Pall a=4 | ange exp (=97)) (Laila 


Curves of P(f,/f; a) for a=0, 
0-2, 0-4 are shown in figure 5; 
that fora=O0is P(f,/f) of V, so it 
appears that the effect of oblique 
incidence is to introduce the 
factor sec J in (40) and to change 
from P(j,/7; 0) to PU. ay: 

If the wave is reflected from 
the region, the excess of equiva- 
lent path in the whole trajectory 
over twice the distance measured 
along OA to the level of reflection 
is 


Eisec LP ah p.- a), 


where 


Figure 5. Excess equivalent path for penetration. 


Pies =4] {= d ek (43) 


»\(1—a In y—dy exp (—y?))# x =a In y)# 


y 
and y, is given by (32). Values of P,(f/f,; a) and y, are shown in figure 6. 

The distance along OA to the level of reflection follows immediately from the 
trigonometry of figure 1 using the value of h, given by (33). 


§7. THE PROJECTIONSORS EEE PATH OFSTHE RA y 
ON THE EARTH’S SURFACE 


From figure 1 it follows that the projection DD’ on the earth’s surface of the 
element of path BB’ =ds of the ray is 
Rdssini R®dhsina secz 


Rik 7S agi(Rempy eens el ay ae (44) 
using (1). 
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Also, the projection on the med 
earth’s surface of the element of Fae 
OA between the heights h and es 
h+dh is 4 


R dh tan’ __ R* dh sin« sec 7’ 
R+h — (R+h)? 


Thus the horizontal projection 
on the earth’s surface of the ray up 
to height A is greater than the 
horizontal projection of the portion 
of the tangent to the ray at ground 1 
level up to height # by the amount 


a k or 4 R2 dh a! : 
sin % == SSL peo 0°2 Or4 0°6 0°8 1:0 
ol B (R+h)? Fie 
46 Figure 6 
Pea ee (46) ire x 
———— Function P,(f / fx; a) giving the excess 
: nee equivalent path for reflection. 
As in (39), this is equal WO ee Function 9, giving the level of reflection. 
2H : 
ee sec / sing 


io) 


1 1 dy 
ss ik ne In y—dy exp (—y?))! (1—alny)! IE In y]? 


Thus for a ray transmitted twice through the region, the projection of the path 
on the earth’s surface is greater than the projection of the portion of tangent to the 
ray at ground level by the amount 


(H/B2)secIsinxQ(f,/f; a; 2H/BR)) see (48) 


24 0 fe In y—dy exp (— y”))# ~ (1—alny )!J y,1—(24/BR) In y}? 


The integral Q defined in (49) differs only from P defined in (41) by the 
factor [1 —(2H/BR)Iny]-? in the integrand. Since, for transmission through 
E region, 2H/BR is of the order of 0:005 the effect of this factor is small, and it is 
found that if Q in (48) is replaced by the function P(/,/f; @) given in figure 5, the 
results obtained are less than 1%, too large. 
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For a ray reflected from the layer, the excess of the projection of its path on the 
earth’s surface over twice the projection of the distance measured along OA to 


the level of reflection is 


(H/B?\secT snaOi(fij,;¢; 20) pi). a) eee (50) 
where 
ig. 
O(z i « a) 
ri is 1 1 } dy 
- if (1—alny—dyexp(—y?))! (1—alny)!J y(1—(2A/BR) In y/? ; 


and y, is defined in (32) and plotted in figure 6. | } 
As before, Q, may be replaced approximately by the function P,(f/f,; @) of 
figure 6; values obtained in this way will be too small, the error being less than 3%. 
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ABSTRACT. A brief account is given of the physical limitations which influence the 
power output of stationary and rotating anode x-ray tubes used in crystallographic work. 


Brief mention is made of various tubes which incorporate rotating vacuum seals in the 


anode assembly. An account is then given of the “ Peristron’’}, a continuously water- 


cooled rotating anode x-ray tube which completely dispenses with the need of rotating 
vacuum seals. Incorporated in the tube is a novel type of cathode assembly which permits 
of rapid adjustment to the focus from outside the vacuum envelope. 


§1. REQUIREMENTS AND LIMITATIONS ON X-RAY TUBES 
FOR CRYSTALLOGRAPHIC PURPOSES 

N registering x-ray diffraction patterns from crystalline and amorphous 
| materials it is not uncommon to require exposure periods lasting many 

hours. Several factors play their part in determining the exposure time, 
only a few of which are under the control of the investigator. To a large extent 
the exposure time is controlled by the nature of the irradiated specimen and on 
the type of diffraction pattern required. 


t From the Greek “ peristrophe ” (spinning round like a top) and “ electron ”. 
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Broadly speaking, there are two direct ways in which the exposure times 
could be reduced by a considerable amount. ‘The first is to increase the speed 
of the recording x-ray film, and, if this could be done without introducing 
undesirable large grain-size effects in the emulsion, it would be an achievement 
of the highest importance. The second method, which is the subject of this 
communication, is to increase the intensity of radiant energy at the source. 
This is the most immediately practical method of obtaining a substantial reduction 
in exposure period. 

With the exception of the Seemann—Bohlin type focusing cameras using 
curved crystal monochromators or Soller slits, the bulk of x-ray diffraction work 
is carried out with cameras utilizing pinhole apertures or with powder cameras 
using narrow slits in their collimator systems. These pinholes or slits must be 
of sufficient size to bathe the specimen completely in radiation, and, at the same 
time, confine the beam to prevent unwanted scattering ‘effects. In general, 
the x-ray beam is admitted into the collimator system through a pinhole aperture, 
or slit, which seldom exceeds 1 mm. across, and which, in special circumstances, 
may be as small as 0-25mm. or less. This means that our first major problem 
in the design of the x-ray tube is to obtain a focal spot of high intrinsic brilliance 
and with dimensions comparable with the diameter of the pinhole aperture. 
Our second problem is to position the focal spot as close to the window of the 
x-ray tube as possible in order to cut down its distance from the irradiated 
specimen. 

Any increase of focal spot dimensions beyond the size of the pinhole aperture 
merely represents so much waste of x-ray energy. Thus a tube with a nominal 
electrical energy consumption of 1 kw. is no better for our special purpose than 
a tube rated at $kw. if only one half of the focal spot is utilized, or if, by the 
construction of the x-ray tube, the distance of the focal spot from the camera is 
unduly large. In the final analysis, the ideal solution would be to reduce the 
size of the focal spot so that it can effectively replace the first pinhole aperture 
of the camera, the lateral definition of the x-ray beam now being produced 
entirely by the specimen and the focal spot together. 

Of the electrical energy dissipated in the focal spot little more than one 
thousandth part is converted into useful x-radiation, the remainder merely 
serves to heat the anode. Thus it is generally unsatisfactory to attempt to produce 
an intense beam of radiation from a small circular focus, for the local heating 
soon melts a crater in the anode. In the majority of sealed-off and demountable 
x-ray tubes of the hot-cathode type, the focus is rectangular in shape and with 
the approximate dimensions 10x1:0mm. Such a focus, when viewed along 
its major axis at a grazing angle of 6° to the anode face, presents a foreshortened 
area of high intrinsic brilliance approximately 1-0 mm. square. 

An x-ray tube with water-cooled copper anode and rectangular focal spot 
as described above, will normally be operated in the region of 40kv. with an 
electron current of 20ma. This total load of 0-8 kw. may be sustained indefinitely 
without damage to the anode, although further considerations with regard to 
the life of the filament may necessitate a somewhat lower tube current. ‘This is 
specially true in the case of sealed-off tubes. ; 

The energy dissipated in the focal spot is thus in the neighbourhood of 
0-08kw/mm?2. Any effort to increase this loading by concentrating the focal 
spot into a still smaller area, or by increasing the electrical input, merely raises 
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the surface of the anode to an unduly high operating temperature. In a tube 
which has been carefully degassed so that there is no instability from ionization 
currents, the temperature limit would appear to be determined by the melting 
point of the anode material. In point of fact, the practical limit is set much 
lower by the recrystallization of the anode material and the fluctuating thermal 
stresses to which it is subjected. For the case of copper, de Graaf and Oosterkamp 
(1938) give an upper limit of 300° c. 

It was realized very early in the history of x-ray tubes that a substantially 
larger output of radiation would be possible if the anode were set in motion 
thereby continuously exposing a fresh, cool region to the electron bombardment 
(Breton 1897). This has been achieved by oscillating the target with a rotary 
movement of a few degrees (Hosemann 1939) and by a slow translatory oscillation 
of several centimetres (Astbury and MacArthur 1945). It has also been achieved 
by a gyratory oscillation which gives the effect of continuous rotation (Du Mond 
and Youtz 1937) and by a true continuous rotation of the anode (Coolidge 1915, 
1930, Elihu Thomson 1914, Bouwers 1929). Most ingenious of all is the method 
in which the entire x-ray tube rotates about the cathode—anode axis while the 
electron stream, deflected by an external magnet, strikes the anode at a fixed 
position in space (Breton 1897, Coolidge 1917). Few of these arrangements, 
are, however, satisfactory for x-1ay crystal analysis work. 

Miller (1927, 1929, 1931, 1939) has made a detailed mathematical analysis 
of the heat dissipation in water-cooled stationary and rotating x-ray anodes. 
Within certain limits, he has shown that for a line focus 


Wax (test) =4-25(T —T,)kA(a, b), team (1) 
; 4-04(T —T,)kav/(vbpe/k) 
Wooae (moving) = ———— = 
max (Moving) = 7 (2a/72R)\(1 +1n4R/a)\/(vbpe/ky’ (2) 
= 4-04(7.— 1 \kav/(obpc/R); 5) ae ea eee (3) 


where W,,., = limiting input in watts, 7 = maximum safe temperature of surface 
Y)=temperature of cooling metal-water interface, a,b=half length (width) 
of focal spot—assumed elliptical, k=thermal conductivity of target material 
(watts/cm/°c.), p=density of target material, c=specific heat of target material 
(watt.sec/gm/°c.), v=relative velocity of focal spot (cm/sec.), R=radius of path 
of focal spot on target, A(a,b)= Lim A, = Lim B,, where A, =4(A,_,+B,, 4) 


N>D N—> 1 
and Aj =a,.By—0. 


Thus, by giving the anode the velocity v relative to the focal spot, we can 
increase the performance of the x-ray tube in the ratio 


Wax (moving) ir 0-95 1ay/(vbpc/k) 
Ve (hese) A(a, 6){1 + (2a/7?R)(1+1n4R/a)s/(vbpc/k)y 4) 


In the case of a copper anode, pc/k™1. Fora focus of dimensions 2a = 0-8 cm, 
26=0-05 cm, A(a,b)=0-15. With an anode of radius R=9-0.cm. which rotates 
at 1370 r.p.m., we easily find that W,,,. (moving)/Wy,ax (rest) =11-3. Equation 
(4) from which the value 11-3 is derived assumes, among many things, that the 
x-ray tube is operated at a steady potential. It does not follow then, that if 
A.C. is applied across the tube the same increase in efficiency is to be expected : 


in practice, the output is raised by a factor more nearly approaching half the 
theoretical limit. 


On a new type of rotating anode X-ray tube 89 


§2. TYPES OF ROTATING ANODE TUBE 

In the case of water-cooled rotating anode tubes we are faced with the dual 
problem of rotating the anode within the vacuum envelope, and, at the same 
time, introducing a continuous flow of cooling water into the anode. A number 
of solutions to this problem have been described in the literature. 

The first water-cooled rotating anode tube, developed by Stintzing (1926); 
was designed primarily for x-ray spectroscopy and was not suitable for crystal- 
lographic work. The x-ray tubes developed by Miller and Clay (Miller 1929 b, 
Clay 1934, Muller and Clay 1939) use a cooling system closely resembling that 
of Stintzing. The anode is in the form of a flat hollow cylinder driven by a 
hollow axle which is coupled to an electric motor and rotated at 2000 r.p.m. 
The flat face of the cylinder, on the side opposite the driving axle, serves as the 
effective anode surface, the focus marking out a continuous narrow track near 
its periphery. Cooling water is injected through a co-axial inlet tube housed 
within the hollow driving shaft and is forced against the inside wall of the anode 
periphery and then returns by way of the hollow shaft to the outlet pipe. Because 
of its fundamental importance the arrangement in Clay’s 5kw. x-ray generator 
is shown diagrammatically in figure 1. 
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Figure 1. Rotating anode x-ray tube with stuffing box vacuum seal. (Clay.) 


Perhaps the most important part of the assembly is the vacuum seal surrounding 
the rotating shaft. This consists of a stuffing-box packed with hat-leathers 
pre-treated with vacuum grease and backed with high-speed vacuum pumps. 
A very similar arrangement is described by Beck (1939). 

A stuffing-box type of seal has a number of drawbacks. It leaks when the 
backing pumps are shut down, and the frictional forces, which produce an 
appreciable wear on the driving shaft, must be overcome by employing a 2 H.P. 
electric motor. 

These drawbacks led Astbury and Preston (1934) to make use of an almost 
frictionless permanent liquid seal. This takes the form of a column of mercury 
of barometric height surrounding the hollow driving shaft of the anode. A 
detailed description of the apparatus and its performance is given by MacArthur 

4). 
evo nae and compact seal is provided by a greased conical vacuum 
joint on the rotating hollow anode shaft. (Fournier, Gondet and Mathieu 
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1937). The speed of rotation is relatively small (60r.p.m.), and this limits 
the energy output. 

A description is given by Linnitzki and Gorski (1936) of two rotating anode 
tubes which eliminate the conventional rotating seals on the driving shafts, 
replacing them with molecular pumps of the combined Gaede—Hollweg or 
Gaede-Siegbahn type. No performance data for these tubes are given. 

From the brief descriptions of the rotating anode x-ray tubes given above, 
one feature common to them all becomes evident. That is, the cooling water 
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Figure 2. Section through x-ray tube. 


Is led to the inner surface of the anode through a hollow driving shaft, which 
in its turn, must pass the wall of the x-ray tube by way of some form of rotating 
vacuum joint. The Peristron (Taylor 1944) described below, is a rotating 
anode x-ray tube using continuous water cooling, and eliminating entirely the 
need for rotating vacuum seals. 


§3. DESCRIPTION OF THE PERISTRON 


Figure 2 is a vertical section through the rotating anode x-ray tube. Motive 
Po to the Totating anode is given by a 3-phase 400-volt induction motor 
aving a specially large air-gap 5/32’ across and therefore consuming only a 
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small fraction of 1u.p. The speed of rotation is 1370r.p.m. The brass base 
of the vacuum envelope of the x-ray tube is accurately spigotted into the cylin- 
drical top of the induction motor which serves as a very rigid support. Brazed 
into the base is a right-angled vacuum pipe 14’’ in diameter. The centre of the 
base houses the ball-races for the driving shaft in a central tubular member, 
the lower end of which is terminated by a waxed-on }” thick glass dome. 
This oe the air-gap of the motor and effectively seals the rotor, at the lower 
end of the driving shaft, within the vacuum env elope. 

Screwed into the upper end of the shaft is the detachable hoilow anode. 
The whole sub-assembly, comprising anode, base-plate, rotor and glass dome, 
removed from the mctor housing, is shown in figure 3. 

The water-cooled brats top of the x-ray tube is vacuum sealed to the base- 
plate by means of a flat neoprene gasket. The flanges gripping the gasket in 
position are turned with a series of fine concentric ridges. These bite into the 

neoprene so that the atmospheric pressure bearing down on the top is quite 
sufficient to maintain a perfect vacuum seal which is both easily demountable 
and trouble-free. ‘The cathode assembly is mounted on the top by means cf a 
demountable plasticine vacuum joint. ‘The window, which is of lithium metal 
1mm. thick, is pressed into position by a screw cap and is easily replaced. 


(1) The anode assembly 


The anode with its associated cooling system comprise the special feature 
of the x-ray tube. It is essentially a hollow cylinder, the base of which is of 
brass 83’ in diameter and 4” thick, with a screwed boss on the lower surface 
for attachment to the driving shaft of the electric motor. The upper anode 
surface is made of pure copper sheet, }” thick and with a slope of 6°, the 
focal spot delineating a track approximately 10mm. wide at a distance of 1” 
from the periphery. In the original construction, the bush at the centre of 
the anode surface and the cylindrical walls were of brass and were brazed to the 
copper anode, but a much simpler construction was evolved in which the bush, 
anode surface and cylindrical walls were spun from a single copper sheet. 

A flat double-spiral of thin-walled copper tubing, 5/16’’ external diameter 
and 1” bore is held stationary within the hollow anode by the spindle which 
passes through the centre bush in the upper surface of the rotating anode. Water 
is introduced through pipes passing through stationary demountable seals in 
the base of the vacuum envelope and passes, in turn, through the spindle into the 
spirals within the anode and out again through the spindle. 

The water-cooling leads to the spirals are bridged together by a brass plate 
which bears on a rubber gasket, thus forming a demountable vacuum seal covering 
the holes where the nipples on the pipes pass through the base. The bridge 
also serves to prevent any twisting of the pipes as the Jocknuts on the vacuum 
seals and pipe nipples are tightened. 

The anode is partly filled with vacuum oil through the large inlet hole drilled 
in the spindle which is then closed by means of a screwed filler-cap. About 
i litre of oil is sufficient. The surface of the oil and its vapour are in active 
munication with the high vacuum within the x-ray tube through the gap 
between the spindle and the “bushing. Thus it is necessary to use an oil w ith an 
- extremely low vapour pressure such as Apiezon ““B”’, Octoil “5S” ora silicone. 
In the present instance Apiezon “B” with its vapour pressure of 10°°mm. 
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has been successfully used. Any possible creep of oil through the gap between 
the spindle and bush is prevented by a screw thread cut into the bush which 
drives the oil back into the anode. A clearance of about 10/1000’ is required 
to prevent the bush seizing on the spindle, since the vacuum oil has very poo: 
lubricating properties. 

In the course of rotation, the centrifugal forces fling the oil against the circum- 
ferential inner surface of the anode. By this means, the heat generated in the 
metal surface of the anode by the electron bombardment is instantly transferred 
to the oil, which in its turn is transferred to the cooling water flowing through 
the stationary spirals. Moreover, the churning motion imparted to the oil 
as it is forced past the spirals provides an extremely efficient method of transferring 
the heat, for the cooling is not dependent on natural convection currents as in 
a static body of liquid. With a flow of water of only 10 litres per minute, and an 
initial temperature of 14°c., the outflow temperature was 20°c. when the power 
dissipated in the anode was 5 kw. 

The design of the x-ray tube is such that the anode is easily accessible for 
cleaning and for replacement whenever different characteristic radiation 1s 
required. Upon lifting away the lid from the rubber gasket which forms the 
vacuum seal with the base-plate, the whole anode assembly is exposed (figure 4). 
To remove the anode, the water-cooling pipe nipples and locknuts on the base- 
plate are unscrewed and then the anode is unscrewed from the driving shaft, 
which is held stationary by means of a flat spanner. The driving shaft, bearings 
and rotor are thus left behind as an integral part of the base-plate and need never. 
be disturbed. ‘The time required to strip down the tube and replace the anode 
is little more than 5 minutes. 

After re-assembly it is important to rotate the anode before commencing 
pumping operations. This prevents any tendency for oil to creep past the 
bush should it foam during the escape of adscrbed air. For this reason a thermal 
delay switch was incorporated in the electrical controls which gave a few seconds 
for the anode to rotate up to speed before allowing the rotary pump to start. 
Once vacuum conditions were established the anode could be rotated or stopped 
at will independently of the pumps. 

The x-ray tube is mounted on a steel stand which houses the vacuum pumps 
and also serves to carry an electrical control panel. "The pumping unit comprises 
Metropolitan Vickers O2 and O3 oil condensation pumps mounted in cascade 
and backed by a DR1 type rotary pump. , The cooling water’from the pumps 
is arranged to pass into the cooling spirals of the anode. This is a more efficient 
arrangement than passing the water in the reverse direction since the temperature 
rise on passing through the pumps is negligibly small. It has the advantage 
over duplicate sets of cooling leads that only one water relay safety switch is 
required in the electrical circuit. 

Directly beneath the platform supporting the x-ray tube are two small air 
compressors. One circulates a stream of cooling air round the stator windings 
and glass dome of the x-ray tube, the other serves to cool the cathode assembly. 


(11) The cathode assembly 


In the majority of demountable x-ray tubes, the provision made for adjusting 
the position of filament and focusing tube with respect to the anode is very crude, 
and any alterations must be made by opening up the vacuum system. ‘This is 
a time-consuming process with the result that an operator is too often content 
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Figure 3. Anode sub-assembly removed from motor-housing. 


Figure 4. Rotating anode x-ray tube with top removed, exposing anode assembly, 
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Figure 6. Transmission 
with stationary ano 


(a) and back-reflection (d) 


de x-ray tube, CuKa+ 8 radiation. 
20 minutes. (c) & ( 
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photographs of 0:0015” copper foil, taken 


40 kvp. 20 ma, Exposure time 
d) Rotating anode tube. 40 kvp. 100 ma. Exposure time 3 minutes, 


On a new type of rotating anode X-ray tube 93 


to let the tube run with an inefficient broad focus. Moreover, it is not generally 
realized that if the tube is focused for running with a chromium anode at 35 ky. 
a change to 90kv. for operating with a silver or tungsten anode will have a marked 
defocusing effect. It was in order to overcome these difficulties that a new ad- 
justable type of cathode assembly was developed. 

In the present arrangement shown in figure 5, the position of the cathode 
assembly relative to the anode is altered by means of a flexible metal bellows 
moved by a nut, the cathode adjuster, situated outside the vacuum system on 
the metal cap sealed to the glass insulator. A second internal bellows is connected 
to the filament holder, which may be moved up or down within the focusing 
tube by means of the filament adjuster nut, which is soldered to a porcelain 
insulator above the cathode adjuster. It is therefore possible to make adjustments 
and obtain the best focusing conditions in a few seconds by observing the changes 
in a pinhole image of the focal spot upon a fluorescent screen—an operation 


which normally takes several hours with the more orthodox type of cathode 
arrangement. 


Figure 5. Cathode assembly. 
List oF PARTS 


-1. Filament assembly 9. Glass insulator 16. Metal-porcelain vacuum 
2. Focusing tube 10. Screwed tube joint 

3. Steel tube 11. Clamping ring 17. Bellows-cap vacuum 

4. Filament lead 12. Cathode assembly joint 

5. Porcelain insulator adjusting nut 18. Flexible filament lead 

6. Flexible bellows 13. Filament retractor nut 19. Porcelain guide 

7. Flexible bellows 14. Soldered vacuum joint 20. Mica insulation 

8. Brass cap 15. Air pipe 21. Wax vacuum seal 


Experiments were made with a number of filament sizes and types of focusing 
tubes. Ultimately a simple combination was found which gave a sharp focus 
of 10mm. x 0-5mm. at 45kv., and 100ma. ‘This consisted of a tungsten spiral 
of 16 turns of 0:33 mm. diameter wire wound on a mandril 0-8 mm. diameter 
into a helix 10mm. long, housed within a steel focusing tube with a simple slot 
24mm. across and 12mm. long. The slot opens out from the surface nearest 
the anode towards the filament at an angle of 45°. ‘The focusing tube stands 
15mm. from the target with the filament surface approximately 1mm. behind 
the front surface of the slot. This latter distance is most critical, for a movement 
of mm. in either direction will produce considerable changes in the sharpness 
of focus. . With the target bevelled at 6° the effective focal spot is a rectangle 
1mm. high by 0:5mm. across. A filament current of 8-5 amp. is sufficient for 
atube current of 100ma. Excellent results have also been obtained with tantalum 
filaments. 

The filament is held in position in the simplest possible manner. ‘The 
straight ends of the filament spiral are pushed into 1mm. diameter holes drilled 
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in the filament holder and held rigidly in position by driving in thin soft iron 
wedges, These have the great advantage of being easily removed when the 
filament is due for replacement. Since an appreciable amount of heat is conducted 
from the filament along the steel cathode assembly to the soldered joints on the 
bellows and to the Apiezon ‘‘W” wax seal between the glass insulator and the 
end cap, some form of cooling is desirable. A stream of cooling air, generated 
by a small compressor mounted below the x-ray tube, is circulated through the 
bellows system via an inlet pipe soldered into the cathode adjuster nut. Such 
cooling only becomes necessary if the tube is to be run continuously for periods 
of over an hour. 


§4. PERFORMANCE OF THE X-RAY TUBE 

The x-ray tube may be operated continuously at 100-120ma. and 45kv. 
without giving any trouble. It has also been run for varying periods at 30kyv. 
at 175 to 200ma. It must be stated that these figures are nominal in that raw 
A.C. is applied from the H.T. transformer directly to the tube and the resultant 
D.C. component in the transformer may cause serious distortion of the wave form. 
In this work the focal spot was 0-5 x 10-0 mm. 

As a basis of comparison, photographs were taken with the same specimen 
in the same camera with the Peristron and with a new sealed-off Victor x-ray 
tube run at its maximum output. Figure 6 illustrates transmission and back 
reflection photographs of 14/1000” copper foil taken with the two tubes using 
unfiltered Cua+-8 radiation. ‘The reduction in exposure time is from 20 minutes 
to 3 minutes. On account of its extremely sharp focal spot, the rotating anode 
tube gives slightly sharper photographs. 
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DETERS 10. THES LEOR 


The Dipole Moment of the Supposed Fundamental 
Magnetic Field due to Rotation 


In his recent important discussion of the magnetic fields of the earth, sun and stars, 
Prof. Blackett (1947) suggested that these fields represent a new and fundamental property 
of rotating matter. His formulation of his hypothesis (reviving and revising one already 
proposed by H. A. Wilson) is at present tentative, and applicable only to bodies symmetrical 
about their axis of rotation; it is that pu, the mass flux due to the rotation, is associated with 
a fundamental magnetic field as if it were an electric current intensity Kpu ; here p denotes 
the mass density, wv the rotating velocity wp (w the angular velocity and p the distance from 
the axis), and K, when the current intensity is expressed in electromagnetic units, is given 
as BG?/2c (G the constant of gravitation, « the speed of light, and 8 a factor of order unity). 

As Blackett indicated, the external fundamental field of a spherical body in uniform 
rotation, when the density distribution is spherically symmetrical, is of dipole character, and 
the following relation holds between the dipole moment P and the angular momentum U: 
P/U=3kK. Thus P/U is independent of the nature of the function p(r) expressing the 
spherically symmetrical density distribution. 

In this note I show, very simply, that P/U=$K for a rotating body when there is 
symmetry about the axis of rotation, whatever the form of the body, or the non-uniformity 
either of por o. 

Let the position of any point of a half meridian section (bounded by the axis) be specified 
by p, the distance from the axis, and z, the coordinate measured parallel to the axis. Let 
dp dz be the area of an element of the section, with sides parallel and perpendicular to the 
axis, at the point p,z. The equivalent electric current associated with this element is 
Kpwp dp dz; it flows round a circle of radius p. By Ampere’s hypothesis its field is the 
same as that of the magnetic induction B for a uniform magnetic shell of strength 
j=Kpwp dp dz bounded by the circle; it is convenient to take this shell as plane; its 
magnetic moment is j times its area (7p) and is parallel to the axis. | As moment vectors are 
additive, and as these moments are all in the same direction along the axis (except in so far 
as w may vary in sign), the total dipole moment P is given by 


P23p4= | | nK pup? dp dz. 


In this expression p and w may be any functions of p and z. 
The mass of the ring associated with the element dpdz is 2mppdpdz, and its speed 
is wp; its angular momentum is therefore 27pwp* dp dz, and thus by integration 


ei 2n| pop*dp dz. Hence P/U=3K, as stated. 


If C is the centre of the current ring associated with the element dp dz (so that C is the 
axial point with coordinate z), the potential of the field of the ring, at a distant point, can 
be expressed as a series of zonal harmonic terms, of which the first is determined by the 
dipole moment of the ring. If the potential is expressed relative to any other point O on 
the axis, this first term is unaltered, but the coefficients of all the higher harmonics are 
changed. When the contributions of all the rings are combined to give the whole potential, 
the first harmonic corresponds to the dipole moment }KU, and in addition there are in 
general higher harmonics which do not vanish whatever the choice of O. These do 
vanish, however, for a spherical body rotating about a diameter, when p and w are functions 
only of the radial distance. rh 

In a paper on “‘ Solar magnetism and the hypothesis of fundamental magnetization by 
rotation”, to be communicated to the Royal Astronomical Society, some examples of these 
higher harmonics are given for a non-uniformly rotating sun. _ 


Queen’s College, Oxford. , S. CHAPMAN. 
26th April 1948. 


Brackett, P. M. S., 1947, Nature, Lond., 159, 658. 
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The Freezing-in of Nuclear Equilibrium 


In a previous paper (Ubbelohde 1947) a method was described for testing whether the 
abundance ratios of stable isotopes might be due to the freezing-in of thermodynamic equili- 
brium. The test was constructed so as to avoid uncertainties about the density of the system, 
and about the mass of the neutron, at the time of freezing-in. ~ From the freezing-in curve 
a freezing-in temperature was calculated of approximately 101° °K. 

Further data have since become available (Mattauch 1946, Stephens 1947). Where the 
masses of the stable isotopes are known to at least four places of decimals, it is possible to test 
whether the additional abundance ratios are in agreement with the freezing-in curve. 
Uncertainties in AA greater than 0:0002 prevent any decision for all the elements with 
Z~>28. No conclusions can yet be reached about their non-thermodynamic origin, on 
account of this gap in the data. 

New values of AA of sufficient accuracy for the test include the following : 

Abundance Abundance 


m+14, , ent Meee m+2 4, by OI os 
m-+1)/m - m--2)/m 

2>He 107 0:986872 AKO) 0:0020 2:00485 
1183 4-00 0:996732 6ONi 0-396 1-99010 

13 0-011 1:003681 82Ni 0-142 1:99978 
EOIN 0-00384 0-997340 84Ni O22 1:99785 
1@ 0:0004 1-00450 
SINi 0:0449 1:00419 


When these values are plotted on the curve of —log ry against AA, previously established 
(Ubbelohde 1947) (see figure), it will be seen that for the isotopes 'B, 1°C, 170, 18O, ®°Ni, 
52Ni, ®4Ni there is good agreement within the limits of uncertainty previously found. 
(To avoid confusion, not all the earlier data are included in the figure.) For the isotopes 
°He, *N and ®'Ni no thermodynamic relation can be found, unless the reactions main- 
taining equilibrium are very notably slower and faster respectively than the average. 


-Log r 


099500 1-00000 100500 1-Olooo 
AA 

Activity of the neutron at freezing-in 

The value of the intercept made by the freezing-in curve when log r=0 is AA=0-99780 
This intercept may be used in conjunction with the freezing-in temperature and an eeiuecd 
value for the mass of the neutron, to calculate the neutron activity at the time of freezing-in 

For convenience, the equilibrium equations of the previous paper are rewritten : 
—log n+log r= —4-7 x 10!#(AA —my)/ T+ 3/4606 — 3 log T—log [(27myk)tk/h3], 2.2... (1) 
where the effects of the difference in mass between isotope pairs, and any differences i i 
multiplicities have been neglected in the last term. tcl 
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Writing log r=0, T=10", AA=0-99780, myn=1-008941 (cf. Stephens 1947), we find 
—log n= —13-87. There seems to be no likelihood of a value of 7 as large as this. The 
numerical result is, however, very sensitive to the exact numerical value taken for the freezing- 
in temperature. It is rather less sensitive to the value assumed for the mass of the neutron. 
For example, if a freezing-in temperature of 10° were used, —log n= +34:9. 

This value of m is very much too small to be probable. Clearly a temperature inter- 
mediate between 10" °x. and 10® °K. would give quite plausible values of n. But the test 
for the thermodynamic equilibrium is best constructed so as to avoid the use of n altogether, 
as indicated in the previous paper. 

It is important to note that the very rapid change in neutron concentration over the 
temperature range between 10" °K. and 10° °K. is imposed by the physical magnitudes -in 
equation (1). This rapid change provides a simple explanation for the suggestion made in the 
previous paper (Ubbelohde 1947) that freezing-in occurred on account not of a deficiency 
in activation energy, but of a deficiency in neutron concentration at around the freezing-in 
temperature. 

With regard to the effect of the mass of the neutron, using the values given above 
AA—m,=—0-01114. The estimated uncertainty in My, is +-0-00002, which would not affect 
the calculations appreciably. The maximum uncertainty in the AA intercept is about 
+-0:00150, and this would not be enough to make a serious difference in n. But at the 
temperature of 10!° °K., the R.M.s. velocity of particles of unit mass (6O=16) is approxi- 
mately 1-5 x 10° cm/sec., i.e. 0-058. Neutrons hotter than the average will begin to-exhibit 
relativity effects, which may to some extent modify equation (1). 

Thanks are due to Mr. G. E. Roberts and Mr. J. K. Brown for drawing my attention to 
the additional data quoted above. 

Queen’s University, Belfast. A. R. UBBELOHDE. 

25th March 1948. 


Marraucu, J., 1946, Nuclear Physics Tables (Interscience Publ. Inc.). 
STEPHENS, W. E., 1947, Rev. Mod. Phys., 19, 19. 
UBBELOHDE, A. R., 1947, Proc. Phys. Soc., 59, 139. 


The Detection of Heating in Liquid Dielectrics 
at Low Frequencies 


The heating in dielectrics which occurs at radio frequencies is well known and is so 
marked that its detection presents no experimental difficulty. It has indeed found many 
applications in industry since the uniform production of heat throughout the substance 
is difficult, if not impossible, to produce by other methods. The heating is very dependent 
on the frequency of the alternating electric field, and at low frequencies, say 104 c/sec., becomes 
so small that its detection and accurate measurement are extremely difficult. 

During the course of some experiments by Prof. E, N. da C. Andrade and myself 
on the effect of alternating electric fields on the viscosity of certain liquids, the dielectric 
heating of certain polar liquids not only showed itself clearly but proved to be capable 
of rough quantitative estimation at frequencies as low as 10* c/sec. In these experiments 
a viscometer is used in which the liquid flows through a narrow rectangular channel between 
massive steel plates (Andrade and Dodd 1946). The viscometer is contained.in an air 
bath immersed in a large water bath whose temperature is controlled so that its variation 
is less than 1/1000 degree over a period of many hours. The time of flow of a certain volume 
of the liquid through the channel is measured by the successive interruptions of two beams 
of light by the liquid falling in one limb of the viscometer. The amplified electric impulses 
from the photoelectric cells upon which the light beams are directed are applied to thigh 
speed relays used to start and stop a synchronous electric clock reading to 0-001 second, 
Times of flow taken at four-hourly intervals do not differ by more than 0-01, second: for 
a total time of flow of over 1000 seconds, showing that evaporation is negligible. sts) 

Alternating electric fields of amplitude 50 kv/cm. and frequencies up to 20,000 c/sec. 
have been applied to the liquid whilst it was flowing through the gap, the times of flow. of 
the liquid being measured not only while the field is applied but also both before and after 
the application of the field. 


PROC, PHYS, SOC. LXI, T 7 
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With monochlorbenzene, times of flow taken immediately after one another are plotted 
in the figure, the field of 1016-0 
amplitude 35 kv/cm., and 
frequency 10‘ c/sec., being 
applied after six measure- 6 
ments without field. The 
steady decrease of the first 
six values has been shown 
to be due to the slight heat- 
ing of the liquid caused by 
the dissipation of the work 
done in forcing it through 
the capillary. A much bigger 


Time of Flow (sec.) 
Ss 
S 


decrease occurs immedi- 4 
ately after the application of p 
the field, indicating that a ae 
new form of heating has | ee ae een eyes ee) ae ier gey (eh [Ry Te 


: : : 
occurred. This heat is Number of Run 


slowly absorbed by the thermostat, and after four or five runs the times fall accurately 
on the curye through the points obtained before the application of the field. ‘This can 
only be dielectric heating as the liquid is an almost perfect insulator, so that the ordinary 
Joule heating is absent. Confirmation of this arises from the fact that no such decrease 
in time of flow was observed immediately after the application of a steady field of 40 kv/cm. 
The increase in the time of flow which occurs whilst the field is being applied is, of course, 
due to the increase in the viscosity of the liquid and will be fully discussed elsewhere. 

The dielectric heating which occurs whilst the field is being applied can be obtained 
by extrapolation and is seen to cause 0-46 sec. decrease in the time of flow of 1015 sec. 
With the liquid in question this represents a mean temperature rise of about 0-04°. 

The purpose of this note is to draw attention to this method of detecting very small 
dielectric heating in liquids. The thermal capacity and surface area of the electrodes, 
although appropriate for the main purpose for which they were designed, are too large 
for the exact estimation of the dielectric heating. However, a rough estimation can be 
made of the power factor at these low frequencies on the following considerations. 

The heat generated in the liquid will be communicated to the steel plates forming 
the viscometer. The velocity of propagation of a simple harmonic temperature wave 
of period T through the metal of diffusivity h is \/(47h/T). Taking h for steel as 
0:15 cm? sec™} and T as 2000 sec., the value of this velocity is about 0:03 cm. sec~! and 
will, of course, be greater for the higher frequency terms in the Fourier series of simple 
harmonic terms which represents the actual applied temperature wave. ‘This means 
that the temperature wave will have travelled the thickness of the steel plates in a time 
of the order of 100 sec. ‘The damping factor 1/(z/hT) is small in this case, being approxi- 
mately 0:1, so that inside the steel at 1 cm. from the liquid/steel boundary the temperature 
excess is 9/10 of the value at the interface. Hence, during the time of flow with the electric 
field applied, the whole of the metal, of thermal capacity 260 cal. deg-!c., may be con- 
sidered to a first approximation as being at the temperature of the liquid. The heat which 
must be given to the liquid and the surrounding metal in order to produce the rise of tem- 
perature noted (0-04° c.) is thus 10:8 calories. 

During the flow the whole surface of the apparatus is losing heat to the air bath, whose 
temperature is 004° c. lower. Taking the surface area as 320 cm? and its emissivity 
as 2:5 X 10-4 cal. cm~* sec“! deg! c., the heat given to the enclosure is about 3-2 calories, 
so that the total quantity of heat generated in the liquid is 14-0 calories in 1015 sec. Since 
the volume of liquid present at any instant in the capillary is 0-143 cm, the mean rate of 
production of heat in unit volume of liquid is 9-6 x 10-2 cal. cm~* sec—1. 

Theoretically the condenser with dielectric may be considered as made up of a pure 
capacity in series with a small pure resistance. The rate of dissipation of energy W, 
per unit volume of dielectric, when an alternating field of amplitude X and frequency v 
is applied to the system is then 1X?K»v cos ¢, where K is the dielectric constant and cos d 
is the power factor. In the case considered, X=35 kv. cm—=1-:17 x10? £.s.U. cm— 


v=10* sec? and K=6, so that W=4-9 cos ¢ cal. cm-? sec—, 
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Equating this expression for W to that obtained by experiment gives the value 2 x 10-? 
for the power factor for chlorbenzene under the conditions holding in the experiment. 

Similar effects have been observed with higher fields and frequencies for monochlor- 
benzene, chloroform and amyl acetate, all of which are polar. 

No change in time of flow greater than 1 part in 100,000 was observed with benzene. 
This liquid is non-polar, and any dielectric heating can only be obtained through the 
distortion polarization of the molecule, and would, therefore, be expected to be very small, 
since the power factor is, in this case, of the order of 2 10-4. 


Physics Dept., University College, London. C. Dopp. 
20th April 1948. 


ANDRADE, E. N. pa C., and Dopp, C., 1946, Proc. Roy. Soc. A, 187, 299. 


The B-N Band-system of Silver Iodide 


In the course of work on the absorption spectrum of silver-iodide vapour in the far 
ultra-violet region (Barrow and Mulcahy 1948) it was noticed that within a rather narrow 
range of pressure the centre of the B—N system at 3200 a. was well developed. The long 
wavelength end of this system has been the subject of careful investigation by Brice (1931), 
but his measurements end with the 0,4 band at 3295 a., and in the region 3170-3295 a. 
only the rather rough data of Franck and Kuhn (1927) are available. We have therefore 
thought it useful to record our measurements, especially since they lead to somewhat 
different values of wg and D’ from those given by Brice. 

The bands were photographed in a first order of a 2:-4-m. grating instrument on Ilford 
Process plates, the linear dispersion on the plate being about 7-4a/mm. A tungsten 
filament lamp with “‘ Cornex’”’ envelope, kindly supplied by the Research Laboratories 
of the G.E.C., was used as source of continuum. ‘The optimum temperature of silver 
iodide was about 900° c. Of the sixteen bands measured, only for the 2,4 and 0,3 bands 
was the silver isotope separation clearly resolved : for the other bands it was assumed that 
settings on the low-) side of the origin were on 1°’AgI heads and on the high-A side on 
109A oT heads. 


Table 1. Band-head data for 1°’AgI 


6 | 31657 (0) 
?3157-9 
(23) 
Dies 1633-671) Italic figures denote wavelengths in air; roman, wave 
3160-28 numbers in vacuo; figures in parentheses are visual 
64:9 relative intensities. 
4 | 31568-7 (2) 
3166-78 
86-4 
3 | 31482-3 (3) 30868-4 (1) 
3175-47 3238-63 
98-2 99-6 
2 17). 206:95 311 77-2, 30768°8 (2) 200-8 30568-0 (2) 
Serer ees 3249-11 3270-45 
111-3 110-1 
1 | 31272-8 (6) 205-7 31067-1 (7) 205-8 30861°3 (3) 
4h 3196-75 3217-91 8239-37 
o, 118-4 119-0 116-7 | 
311544 (5) 206°3 30948-1 (9) 203°5 30744-6 (10) 204-2 30540-4 (8) 
: zee yn 3251-67 3273-41 
ae ST er cae 3 4 


vu -> 
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‘A correction was applied to the latter, so that all the figures in table 1, which summarizes 
the band-head data, relate to '°’AgI. 

The results, particularly for the v’’=0 progression, provide new information about the 
vibrational levels of the excited state. This is given in table 2 in the form of energies of 
levels v’ in excess of v’=0 (for !°7AgI), calculated both from the bands in table 1 and from 
Brice’s figures for !°7AgI and 1%AgI. In the calculations, Brice’s expression for Gy has 
been used. The agreement is very satisfactory, and shows that the expression for Gy 
deduced for bands with v” ranging from 4 to 17 holds well over the range 0 to 17. From 
the mean results of table 2, values of AG}, 41 were determined and plotted as in the figure. 


Course of the vibrational! intervals in che B state of 1°’AgI. 


We obtain w(=123-5, De=570, Do=510, which differ significantly from the figures 
given by Brice, namely, #=131-3, De =769 cm~. 


Table 2. Energies of levels v’ in excess of v’’=( (cm) 


v’ Present work Brice Mean AG’ 
6 B16 5:2) anal) — (31657) 

5 31633°6 (1) — 31633°-6 (23) 

4 568-7 (1) 31568-6 (4)* 568°6 65-0 
3 482-0 (2) 482-9 (14) 482°8 85:8 
2 383-3 (5) 382-3 (12) 382-6 100-2 
1 272-1 (3) 272-0 (8) 272-0 110-6 
0 153°9 (5) 15339) 153-9 118-1 


* Number of band-heads from which the values were derived, 


The unusual character of the B state (Mulliken 1937) is emphasized by these results. 
Thus, from the distribution of intensity among the bands, it appears that 7, is somewhat, 
but not very much, greater than 7’; yet the force-constant in state B is only 36% of its 
value in the ground state, while the vibrational levels converge so rapidly that Dé is only 
3% of the dissociation energy of the ground state. 

Physical Chemistry Laboratory, 
C ». Oxford. 
~°" 22nd April 1948. 


R. F. Barrow. 
M. F. R. Mutcany. 


Barrow, R. F., and Mutcany, M. F. R., 1948, Nature, Lond., in the press. 
Brice, B. A.,1931, Phys. Rev., 38, 658. 

FRANCK, J., and Kun, H., 1927, Z. Phys., 43, 164. 

Mu tuixen, R. S., 1937, Phys. Rev., 51, 310. 


Io! 


DISCUSSION 


on paper by R. Wet entitled ‘‘ The Variation with Temperature of Metallic 
Reflectivity’ (Proc. Phys. Soc., 1948, 60, 8). 


Dr. D. J. Price. The equation preceding (15) has the obvious factor of 2x—a leading 
to a root having a value $a. This root leads to A; being of the magnitude 0-1 » for most 
metals. _ If, as stated, the cubic had only one real root, the above would be the only solution. 
This, however, is also false. It is probable that the complicated substitutions and expansions 
have given stationary points for that portion of the reflectivity curve giving negative values 
for the reflectivity. (These are, of course, without physical significance.) 

From purely physical grounds, the conclusion of (16) may be seen to be unsound, since 
A, is the limit at which the Hagen-Rubens relationship becomes valid to a high approxi- 
mation, and, therefore, the anomalous behaviour associated with the appearance of an X-point 
cannot occur here. 

Taking all this into account, it would appear then that the only region in which an 
X-point may be found is in wavelengths small enough for the contributions of the bound 
electrons to be considerable. Recent work, to be published shortly, shows that this is indeed 
the case, and that a metal containing free electrons only can have no X-point. 


AuTHoR’s reply. JI am very grateful to Dr. Price for the comments he has made in his 
discussion. It is true that one solution of my equation (15) is x=a/2. As Dr. Price himself 
points out, this gives rise to an X-point at A=0-1 4; this wavelength is so far in the 
ultra-violet region that my original assumption of free electrons being responsible for the 
mechanism of reflection cannot be maintained. ‘Therefore this solution can be rejected. 
On dividing (15) by 2x—a there remains the quadratic 2x?—(3a—2y)x-+-ay=0, which 
gives 


10 2y? 
X=4}3a—2y+ ViGat—20ay-+4y5)}4| 302943 (1 oe +- at )| Bre (49) 


The positive sign gives an X-point which falls still into the region where the effect of bound 
electrons predominates, and only the negative sign yields a value comparable with those 
obtained experimentally by Dr. Price. Since y/a™1/700 it follows from (1) above that 
x=4(—2y+10y/3 ... .)=y/3, which yields my equation (16). Since (16) has been tested for 
nickel by computing values for A, from katoptric data, and the x predicted by this equation 
is confirmed experimentally, I fail to see why the conclusions drawn from this formula should 
be unsound. Nor can I understand how Dr. Price arrives at the conclusion that stationary 
points in dR/dv presuppose negative values for the reflectivity. However, experiments are 
in progress which are designed to test the temperature variation of A,; results will be reported 
on completion of this work. 


DIsSeUSsloN 


on paper by G. F.J. Gartick and A. F. Grsson entitled “The Electron 
Trap Mechanism of Luminescence in Sulphide and Silicate Phosphors ”’ 
(Proc. Phys. Soc., 1948, 60, 574). 


Dr. H. A. Kuasens. Lord, Rees and Wise (Proc. Phys. Soc., 1947, 59, 473), have 
explained their decay results also by assuming the electrons to be trapped near excited 
centres, but they allow for the possibility of such electrons escaping to the conduction band 
and combining with other excited centres. Does Dr. Garlick allow for this process ? 
Also, how does he explain that Cu introduced into zine sulphide activated with Mn prolongs 
the afterglow of the Mn emission ? 
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Autuors’ reply. The assumptions made by Lord, Rees and Wise in their paper are not 
necessary to explain our experimental results. They are based on very restricted phos- 
phorescence decay measurements on one specimen of zinc sulphide activated by zinc and 
silver with a trace of copper. No independent study of the trapping states was made using 
thermoluminescence measurements. It would be difficult to distinguish the effects | of 
certain traps being related to different emission centres in such a phosphor since the emission 
spectra of the two different centres formed by zinc and silver are almost identical. 

With respect to the effect of copper in prolonging the afterglow of zinc sulphide phos- 
phors activated by manganese, it must be assumed that electrons from traps associated with 
the copper impurity cause emission in manganese centres when they escape. However, it iS 
not essential that the transit between centres should take place via the conduction levels of 
the matrix. Such a process, involving entry of electrons into the conduction levels, would be 
temperature dependent, whereas in some specimens of this type of phosphor we have found 
no marked temperature effect. Since electron traps have a large effective size (Garlick and 
Gibson, Proc. Roy. Soc. A, 1947, 188, 485), “overlap” of trapping states belonging to 
different centres is likely which may allow electron transfers. This ‘“‘ overlap ”’ of traps can 
also explain the difference in the copper concentrations for optimum afterglow and optimum 
emission efficiency in zinc sulphide phosphors. 
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Science and the Meanings of Truth, by MarTIN JOHNSON. Pp. 179. (London: 
Faber-& Faber, 1946.) “12s. Gd2 net: 


This book continues the line of thought which Dr. Johnson has followed in two earlier 
works, the emphasis here being laid on the meanings of “ truth’’ in modern physics. 
(Somewhat unfortunately the words “ truth” and “ truthfulness”? are used indiscrimin- 
ately : the reader who is accustomed to understanding by the latter the laudable but limited 
quality exhibited by the youthful George Washington must learn to give it a wider meaning 
here.) In the oider physics causal explanations of phenomena in terms of inferred entities 
resembling phenomena were sought through hypotheses, and a hypothesis was true if it 
explained the facts without exception. In modern physics this aim is abandoned and it is 
recognized that physics reveals only a structure characterizing entities about which we 
know nothing. The word “ truth ” is thus left without a meaning, and the author suggests 
that it should be identified with ‘“‘ communicability ’’ and “‘ coherence ”’ as characteristics of 
the structure. The arguments for this view occupy the major part of the book. The 
remainder is concerned with the philosophical problem of the relation between the structure 
and our sense experience and finally with the relation between the scientific and non- 
scientific realms of experience. For reasons of space, only the first part, which most 
directly affects physicists, will be discussed in this review. 

The book reveals, as would be expected, an intimate knowledge of the present situation 
in physics and allied departments of thought. It is not easy reading but it repays the 
necessary effort. The fundamental problems aroused by relativity and quantum theory 
are considered in relation to the modern branch of logic known as “‘ logistics ”’, especially 
the calculus of propositional functions, but strangely enough the post-logistics develop- 
ments, such as those of the metamathematicians and the intuitionists, which transcend 
logistics almost as thoroughly as that transcended scholastic logic and are much more vital 
for the subject treated, are not mentioned. On the question whether the laws of nature are 
empirical or rational in origin, Dr. Johnson’s view seems to be that they are both. As a 
disciple of Milne he thinks that they can be derived from “ our intuitive appreciation of the 
passage of time’, which is taken as empirical, without further observation, by the use of 
reason alone. At least, I think that this is his view. 

The core of the book, however, is the thesis that ‘‘ communicability ” is the dominating 
criterion of truth in physics. Dr. Johnson has emphasized the importance of communi- 
cability in each of his previous works, and it is evidently invested in his mind with profound 
significance. Unfortunately, however, he uses the word in so many senses that it ceases to 
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have any consistent meaning. In one place a thing is communicable if it can be communi- 
cated from one observer to another. ‘That is legitimate, but shortly afterwards it is implied 
that the two observers must agree about what is communicated, and that is not at all implied 
“a the proper use of the word. Yet again, a scientific proposition, to be true, must be 

communicable from one experimental situation to another’, by which is meant that a 
conception, to be true, must be applicable to different, originally independent, fields of 
vestigation. For example, whenever a wave-like property, such as a measurable fre- 
quency or a diffraction pattern, appears in a communicable set of experiments, then this 
particular set achieves communication in turn with other situations which have been 
characterized by wave-like properties. This state of affairs is also called ‘‘ coherence ”’. 
tf this variation of meaning merely caused confusion it would be a remediable nuisance, but 
it is in fact a fatal flaw because a situation, having been endowed with communicability in 
one sense, is thereupon assumed to possess it in others, and so a fallacy is introduced. It is 
easy to prove that black is white if one is allowed to call them both grey. 

The most striking instance occurs in the treatment of the Lorentz transformation. It 
is a fact that this transformation makes different observations (such as those of the orbital 
motion of the Earth and the non-displacement of fringes in the Michelson-Morley experi- 
ment) coherent. Hence it establishes ‘‘ communicability ’? and so is assumed to demand 
the existence of a host of observers to make the communications which enable them to 
agree with one another. All this is simply false. 'The Lorentz transformation has nothing 
to do with different observers ; it concerns the way in which a single observer changes his 
coordinate system. This is not opinion, it is obvious fact. Lorentz did not derive the 
transformation so that he could agree with the Man in the Moon. He derived it so that 
different observations by what was effectively one observer—an earthbound one—could be 
reconciled. Hence, since virtually only one observer is concerned, the question of ‘‘ com- 
municability ” in the proper sense of the word does not arise. ‘The position can perhaps 
be seen more clearly in relation to a simpler change of coordinate systems—one in which we 
merely change the origin of coordinates. The statement that two places are respectively 
10° and 15° east of Greenwich is rejected as a fundamental statement not because it is not 
communicable ; it is communicable and, furthermore, both parties to the communication 
can agree about it, and the statement preserves the coherence of all geography. But it is 
rejected because Greenwich is an arbitrary starting-point and the simpler statement that the 
places differ in longitude by 5° is preferred because it is free from such arbitrariness. ‘That 
has nothing to do with different observers ; it holds good if there is only one observer in the 
world. With the Lorentz transformation the position is exactly the same, with velocities 
substituted for positions. That also has nothing to do with a multitude of observers ; it is 
concerned with the various coordinate systems that a single observer may choose. The 
illusion that different observers, and communicability, are involved is the result of the most 
unfortunate device that early expositors adopted in order to explain the theory of relativity 
to the man in the street. He was assumed not to be familiar with coordinate systems but to 
have a natural tendency to think himself at rest (no man of science has done so since shortly 
after the time of Copernicus). The second coordinate system was therefore personified into 
a second observer who also thought himself at rest. 

All this, I repeat, is not opinion but simple, undeniable fact. As in his earlier works, 
Dr. Johnson overlooks it and bases his argument on the supposed actuality of imaginary 
observers. He is not a dishonest thinker, and one can only conclude that he is so charmed 
by Milne’s theory of kinematical relativity, his devotion to which is manifest throughout his 
writings, that he cannot bring himself to discard the fictitious situation on which Milne 
himself bases it. For if—acknowledging as he does, the Lorentz transformation to be “ the 
classic instance’? of what he calls the “‘ attainment of communicability’’—he takes that 
tcansformation in its legitimate, unperverted form, he will find that he cannot proceed to 
the conclusion he so much desires, and the whole theory of kinematical relativity will be 
exposed as a baseless structure. ‘‘ Communicability ”? makes a host of observers a necessity 
for the possibility of existence, not merely the practicability of prosecution, of science, and 
this is fundamentally false. A single mind, with unlimited intelligence, mechanical skill 
and freedom of movement, and sufficient time, could derive and verify the whole of our 
physical knowledge without anything with which it could “‘ communicate’. If any one 
doubts this he has only to try to name a piece of knowledge not so obtainable. Any 
philosophy which is inconsistent with this fact stands self-condemned, 
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What makes this characteristic of Dr. Johnson’s work the more regrettable is that he is 
one of the few writers who, in this age of specialization, have the breadth of vision and the 
depth of thought necessary to grasp the meaning of the scientific movement as a whole and 
estimate its significance. He has something to say and has the ability to say it. It is 
profoundly to be hoped that he will have the courage to discard fanciful speculations and 
reform his contribution to the philosophy of our time so that it can rest on foundations true 
in the George Washington sense of the word. HERBERT DINGLE. 


Magnetic Materials, by F. Brattsrorp. Pp. ix+156. (London: Methuen 
& Co. Ltd., 1948.) 6s. 


Dr. Brailsford, as a member of the Research Department of the Metropolitan-Vickers 
Electrical Company, has had ample opportunity of first-hand acquaintance with the many 
new magnetic materials which have found, and are still finding, their way into all kinds of 
laboratories and workshops. It is natural that he should open his survey of these materials 
and their properties with a short, but adequate, introduction to modern theories of ferro- 
magnetism, the properties of ferromagnetic single crystals and the various factors which 
determine the magnetic behaviour of metals and their alloys. This introduction, which 
takes up more than half of the booklet, is excellently illustrated with clear diagrams. 

The materials themselves are studied under three headings: Iron and Silicon-Iron 
Alloys, Nickel-Iron and other Alloys, and Permanent Magnet Materials. All students 
should find these sections very helpful. Lists of original references are given, and it should 
therefore be easy to follow up particular points of interest. The author has certainly packed 
a great deal of useful information into a very small space, mainly by the wise use of graphs. 
The book can be unreservedly recommended to all interested in the subject. 

The reviewer was struck by the use of the term retentivity on page 135 to denote residual 
inductance or remanence, for which Ewing used the term retentiveness. As far as can be 
ascertained, the word retentivity has hitherto been used to denote residual intensity of 
magnetization. L. F. BATES. 


Précis de Physique Générale, II: La Chaleur, by P. and A. Mercier. Pp. 105. 
(Neuchatel: Editions du Griffon, 1947.) 12 fr. (Swiss). 


One of the authors is an engineer, former dean of the scientific faculty at the Collége de 
Geneve, and the other is professor of theoretical physics at Berne. Their joint work 
should, therefore, present a good balance between the theoretical and the practical outlooks. 
They have, in fact, written a very good summary of the fundamentals of classical heat theory. 
Classical thermodynamics is presented in a very sound manner, with a very clear treatment 
of the basis of the second law and the definition of absolute temperature, but there is no 
reference to modern quantum mechanics and none to statistical mechanics. Of the subjects 
dealt with, it is interesting to note that whilst internal energy and entropy are introduced, the 
only functions derived from them is Gibbs’ function, here called the thermodynamic potential. 
No mention is made of the function known as heat content or enthalpy, nor of Helmholtz’s 
free energy. 

The book is thus very good for those who desire a brief development of thermodynamics, 
with the practical side of thermometry and calorimetry, but it needs supplementing for 
such matters as heat conduction and radiation. Ole ay, 
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